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Abstract 



Recent proposals of large and infinite extra dimensions triggered a strong research ac- 
tivity in theories in which our universe is considered as a sub-manifold of some higher- 
dimensional space-time, a so-called 3-brane. In this context, it is generally assumed that 
some mechanism is at work which binds Standard Model particles to the 3-brane, an 
effect often referred to as the localization of matter on the brane. Gravity, however, 
is allowed to propagate in general also in the extra dimensions. As demonstrated by 
Randall and Sundrum in 1999, it is also possible to localize gravity itself on a 3-brane. 
In the setup they proposed, the 3-brane is realized as a singular domain wall separating 
two patches of 5-dimensional anti-de-Sitter (AdS 5 ) space-time. The potential between 
two test masses on the brane at distances larger than the AdS 5 -radius was shown to be 
the usual 4-dimensional Newtonian 1/r potential with strongly suppressed corrections. 

The model of Randall and Sundrum, usually referred to as the Randall-Sundrum II 
setup, constitutes the center of interest for this thesis. The main goal of this work is 
to find possible generalizations to higher dimensions of the simple setup considered by 
Randall and Sundrum. One of the motivations for such a generalization is that a realistic 
theory should possibly be able to explain the chiral nature of 4-dimensional fermions on 
the brane. One way to explain chiral fermions from higher dimensions is to consider 
3-branes identified with the cores of topological defects located in a higher- dimensional 
transverse space. Naturally a richer topological structure of the field configuration in 
transverse space provides the possibility of a more realistic spectrum of chiral fermions 
localized on the 3-brane. 

After two introductory chapters on extra dimensions and non-factorizable geome- 
tries which are relevant for the Randall-Sundrum II model, we briefly discuss basics of 
topological defects in the following third chapter. In the rest of the third chapter we 
consider various solutions to higher-dimensional Einstein equations coupled to a series of 
physically different sources and discuss their properties of localization of gravity. Due to 
their asymptotic nature, these solutions are only valid far from the cores of the defects 
in transverse space. Therefore, it seems reasonable to complement the considerations 
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by presenting a particular numerical example of a solution to the Einstein equations 
coupled to a set of scalar and gauge fields: this solution describes a 3-brane realized 
as a 't Hooft-Polyakov monopole residing in the 3-dimensional transverse space of a 
7-dimensional space-time. 

The last chapter of this work is dedicated to the study of a modification of the 
original Randall- Sundrum II model of another kind. The motivation is given by the 
geodesic incompleteness of the latter scenario with respect to time-like and light-like 
geodesies. We will describe a model which resembles the Randall-Sundrum II model 
with respect to its properties of gravity localization but with the advantage that the 
underlying space-time manifold is geodesically complete. Parts of the calculations related 
to the properties of gravity at low energies in this model are rather technical in nature 
and we therefore preferred to assemble them in several appendices. We finally add some 
concluding remarks and discuss possible further developments. 
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Chapter 1 
Introduction 



The discovery of the hidden SU{2) x U(l) gauge symmetry by Sheldon Glashow, Mo- 
hammed Abdus Salam and Steven Weinberg which provides a unified description of the 
physics of electromagnetic and weak interactions between elementary particles is com- 
monly considered not only as one of the major contributions to our understanding of 
the physics of elementary particles but also as an outstanding achievement in theoretical 
physics in the 20^ century as a whole. Together with the theory of the strong inter- 
actions of elementary particles, quantum chromodynamics, based on the gauge group 
SU(3), the Glashow-Salam- Weinberg theory constitutes the so-called Standard Model of 
the interactions of elementary particles, the pride of particle physicists today due to its 
impressive agreement between observations and predictions. 

Despite of the enormous successes of the Standard Model, it is Weinberg himself who 
starts his by now famous review article Ij of the so-called cosmological constant problem 
with the lines: 

Physics thrives on crisis. 

The crisis Weinberg describes is not only a result of the puzzling and yet unexplained 
mismatch between the observed and predicted values of the cosmological constant of 
around 120 orders of magnitude, but also stands for numerous unsuccessful attempts to 
elucidate the quantum structure of the fourth fundamental interaction in nature, gravity. 

Since its introduction by Einstein in 1916, the classical theory of general relativity 
remains the main theoretical framework for our understanding of the expanding universe 
and its history. On the other hand, general relativity so far withstood all attempts to 
be combined with the principles of quantum mechanics, hence we presently still lack a 
successful theory of quantum gravity. 
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From today's point of view, string theory appears to be the most promising candidate 
for achieving the goal of a unified description of all fundamental interactions since the 
graviton is part of the field contents of all known theories of strings. A consistent 
formulation of string theory requires that the dimensionality of space-time is equal to 10 
(or 11 in the case of M-theory). The reduction of a 10-dimensional string theory to the 
4-dimensional Standard Model in the limit of low energies is complicated by the wealth 
of possibilities of compactifying the additional dimensions. 

Even though the requirement of a 10-dimensional space-time in string theory was 
unexpected, the idea of extra dimensions was not new. It was first proposed by Nord- 
strom j2] and further developed by Kaluza and Klein E] . 

Later, inspired by the success of Kaluza in the unification (on a classical level) of the 
gravitational and electromagnetic interactions (the only two known interactions of his 
times) many physicists applied the idea to non-abelian gauge theories, supergravity and 
string theories in the search for a realistic Kaluza-Klein theory jSJ [HI 13 El El Ell HH E2] ■ 

The idea of extra dimensions also fell on fruitful grounds outside the usual Kaluza- 
Klein approach. During the last few decades, theoretical physics discovered the potential 
of extra dimensions in addressing many outstanding problems of physics: the cosmolog- 
ical constant problem [121 EH EH], the origin of chirality (THl 113 ESI EH1 12U| - the strong 
CP-problem |2*T) 12*2*] , the hierarchy problem [23J EH EH EI3 EZ] , and the horizon-problem 
[2H] to mention only some of the most prominent ones. 

In spite of the long history of extra dimensions, a boost in research in this area was 
triggered only recently. In the search of a solution to the hierarchy problem it was realized 
that, from a phenomenological point of view, the existence of large extra dimensions from 
about 1mm to lfm is not excluded by observations [2*41 I25j . as long as one assumes a 
mechanism of localization of the Standard Model fields on a 4-dimensional sub-manifold 
in the extra dimensions. The reason is that the validity of Newton's law of gravity 
is established only at distances larger than about 0.1mm j2Hl- Although models with 
extra dimensions in the above distance range suffer generically from the effects of light 
Kaluza-Klein gravitons, the corresponding constraints from particle physics [HO]- [US] and 
astrophysics and cosmology [HI)] -[HE] are not stringent enough to rule them out. 

Moreover, inspired by the discovery of solitonic solutions in string theory [39J, so- 
called D(irichlet)-branes, it was pointed out that even an infinite extra dimension can be 
in complete agreement with 4-dimensional observations [40] . In particular it was shown 
that under the assumption that Standard Model fields were bound to a so-called 3-brane 
in an 5-dimensional bulk space-time characterized by a non-fact or izable geometry, it is 
possible to reproduce standard 4-dimensional Newtonian gravity on the brane. Even- 
tually, a similar setup involving two branes using the same idea of a warped geometry 
provided a solution to the hierarchy problem different from the one mentioned above [2H] . 
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The main advantage of this latter approach is that it does not predict light Kaluza-Klein 
excitations of the graviton. 

The fact that these new ideas of large and infinite extra dimensions make use of the 
concept of a 4-dimensional sub-manifold, a 3-brane, embedded in a higher- dimensional 
space-time resulted in an immense effort in the study of properties of models of this kind. 
The investigations in this area of brane- world scenarios range from questions of stability 
jlU H21 EH] , general space-time structure jJHEHl, properties of gravity |4H Ho] W7\ US], 
black holes on the brane cosmological aspects (related to a brane motion in the 

bulk) [53]- [7H] including brane inflation [IH]-[E2], to gauge invariant formalism to study 
brane perturbations (HH1 EH EH] , etc. 

This work constitutes a contribution to the above research activity related to gener- 
alizations of the above mentioned Randall-Sundrum II setup [40] to either higher dimen- 
sions or to models with different spatial topology. 

The structure of this work is the following: in chapter|2]we review the subject of non- 
factorizable geometries and in particular the example of the Randall-Sundrum II setup 
40 which lies in many respects at the center of our considerations. The developments 
of the following chapters will be in one or another way related to the above mentioned 
setup. 

The subject of chapter El illustrates another example of a fruitful symbiosis of the 
a priori unrelated subjects of topological defects and extra dimensional physics. The 
physics of topological defects (or topological solitons) is extremely rich and it is clearly 
out of the scope of this report to provide a complete treatment. We refer here to excellent 
textbooks, relevant articles and reviews [HB] - P"D4j and references therein since our main 
goal will be to provide basics ideas and concepts that will be of relevance for this work. 
Based on the idea of identifying a 3-brane with a topological defect residing in the 
transverse space of some higher-dimensional space-time [Tfi] . we present a discussion of 
solutions to higher- dimensional Einstein equations with different kinds of source terms. 
Attention will be paid in particular to the localization of gravity in models of this kind, 
where by localization we mean that (at least) the effective 4-dimensional Newtonian 
theory of gravity is reproduced on the brane. The solutions for the space-time metric 
presented in this chapter will be of an asymptotic nature with the meaning that they will 
constitute solutions to the Einstein equations only for large distances from the topological 
defects. Under certain conditions the underlying field theory allows to consider a singular 
defect limit in which the size of the defect tends to zero in such a way that the above 
asymptotic solutions become complete solutions when supplemented by appropriate fine- 
tuning relations. In order to support the analytical evidence for the existence of a 
particular set of these solutions (involving generic p-form fields in the bulk) we present 
a specific numerical example in chapter a 't Hooft-Polyakov monopole residing in the 
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3- dimensional transverse space of a 7- dimensional space-time. After a short discussion of 
't Hooft-Polyakov monopoles in Minkowski space-time, we present numerical solutions of 
the 7-dimensional case of interest. Moreover, we identify the regions in parameter space 
leading to the localization of gravity together with the corresponding fine-tuning relations 
linking various brane tension components. Finally, we present a specific 5-dimensional 
brane-world model in chapter with similar phenomenological properties (in particular 
similar large distance corrections to Newton's law) as the Randall- Sundrum II setup. 
The main difference to the latter is that the global topology of the spatial sections of its 
metric is compact. In contrast to the Randall-Sundrum II model in which the 3-brane 
is characterized by a flat, induced metric, the induced metric of the model we present is 
given by the metric of Einstein's static universe. The main advantage of this modification 
is that the given setup does not suffer from being geodesically incomplete, a problem of 
the Randall-Sundrum II model. We close by drawing conclusions in chapter 15.51 



This thesis is based on the following publications: 

T. Gherghetta, E. Roessl and M. E. Shaposhnikov, "Living inside a hedgehog: Higher- 
dimensional solutions that localize gravity", Phys. Lett. B 491, 353 (2000) [arXiv:hep- 
th/0006251]. 

E. Roessl and M. Shaposhnikov, "Localizing gravity on a 't Hooft-Polyakov monopole in 
seven dimensions", Phys. Rev. D 66, 084008 (2002) [arXiv:hep-th/0205320]. 

A. Gruppuso, E. Roessl and M. Shaposhnikov, "Einstein static universe as a brane in 
extra dimensions", [arXiv:hep-th/0407234]; accepted for publication in JHEP. 
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Chapter 2 

Non-factorizable geometries 



2.1 Introduction 

In this introductory chapter we will discuss briefly the relevance of non-factorizable ge- 
ometries to brane-world scenarios. After collecting basic formulae for the metric, the 
connection coefficients, the Ricci tensor and the Ricci scalar of non-factorizable geome- 
tries and presenting some motivations for their study, we will briefly discuss the Einstein 
equations with a cosmological constant term in vacuum. Then we will study the two 
examples of the Randall-Sundrum brane-world models. Our motives for the last choice 
are two-fold: first, as we already emphasized, it is to a good part the Randall-Sundrum 
models which renewed and reinforced the interest in non-factorizable geometries. Second, 
because it will serve as a reference for later setups studied in chapters 0HS1 in appendix 
El we consider a simple generalization of the metric ()2.1|) given in section 12721 below. 



2.2 Metric, connection coefficients and curvature of 
non-factorizable geometries 

The concept of non-factorizable product manifolds was first introduced by Bishop and 
O'Neill in 1969 [105J and later applied to general relativity and semi-Riemannian geom- 
etry |lU6t I107j . Many of the well known metrics of general relativity have this form. 
For example the simple Schwarzschild solution of general relativity in four dimensions 
can be considered as a non-factorizable product of a two-sphere S 2 and a 2-dimensional 
space-time represented by the Carter-Penrose diagram for a black hole [108 . For recent 
work in this area we refer to 109, 1110[ llllj . We will point out later the relevance of 
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non-factorizable geometries in theories with extra dimensions. 

By a non-factorizable metric we mean a metric which locally can be written in the 
following form: 1 



ds 2 = cjmn dx M dx N = a (x a ) g^ (x p ) dx p dx v + g ab (x c ) dx a dx" . (2.1) 

The function a (x a ) is supposed to depend only on the "internal" coordinates x a and is 
sometimes called a "warp" factor. Moreover, geometries described by metrics of the kind 
()2.1j) are often referred to as warped geometries. We will assume that the dimensionality 
of space-time is given by d\ + g?2, where d\ and d<i denote the dimensionalities of the 
spaces described by g^ u (x p ) and g a b (x c ), respectively. From the explicit expression (|2.1|) 
it is straightforward to derive the following relations for the connections Tm N , T c ab 

K b = K b , f^ = -i^vv, f; b = ^v b a, f^ = r^, (2.2) 



Ricci tensors Rmn, Rfiu, Rab 



RflV RfiV r\9lLV 



I a 



(2.3) 



^1 .=, ^ dl ~ 

2 a 4 cH 

i^a = 0, (2.5) 



R a b = Rab-^-V a V b a + T ±V a aV b (T, (2.4) 



and the Ricci scalars R, R and R corresponding to the metrics c/mn, g^v and g ab - 

R = R + -R - ^ V a V V - rfl( ^~ 3) V a ,VV. (2.6) 

In the upper formulae V a denotes the covariant derivative associated with the metric 
<jab {% c )- Note that in the case a = 1 we recover in eq. ()2.6)1 the well-known property of 
additive Ricci scalars for direct product spaces. 



1 Throughout this report we will consistently use the sign conventions of ref. [112J for the metric and 
the Riemann and Ricci tensors. The metric signature is mostly plus — h • • • +■ The expression for the 
Riemann curvature tensor is Rp 7 s — ^BSy ~ ^p-y 8 + ^"7-^/35 — ^aS^0-y an d the Ricci tensor is defined 
by R af} = R° aafi . 
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2.3 Einstein equations with a cosmological constant 
term 



One of the early motivations of considering a metric of the form (j2.1|) was to provide a 
solution to the cosmological constant problem. The original hope was that the higher 
dimensional Einstein equations with a cosmological term might give a dynamical ex- 
planation for the smallness of the observed 4-dimensional cosmological constant. The 
d\ + d.2 dimensional Einstein equations with a cosmological constant and an otherwise 
empty space-time are given by: 

Rmn — 2 9mn R + A g MN = . (2.7) 
Using the above relations ()2.3|) - (j2.fij) in ()2.7|) gives after some algebra: 

Rnv — 2 9tiv R + Aphy S 9tiv — ) (2-8) 
R»b = ^-A o A 9ab + d! ( —V a V b a- — V a aV b a J , (2.9) 

These equations were first studied in [TnjIIlllISllHj, in order to provide an alternative way 
to the usual Kaluza-Klein compactification of higher dimensions and, as we mentioned 
above, a way to address the cosmological constant problem. In ^3] solutions of the 
above equations (l2~%J) - (l2~TU|) are given in the case of two extra dimensions (d 2 = 2). 
Since Ap^y S can be considered as an arbitrary integration constant, the solution of the 
cosmological constant problem relies on an additional explanation why a vanishing Ap^y S 
should be preferred. In ^5] it is demonstrated that a similar continuous spectrum of 
solutions exists in the case of more than two internal non-compact dimensions. 

In the next two sections we discuss the two brane-world models of Randall and 
Sundrum [2S],|ini as well known examples of warped geometries in physics of extra 
dimensions. 



2.4 The Randall- Sundrum I model 

In 1999 Randall and Sundrum proposed a solution to the hierarchy problem based on a 
single spatial extra dimension (251- In their model, the observable 4-dimensional universe 
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is considered to be a subspace of a 5-dimensional space-time, referred to as a 3-brane. 
The metric that characterizes the model takes the following form: 

ds 2 = cjmn dx M dx N = e - 2krc(t> T}^ dx^ dx v + r 2 d<j) 2 . (2.11) 

Upper-case latin indices M , N range from to 4. Greek indices denote the coordinates 
of the four familiar dimensions and range from to 3 while <j) G [0, tt] denotes the 
extra (compact) dimension (of "size" irr c ) given by the quotient S 1 /Z 2 . T]^ is the usual 

Minkowski-metric in 4-dimensions with signature h ++. The line element ()2.11|) 

describes a certain region of AdS^-space-time and k can be identified with the inverse 
AdS-r&dius. There are actually two 3-branes in this setup located at the boundaries of 
the interval [0, tt]: a negative tension brane located at <fi = it with the usual standard 
model fields localized on it and a positive tension brane at = needed for overall 
consistency. The action functional can be written in the following way: 

TT 

S = J d 4 x J d<p y/=g (-A + 2 M 3 R) (2.12) 

— TT 

+ J d 4 x^-g vis (£ vis - cr vis ) + J d A xy/-g hiA (C Md - a hid ) 

In the last expression A and M denote the 5-dimensional cosmological constant in the 
bulk space-time and the fundamental scale of gravity. g vis (g hid ), £ vis (C Md ) and cr vis (a hid ) 
represent the induced metric, Lagrangian and tension of the visible (hidden) brane. 



As shown in detail in [2^], the metric (jZHTj) solves the Einstein equations following 
from ()2.12j) (in the sense of distributions) provided the following relations between the 
brane-tensions, the cosmological constant, the fundamental scale and the inverse AdS- 
radius hold: 



o"hid = — o"vi S = 24 M k , K = -2AM 6 k\ (2.13) 

From the action (J2.12j) it is straightforward to relate the fundamental scale M to the 
4-dimensional Planck-scale M 2 X with the result: 

M 3 

M 2 1 = — (l-e- 2kr ^) , (2.14) 

K 

which shows that M P1 stays finite in the limit r c — > 00 contrary to the case of a product- 
space compactification. 

The way the model solves the hierarchy problem can be summarized as follows: 
physical mass parameters as measured by an observer on the visible brane are generically 
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related to the fundamental mass parameters via a renormalization of the induced metric 
g vis . For kr c rs 50 the exponential suppression of the metric at the position of the 
visible brane implies that any fundamental mass parameter of the order of the Planck 
mass corresponds to a physical mass in the TeV regime from the point of view of an 
observer on the visible brane. The predictions of the model are very different from the 
predictions of models with large extra dimensions introduced in (2HI2S1- The constraints 
coming from particle physics [3U]-|III3; astrophysics and cosmology [SS]-|SH] restricting 
large extra dimensions do not apply in the case of the Randall-Sundrum I model, since 
the splittings of the Kaluza-Klein spectrum of gravitons are roughly given by the weak 
scale. Moreover, their coupling to matter is also expected to be set by the weak scale 
and not by the Planck scale. A somewhat unpleasant feature of this model is that the 
fields of the Standard model are localized on a 3-brane with negative brane tension. 

2.5 The Randall-Sundrum II model 

In a related article [40], Randall and Sundrum demonstrated that even an infinite extra 
dimension is not in contradiction with present day observations of gravity at low energy 
scales. The main motivation was to prove that the Newtonian force law is effectively 
4-dimensional for two test masses placed on the 3-brane. In order to achieve this "local- 
ization of gravity" , it is necessary to assume that the Standard Model fields are localized 
on the brane at <ft — (in the coordinates of the last section). The brane at <fi = it serves 
merely as a "regulator" brane for imposing boundary conditions and will eventually be 
removed from the setup by taking the limit r c —>■ oo. 

2.5.1 Localization of gravity 

In order to work out the force between two test masses on the brane we need to determine 
the spectrum of Kaluza-Klein excitations of the 4-dimensional graviton h^ u defined by: 

ds 2 = (e- 2 k M 7? M1 , + V) dx" dx v + dy 2 , (2. 15) 

where y = r c (p in the notation of the last section. 2 Following [3U| we use the gauge 
d^h^ = h 1 ^ = in order to obtain a single perturbation equation for the metric fluc- 
tuations h^ u from the linearized Einstein equations in the background (J2.ll)) . After 

2 Note that even thought in reference (3^ Randall and Sundrum stick to the topology S 1 /Z2 at the 
beginning, they allow y to range from (— oo, oo) later in the paper, in order to avoid a semi-inhnitc extra 
dimension. 
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neglecting the tensor structure of an expansion in 4-dimensional plane waves gives: 



-^e 2k ^- 1 -d 2 -2k5(y) + 2k 2 



i/>(y) = 0. 



(2.16) 



We now transform the above equation into a Schrodinger equation by using the change of 
independent and dependent variables z = sign(y) (e fc l y l — l) /k, i>(z) = ip(y) e fc ' y '/ 2 with 
the result: 

l„o 15 k 2 3 k 



--di + * - — 8(z) 

2 2 8 (A;|^| + 1) 2 2 1 ; 



= m ip(z) . 



(2.17) 



From the potential term in the above equation one can infer the main properties of 
the Kaluza-Klein spectrum in the limit r c — > oo. The 5-function allows for a single 
normalizable bound state $o(z) = k^ 1 (k\z\ + 1)~ 3 ^ 2 corresponding to m 2 = 0. Next, 
since the potential vanishes in the limit z — > oo there will be a continuous spectrum for 
all m 2 > 0. Finally, one expects a suppression of the wave-functions of the continuum 
in the vicinity of the brane (z = 0) due to the potential barrier. The correct continuum 
wave-functions are given in terms of Bessel functions: 



A + 



k 



AYo 



m 



z + 



k 



+ BJ 2 



m 



z + 



k 



(2.18) 



The constants A and B can be determined by going one step back and using the regulator 
brane to impose the following boundary condition on ip(z) at z = and z = z c = 



d z ip(z) 



3k 



Mz) 



0. 



(2.19) 



2(kz + l) 

In this way the spectrum becomes discrete again and at least in the most interesting 
case of small masses m k the correct normalization of the modes can be deduced. 
As usual, masses in the range m k will determine the behavior of the potential in 
the distance regime r ^> 1/k. Going back to the infinite extra dimensions, we can now 
assemble all contributions coming from the zero-mode and the continuum modes in the 
potential between two test particles on the brane at z = 0: 



V(r) 



m 1 m 2 
(jjv r 



dm 



G 



N mi m 2 e mr m 



k 



= G 



mi m 2 



N~ 



r 2 A; 2 



(2.20) 
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with Gn = k/M 3 following directly from (|2.14j) . We see that it is the zero mode that 
reproduces the conventional 4-dimensional Newtonian behavior while the collective con- 
tributions of the higher modes is manifest through a correction term which is strongly 
suppressed at distances r larger than 1/k. 

We conclude this section with a few remarks about the above calculation. First, 
the result (|2.2(J|) applies only in the regime r S> 1/k, since the normalization constant 
has been worked out correctly only in the regime of small masses. Second, due to the 
incomplete treatment of the tensor structure of the gravitational perturbation we expect 
()2.20|) to represent the correct potential only up to numerical factors in front of the 
l/(r 2 k 2 ) -correction term inside the parenthesis on the right hand side. Indeed several 
careful studies have shown that the factor of 1 in (J2.2U|) has to be replaced with 1/2 in the 
case of the massless scalar in the Randall-Sundrum background, whereas a full treatment 
accounting for the tensorial nature of the graviton and the effects of "brane-bending" 
gives a factor of 2/3, see for example [4^ I48| Hi 31 I114j and |115j and references therein 
for a recent discussion. We also refer to appendix IU1 where we redo the above calculation 
bypassing the complications of the determination of the Kaluza-Klein spectrum and 
wave-function normalization by directly solving for the (scalar) Green's function. 

The singular treatment of branes in the models of Randall and Sundrum, manifest 
through the appearance of delta-functions in the stress-energy tensor components, con- 
stitutes a welcome simplification when it comes to solving Einstein equations. However, 
from a field theoretical point of view, branes can be considered as finite energy solutions 
to some field equations with an energy- density localized in some region of space of a 
given finite thickness 6. On the one hand, in the case of a singular brane, a complete 
solution to Einsteins equations can be obtained in a straightforward way by imposing 
the Israel-junction condition |116j at the boundary between the two regions separated 
by the brane. On the other hand, the treatment of the non-singular, thick-brane case 
is in principle based on the complete solution of the coupled Einstein and matter field 
equations and is therefore in general a much more involved task with the result that 
analytic solutions are rare. As we will see in the following two chapters, the idea of re- 
placing a singular brane by a thick domain wall can be generalized to higher- dimensional 
cases. Depending on the dimension of space-time, we will have to replace domain walls 
by various other topological defects like strings or monopoles. 
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Chapter 3 



Topological defect solutions of 
higher-dimensional Einstein 
equations 



3.1 Introduction 



Topological defects like domain walls, strings and monopoles arise in a variety of different 
areas of physics starting from condensed matter physics, where particular examples of 
defects can be observed experimentally, over particle physics, astrophysics and cosmology 
to even cosmology involving extra dimensions. Their physics is enormously rich and 
consequently the attention devoted to the study of topological defects in the past and 
present is considerable. Before embarking on a discussion of the relevance of topological 
defects in theories with extra dimensions we would like to review briefly basic properties 
and examples in the next section. It is clear that our presentation is again not more 
than an introduction to this vast field and in no sense complete. We start by giving a 
definition for what we will mean by topological defect (or topological solitons), enumerate 
the most prominent examples and discuss general common physical properties. After a 
short sketch of homotopy theory, together with the statement of the two fundamental 
theorems, we are ready to discuss the by now standard classification scheme of topological 
defects going back to Kibble [HE] • We will conclude this introduction by shortly discussing 
the role of particular topological defects in condensed matter physics, astrophysics and 
cosmology. 
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3.1.1 Definitions and physical properties of topological defects 



Topological defects or equivalently topological solitons are particular cases of a larger 
group of objects called solitons or better solitary waves. The last two expressions are not 
used in a standardized way in the literature. We decide to use the definitions given in 
[H7] . On the one hand, a solitary wave is a classical solution to any in general non-linear 
field equation (or sets of non-linear field equations) with the properties that 

1. the total energy of the given field configuration is finite, 

2. the energy density e(t, x) is nonsingular and localized in space. Moreover, for 
time-dependent solutions the shape of e(t,x) propagates undistorted in time. 

A soliton, on the other hand, is a solitary wave with the additional property that the 
underlying nonlinear dynamics involved in a scattering process of two or more solitary 
wave packets asymptotically (for large times) restores the forms of the incoming wave 
packets. Obviously, this definition is very restrictive and apart from a few examples the 
vast majority of known solitary waves does not satisfy this criterion. 1 With the above 
definitions, we can now come back to topological defects and identify them with solitary 
waves or solitons whose stability is guaranteed by means of a topological argument. More 
precisely this means that the underlying field theory allows for a topologically conserved 
(and generically quantized) charge in such a way that field configurations with a given 
charge Q are forbidden to evolve into configurations characterized by charges different 
from Q. 

As we will see in the next section, a mathematical classification of topological defects 
is based on specific topological properties of the vacuum manifold of the underlying field 
theory. A more phenomenological (and related) approach classifies defects according to 
their spatial dimensionality. Therefore, one distinguishes 

• O-dimensional, point-like defects or monopoles, 

• 1-dimensional, line-like defects or strings, 

• 2-dimensional, surface-like defects or domain walls. 

Somewhat out of this scheme lie textures which are not localized in space and hybrids of 
the above mentioned defects. For example one can consider theories in which monopoles 
appear to be attached to the endpoints of strings or in which domain walls are bounded 
by strings. 

1 It is mainly for this reason that it became common practice to use the term soliton instead of solitary 
wave even in cases where the above criterion is not satisfied. 
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Some observed physical examples in condensed matter physics comprise magnetic 
flux lines in a type II superconductor, quantized vortex lines in superfluid 4 He, textures 
in 3 He, line defects and dislocations in crystalline substances, various defects in different 
phases of liquid crystals and domain structures in ferromagnets. The formation of all 
kinds of defects is believed to be triggered via rapid phase transitions in connection with 
spontaneous symmetry breakdown by the so-called Kibble-mechanism [86, • The main 
idea is that below a certain critical temperature the effective Higgs potential will develop 
degenerate minima. Random quantum fluctuations of the Higgs field will decide which 
vacuum will be chosen in a given volume. If two neighboring regions in space turn out to 
fall into different minima, then a defect will form at the boundary between them. The 
typical length scales of these structure will be set by the correlation length beyond which 
the Higgs phase is completely uncorrelated. An example for the above situation is the 
magnetic domain structure in a ferromagnet. 



3.1.2 Homotopy theory and classification of topological defects 

In this section we intend to assemble basic definitions and results of homotopy theory as 
far as they are of relevance to the classification of topological defects. Our presentation 
is based on [88| . For more complete discussions of the subject we refer to [EH1 HUB EH E2] ■ 
We start with the definition of homotopy groups. Let us consider the set C x (Ai) of all 
closed loops with given base point x G Ai in a given topological space Ai, that is the 
set of all continuous maps / from the circle S 1 into Ai with x in the image of /. We 
can define an equivalence relation on Ai in the following way: / and g G C x (Ai) are 
said to be homotopic if they can be continuously deformed into each other (within Ai) 
keeping the base point x fixed. Now we have / ~ g if / and g are homotopic and we 
denote [/] the equivalence class or homotopy class corresponding to the loop /. This 
set of equivalence classes {[/], [<?],■■ ■} can be given a group structure: the inverse of [/] 
is defined by [f]^ 1 = [f^ 1 ] with f~ x being the loop obtained from / by reversing the 
order in which it is traversed. The group multiplication of [/] and [g] is realized by the 
operation of connecting two loops according to [/] ■ [g] = [fog] and the identity class 
[I] is the class of loops continuously contractible to the point x. The resulting group is 
called the first fundamental group or the first homotopy group of the topological space 
Ai with base point x and is denoted by tti {Ai ,x). It can easily be seen that for a 
connected space Ai two base points x,y G Ai give rise to isomorphic first homotopy 
groups 7v i (Ai ,x) = n 1 (Ai ,y) so that it is justified to drop the reference to the base 
point altogether and simply write m (Ai). 

In a similar manner it is possible to define homotopy groups of Ai of higher order as 
equivalent classes of maps from n— spheres into Ai. For example the n homotopy group 
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of the topological space At, denoted by 7t n {At) is the set of homotopically equivalent 
classes of maps form the n-sphere into the manifold At. 2 An important difference between 
the groups 7Ti {At) and i\ n {At) (for n > 2) is that the former is not necessary abelian 
while the latter always are. Finally, it is useful to introduce 7r {At) as the homotopy 
classes of maps from S° (a point) into At. This set labels the disconnected components 
of the space At. However, in general it cannot be given a group structure. A trivial 
7r {At) is equivalent to the statement that At is a connected space. 

We are now prepared to state (without proof) the first two fundamental theorems of 
homotopy theory. 

First Fundamental Theorem 

Let G be a connected and simply-connected Lie group {n {G) = i\\ {G) = I), and let if 
be a subgroup of G. We then have 



Second Fundamental Theorem 

Let G be a connected and simply-connected Lie group {n {G) = it\ {G) = I), and let if 
be a subgroup of G. Let if be the component of if connected to the identity. We then 
have 



The usefulness of the above two theorems follows directly from the classification 
schema for topological defects which we set out to describe in the following. Topological 
defects are thought to be produced during phase transitions in the early universe. Let us 
suppose that in the high temperature phase the theory possesses a symmetry described 
by a connected, simply-connected continuous Lie group G. At the onset of the phase 
transition, when the temperature drops below a critical value, the Higgs potential devel- 
ops non-trivial minima and the symmetry is spontaneously broken down to the subgroup 
if of G. Let At be the vacuum manifold of the theory and $ a given vacuum expectation 
value in At. Then D{h)&o = $o f° r all /i e if by definition of the unbroken subgroup 
if. On the other hand, any g G G with g ^ if will give rise to another possible choice of 
vacuum distinct from $o- More than that, all different distinct elements in the vacuum 
manifold At can be obtained by the action of distinct cosets g if of G with respect to 
the unbroken subgroup if on <3> . Therefore, the vacuum manifold can be identified with 
the coset space At = G/H. 

2 Note that we immediately dropped the dependence on the base point which again is possible only 
for a connected space. 



TTi {G/H) = 7T (ff ) . 



(3.1) 



7T 2 {G/H) = TTi (ff ) • 



(3.2) 
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We are now in a position to state the classification of topological defects established 
by Kibble in 1976 [86J. Accordingly, domain walls occur during phases transitions ac- 
companied by spontaneous symmetry breaking if the vacuum manifold has disconnected 
components, that is ttq (G/H) ^ I. Strings occur if the vacuum manifold contains un- 
shrinkable loops that is 7i\ (G/H) ^ I. Analogously we find monopoles for ix?, (G/H) ^ I 
that is unshrinkable 2-spheres. A somewhat special kind of defects that can also be in- 
corporated in the above classification scheme is so-called textures. While the Higgs field 
of domain walls, strings and monopoles leaves the vacuum manifold inside and around 
the core of the defect, the Higgs field for textures always takes values inside the vacuum 
manifold. Therefore, only gradient energy contributes to the total energy of textures. 
Textures arise in theories in which the vacuum manifold allows for non-shrinkable 3- 
surfaces that is for non-trivial 713 (G/H). We summarize the classification of topological 
defects according to the topology of the vacuum manifold of the underlying field theory 
in the following table 13.11 



Topological Defect 


Classification 


Domain walls 
Strings 
Monopoles 
Textures 


7T0 (M) 
TTl (M) 
7l 2 (M) 
7T 3 (M) 



Table 3.1: Homotopy classification of topological defects. A4 = G/H. 



By making use of the first fundamental theorem ()3.1|) . we see that for connected and 
simply-connected groups G the question of whether a given theory allows for topological 
string solutions can be reduced to the question of whether the unbroken subgroup H has 
disconnected components or not. In a similar manner (and under the same assumptions 
about G) the condition for the existence of monopole solutions 7r2 (G/H) ^ I is equivalent 
to 7Ti (H ) 7^ J, where we recall that H was the component of H connected to the identity. 

Since the discoveries of topological defect solutions around 30 years ago in simple 
relativistic field theoretical models [HHJ EH EE1 ESI EH EEj, the field underwent an immense 
boost in research activities in different directions. Apart from their direct relevance in 
condensed matter physics where corresponding non-relativistic analogs of the above- 
mentioned models describe magnetic flux lines and vortex lines in superconductors and 
superfluids, the rich potential for new physics was immediately recognized also in the 
areas of particle physics, astrophysics and cosmology. Research interests ranged from field 
theoretical, classical and quantum aspects of topological defects over revealing detailed 
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defect structures, extension of the classes of known defects to more and more complicated 
field theoretical models including supersymmetry versions, dynamics of interactions of 
defects accompanied by particle productions, the study of hybrid defects etc. Also the 
gravitational properties of defects were extensively studied for example the gravitational 
lensing effects of strings, emission of gravitational waves from oscillating string loops, 
formation of black holes, or the gravitational fields of collapsing textures to mention 
only a few. 

Interest in cosmology ranged from the formation and evolution of domain structures or 
string networks in the universe in connection with structure formation, predicted power 
spectra, imprints in the temperature fluctuation of the cosmic microwave background 
radiation, etc. Our list is by no means complete and we refer to excellent textbooks and 
review articles for more exhaustive descriptions [HI EH IH1 CDOl DHH CHS QUI EH and 
also for lists of relevant articles. 

We now turn to our main interest for studying topological defects, namely to point 
out the role they play in cosmological models in connections with extra dimensions. As 
we will see in the next section, there are numerous motivations for this symbiosis. 

3.2 Topological defects in higher dimensions 

The picture of identifying our 4-dimensional universe with a topological defect in a higher 
dimensional space-time has emerged in an attempt to find alternatives to Kaluza-Klein 
compactifications [HH H31 EH [231 12U [2E1 IHH EH EEE1 EH- In this connection, 
the 4-dimensional sub-manifold of the higher- dimensional bulk associated with our world 
is usually referred to as a 3-brane. While in theories of this type, Standard Model 
particles are thought to be confined to the brane, gravity is free to propagate in the 
extra dimensions, being the dynamics of space-time itself. One way of achieving such 
a confinement is to use string theoretical arguments [39|. Other ways use various field 
theoretical effects. In 16J it was shown that in the absence of gravity a domain wall 
residing in a 5-dimensional Minkowski space-time can explain the necessary trapping 
of Standard Model particles in the vicinity of the brane, making the extra dimension 
inaccessible to particles below a certain threshold energy. Moreover, it was shown that 
this mechanism can give rise to 4-dimensional chiral fermions on the brane. Similar 
mechanisms of confining particles to 4-dimensions had also been envisaged in theories of 
pure gravity |118j . In |119j a possibility of obtaining massless gauge bosons confined to 
domain walls has been demonstrated. Lattice simulations on related topics can be found 

in dmrm]. 

As already discussed in chapter |2l also gravity itself can be localized on the brane in 
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a simple (singular) domain wall setup in 5 dimensions involving a warped geometry [UJ| > 
even in the case of an extra dimension of infinite size. This last work motivated also the 
generalization of the simple model of Randall and Sundrum to more involved possibilities 
and to higher-dimensional models. In |122j a non-singular supergravity version of the 
Randall- Sundrum II model was considered. Motivated by solving the hierarchy problem 
along similar lines as in the Randall-Sundrum I model and based on their earlier works 
in 4-dimensions |123| , the authors of j2Z] considered a brane made from a global cosmic 
string defect residing in a six- dimensional space-time. The solution corresponds to a 
vanishing cosmological constant in the bulk and is characterized by a singularity at 
a finite distance from the string. A similar but singularity-free static solution with a 
negative bulk cosmological constant based on the compactification on a global string 
was found in |124j . More general constructions in six dimensions have been considered 
in |125j . Then, in |126j it was shown that similar to the Randall-Sundrum II case in 
5-dimensions, gravity can be localized on a global string like defect in six dimensions. 
The corrections to Newton's law were also calculated and shown to be highly suppressed 
at large distances. See also |127j for a related discussion of the physical properties of 
the thin-string limit and j!28j for a more detailed study of the localization of gravity 
in this model. A discussion of the localization of matter fields of different spin in the 
above background of the local string can be found in |129j . Further generalizations of 
these scenarios to so-called "hedgehog" scalar and p-from field configurations were given 
in [130j. p^T| and |132| . In the latter paper, the corrections to Newton's law were given 
for theories with n-transverse extra dimensions. After finishing this short overview of 
existing literature, we will come back to |132| in the rest of this section in more detail. 
While |13Uj and |132j concentrate on the case of a spherical symmetric transverse spaces, 
in |133j more general cases of Ricci-flat transverse manifolds or homogeneous spaces 
with topologically non-trivial gauge field configurations were included. The motivation 
for considering more complex structures of transverse spaces lies within the possibilities 
of obtaining a more realistic spectrum of 4-dimensional particles localized on the brane, 
see e.g. |134[ I135j for a construction involving three families of chiral fermions localized 
on the brane which in this setup is identified with an Abrikosov-Nielsen-Olesen vortex 
in 6 dimensions. Another attempt to construct a realistic theory of extra dimensions 
in the Kaluza-Klein approach can be found in 136 , with the interesting observation 
that the identification of extra components of a gauge field with the Higgs field provide 
a possible realization of the spontaneous electroweak symmetry breaking. In the latter 
model, the existence of a topological background in the extra dimensions is crucial for 
explaining both, the symmetry breaking of the Standard Model and chiral fermions in 
4 dimensions. More recently (see |137j and references therein) similar work realized a 
purely gravitational localization of three generations of fermions on a brane. 
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As already mentioned in the introductory chapter, it was realized early [TBJ H7J CHI 
UJU I2IH that the topological nature of the field configuration in transverse space is crucial 
for obtaining chiral fermions in 4-dimensions. For a recent discussion of the localization 
of chiral fermions in brane-world models see |138l 139 . 

We should also mention another motivation for studying topological defects in con- 
nection with brane-world scenarios. In field theories in flat space-time which allow for 
topological defect solutions, the conservation of a topological charge guarantees the sta- 
bility of the configuration. When gravity is included, it will of course not be possible 
to conclude stability by topological arguments only, but at least in cases where gravity 
is weak one can hope that "gravitating" topological defects can provide good potential 
candidates for stable configurations. 

Also some numerical work exists in this general context. See e.g. |14(Jj for a discussion 
of a gravitating abelian vortex in 6 dimensions and |141j for a gravitating 't Hooft- 
Polyakov monopole in 7 dimensions. The latter case will be the subject of chapter HJ 
Other related works include for example the case of a global monopole in 7 dimensions 
|142j or a texture in 5 dimensions [143 . In contrast to the corresponding 6-dimensional 
case of the global string [27], the geometry of |142j is singularity-free but the graviton 
zero-mode is not normalizable. 

As announced earlier in this section, we are now going to generalize the results of 
ref. |126j to the case of transverse spaces of dimensionality n greater than 2. We are 
restricting ourselves to transverse spaces with spherical symmetry and as we shall see 
in the following, their properties are qualitatively different for n < 2 and n > 3 extra 
dimensions, a fact which complicates the generalization. In contrast to the case with 
n > 3, for n = 1 the extra space is flat while for n = 2 the extra space can be curved, 
but is still conformally flat. 

We will consider three different possibilities: we will call the first possibility a strictly 
local defect and mean the situation where the stress-energy tensor of the defect is zero 
outside the core or, for the more realistic situation of a non-singular brane of finite 
thickness, exponentially vanishing outside of the defect core. In this case, Einstein's 
equations do not allow a simple generalization of the geometries for n = 1 and n = 2 
leading to the localization of gravity. 

The second possibility is related to the so-called global defects. In this case one 
assumes a multiplet of scalar field, with a global symmetry e.g. O(N) which is spon- 
taneously broken. Outside the core of the defect this field may have a hedgehog type 
configuration, which gives a specific contribution to the stress-energy tensor. As men- 
tioned above, this case was studied in fZT\ I124j for n = 2 and in |13Uj for higher di- 
mensions, where the solutions with an exponential warp factor were found. We compute 
the corrections to Newton's law in this case and define a generalization of these metric 
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solutions. 

The third possibility is related to configurations of the monopole type. In this case, 
outside of the defect, the stress-energy tensor is dominated by the energy of a magnetic 
field or its generalization to higher dimensions, p-form abelian gauge fields. We consider 
different spherically symmetric configurations and define those which lead to gravity 
localization on the 3-brane and a regular geometry in the bulk. We also discuss the 
corrections to Newton's law for these solutions. 



3.3 Einstein equations with a 3-brane source 

In D-dimensions the Einstein equations with a bulk cosmological constant and a 
stress-energy tensor Tab are 

RAB ~ ^9AbR + J^^gAB = J^2 Tab ' ( 3 - 3 ) 

where Mp is the reduced D-dimensional Planck scale. We will assume that there exists 
a solution that respects 4d Poincare invariance. A D-dimensional metric with D = A + n 
satisfying this ansatz for n transverse, spherical coordinates with < p < oo, < 
. . . , 8 2 } < 7r and < Q x < 2n is 

ds 2 = a(p)g liu dx fl dx u + dp 2 + -f(p)dn 2 n _ 1 , (3.4) 

where as usual we follow the conventions of |112j with the metric signature of given 
by (—,+,+, +). dQ 2 l _ 1 is defined recursively as 

dn 2 ^ = de 2 n ^ + sin 2 e n ^dn 2 n _ 2 , (3.5) 

with dttl = 0. We assume in the following that the 3-brane is located at the origin p = 
of the above coordinate system ()3.4|) and that its source is described by a stress-energy 
tensor T A B with nonzero components 

n = ^/o(p), TP p = f p (p), and T% = f e (p). (3.6) 

Here we introduced three source functions / 0) fp-, an d fe which depend only on the radial 
coordinate p and by spherical symmetry all the angular source functions are identical. 
By definition we have 9 = n _i. Using the metric ansatz ()3.4|) and the stress-energy 
tensor (|3.6|) . the Einstein equations become 

— + -{n - 1 -- + « (n - 1 (n - 4 )± J + - n - 1 -L - — n - 1 n - 2 
2cr4 (T78 7^2 7 27 

M l^»-fM] + ^fl, (3.7) 
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-— + (n - 1 -J- + - (n -In- 2 JL - — n -In- 2 
2 (T cr 7 8 7 Z 27 

= AMI + (a. 

2- + -— + (n - 2 -J- + - (n - 2 (n - 5 ^ + - (n - 2 J- - — (n - 2 (n - 3 
o 2 <y a 7 8 7 Z 2 7 27 

1 - . n / \ n 2Ap/ l y S 1 



where ' denotes differentiation with respect to p and the Einstein equations arising from 
all the angular components simply reduce to one angular equation (j3.9j) . The constant 
Aphj/s represents the physical 4-dimensional cosmological constant, where 

< ) -^ (4) + ^ £ ^ = 0. (3.10) 

The system of equations (|3.7|) - (J3.9j) for fi = was first derived in [T3] and describes the 
generalization of the setup considered in |T3J HOI H26| H3Uj , to the case where there are 
n transverse dimensions, together with a nonzero cosmological constant in 4-dimensions 
and stress-energy tensor in the bulk. If we eliminate two of the equations in ()3.7j) - (j3.9j) 
then the source functions satisfy 

/; = 2-(/ - f p ) + —^(fe - f P ) , (3.11) 
a ^7 

which is simply a consequence of the conservation of the stress-energy tensor W M T% = 0. 
In general, the Ricci scalar corresponding to the metric ansatz ()3.4j) is 

a" a' 2 a 1 i 7" 1 V 2 

R = -4 2(71 - 1 -- " (n - 1 — " t(" -1 

CT (J 2 a 7 7 4 7 2 

+ (n-l)(n-2)i + i%^. (3.12) 

7 cr iWp 

The boundary conditions at the origin of the transverse space are assumed to be 

U = °> (V7~)'U = 1 and 7 | p=0 = 0, (3.13) 



cr 



which is consistent with the usual regular solution in flat space. We also set c(0) = A, 
where A is constant. We can now define various components of the brane tension per 
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unit length as 



2 (n-l)/2 , 



fiifi) 



(3.14) 



fii = - / dpa j ( 
o 

where z = 0,p,9. Using the system of equations ()3.7|) - (J3.9J) we obtain the following 
relations 



(n - 2)/z - n p -{n- l)p 6 



2A D - 



{n + 2)A phy s Mp 
M 2 P 



n+2 



cr 2 V7" 



-1 



dp (3.15) 



and 



oo 

4/i + /i p - 3/i e =(n + 2)M£ +2 A 2 5„ 2 + (n + 2)(n - 2)M£+ 2 / a 2 y/r^dp 

o 

oo 

~ 2A dJ a 2 \/T- l dp. (3.16) 
o 

In the derivation of relations (|3.15J) and ()3.16|) we used the boundary conditions ()3.13|) 
together with the assumptions that all integrals converge and that furthermore 



lim oo \[ 



p— >oo 



7 



71— 1 



lim a 2 v / 7"- 2 (v / 7) / = 0. 



p— >oo 



(3.17) 



Since we are mainly interested in solutions that localize gravity, these assumptions are 
physically reasonable. As we will see below, in most of the cases where the condition of 
localization of gravity is satisfied also the above assumptions are valid. The equations 
()3.15|) and ()3.16|) are the general conditions relating the brane tension components to the 
metric solution of the Einstein equations (|3.7j) - (j3.9|) and lead to nontrivial relationships 
between the components of the brane tension per unit length. Furthermore, by analogy 
with the solution for local strings we can identify ()3.15|) as the gravitational mass per 
unit length and (j3.16|) as the angular deficit per unit length. Thus the source for the 

3- brane, in general curves the transverse space. 

From the Einstein term in the D-dimensional Lagrangian we can obtain the effective 

4- dimensional reduced Planck mass. Using the spherically symmetric metric ansatz (|3.4jl . 
the 4-dimensional reduced Planck mass is given by 



Mp = A n - X M n D 



n+2 



dp a 7 



(n-l)/2 



(3.18) 
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where A n -\ = 2tt^ /r[~] is the surface area of an n — 1 dimensional unit sphere. We are 
interested in obtaining solutions to the Einstein equations ()3.7j) - (j3.9j) such that a finite 
four- dimensional Planck mass is obtained. This leads to various possible asymptotic 
behavior for the metric factors a and 7 in the limit p — > 00. Below, we will concentrate 
only on the case of vanishing 4-dimensional cosmological constant, A p h ys = 0. 



3.4 Strictly local defect solutions 

First, let us assume that the functions fi{p) vanish sufficiently rapid outside the core of 
the topological defect. In order to obtain a finite 4-dimensional Planck scale, we require 
a solution of the system of equations (J3.7j) - (J3.9j) for which the function a^' 1 ^ 2 goes 
to zero when p — > 00. For n = 1 and n = 2 the solutions are known to exist, see 
|4*U] and |126j correspondingly. However, when n > 3 the structure of the equations is 
qualitatively different because of the presence of the I/7 terms in (|3.7j) - ()3.9|) . Thus, 
there is no simple generalization of the solutions found for n = 1 and 2. 

To neutralize the effect of the I/7 term, one can look for asymptotic solutions for 
which 7 is a positive constant. However, one can easily check that the system of equations 
()3.7j) - (j3.9p does not allow a solution for which 7 tends to a constant when p — ► 00, and 
a is a negative exponential. 

Alternatively, we can assume that there is an asymptotic solution for which 7 — > 00 
but a tends to zero faster than ^(™ -1 )/ 2 and omit the troublesome I/7 term from the 
equations of motion. In this case the set of equations (|3.7|) - (J3.9j) can be simply reduced 
to a single equation, as in 6-dimensional case [T^L 126J: 

z» = -^±, (3.19) 
dz 



where the potential U(z) is given by 



(n + 3) A D ^ 
4(n + 2) M r D 



U(*) = ^k^»f- (3.20) 



With this parametrization the metric functions er(p) and 'j(p) can be written in terms of 
z(p) as 



a = | jZ /|(2-^/(n+2)(n-l))/(n+3) | z | (2+ y/(n+2)(n-l))/(n+3) {3.21) 
7 = | z /|6/(l-n+2 v /(n+2)(n-l)) | jZ |6/(l-n-2 v /(n+2)(n-l)) _ ^3^2) 
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Solving equation (|3.19|) with the potential (|3.2(J|) gives the general solution 



z(p) = die 



\(n+3)cp _|_ ^ e j(n+3)cp 



(3.23) 



where c 2 = —8An/((n + 2)(n + 3)M^ +2 ), di, d 2 are constants and we take A^> < 0. This 
solution is a generalization of the pure exponential solution considered earlier and in 
ref. |130j . see below. In this picture we can think of particle motion under the influence 
of the potential ()3.2()j) with position z(p) and "time" p. 

Since 1 - n + 2 a/0 + 2){n - 1) > and 1 - n - 2yJ(n + 2)(n - 1) < 0, the metric 
factor 7 can be large in two cases. In the first case z(p) is zero for some po- However, this 
point is only a coordinate singularity (the Ricci scalar is regular at this point) and the 
metric can be extended beyond po, leading then to an exponentially rising solution for 
both cr and 7. This, unfortunately, is not interesting for compactification. In the second 
case both z' and z are non-zero and increase exponentially for large p. Thus, there is no 
possibility of a finite Planck scale in this case either. 

Similarly, we were unable to find solutions in the reverse case when 7 vanishes at 
infinity. Moreover, even if such solutions were to exist, they would likely lead to a 
singular geometry (naked singularity), because the Ricci scalar contains a I/7 term, see 
eq. (J3.12)) . This was indeed shown to be the case for solutions with regular geometries 
at p = in 

3.5 Bulk scalar field 

3.5.1 Global topological defects 

The other possibility is to consider defects with different types of "hair" i.e. with non- 
zero stress-energy tensor outside the core of the defect. We start with global topological 
defects. This case has been extensively studied in |1241 ^30 3 and we list here a number 
of explicit solutions (see also |129| 11441 H45j ). For simplicity we will restrict ourselves to 
the case of vanishing 4-dimensional cosmological constant, A p h ys = 0, and assume that 
outside the core of the defect we have the following asymptotic behavior for a: 



a(p) = e 



-cp 



(3.24) 



Consider n scalar fields a with a potential 



(3.25) 
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where v has mass dimension (n + 2)/2. Then the potential minimum is at a a = v 2 . 
The defect solution has a "hedgehog" configuration outside the core 
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( f ) a (p)= v d a , (3.26) 

d a being a unit vector in the extra dimensions, d n = cos # n _i, gP" 1 = sin Q n -\ cos # n _2, 

The scalar field gives an additional contribution to the stress-energy tensor in the 
bulk with components 

„,2 

(3.27) 
(3.28) 
(3.29) 

Now, for Ad < and v 2 > (n — 2)M^ +2 the following solution leads to the localization 
of gravity and a regular geometry in the bulk |124[ I13()j 

(3.30) 

The transverse geometry of this solution is that of a cylinder, R + x S n ^i with R + being 
the half-line and 5 n _i being an n — 1 sphere. 

When v 2 = in — 2)M]^ +2 the I/7 terms are eliminated from the system of equations 
(I3.7|) - (|3.9|) and the exponential solution to the coupled set of equations (|3.7|) - (|3.9|) can 
then be found with 




8(-Apj 
{n + 2){n + 3)M r D 



1 (p) = R 2 a(p), c= x l - ■ ^- ^— ^ , (3.32) 



where i?o is an arbitrary length scale. As expected, the negative exponential solution 
()3.24|) requires that A D < 0. Notice that the exponential solution ()3.32|) only requires the 
"hedgehog" scalar field configuration in the bulk for transverse spaces with dimension 
n > 3. No such configuration is needed for the 5-dimensional an d 6-dimensional 
cases |126j . which only require gravity in the bulk. 
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The Ricci scalar corresponding to the negative exponential solution with vanishing 
four-dimensional cosmological constant is 

c 2 e cp 
R = -in + 3)(n + 4)- + (n - l)(n - 2W , (3.33) 

4 it 

which diverges when p = oo. Thus we see that for n > 3 the space is no longer a 
constant curvature space and in fact has a singularity at p = oo. This is also confirmed 
by looking at the other curvature invariants, RabR AB and RabcdR ABCD ■ Only for the 
5d and 6d cases do we obtain a constant curvature anti-de-Sitter space. The appearance 
of a singularity is similar to the case of the global-string defect [27]. 

3.5.2 Corrections to Newton's Law 

For the solution (j3.30|) . (j3.31j) the corrections to Newton's law are parametrically the same 
as for 5d case, since 7 is a constant. On the other hand, the singular solution ()3.32|) 
will ultimately require that the singularity is smoothed by string theory corrections 
(perhaps similar to the nonsingular deformations considered in |146| ). Assuming that 
this is the case, then the corrections to Newton's law on the 3-brane can be calculated 
by generalizing the calculation presented in [4"U1 1126j (see also |147j ). 

In order to see that gravity is localized on the 3-brane, let us now consider the 
equations of motion for the linearized metric fluctuations. We will only concentrate on 
the spin-2 modes and neglect the scalar modes, which need to be taken into account 
together with the bending of the brane [IE]. The vector modes are massive as follows 
from a simple modification of the results in ref. |148j . For a fluctuation of the form 
h^ix.z) = ^(z)h IJlU (x) where z = (p,0) and d 2 h fiu (x) = mfyh^x), we can separate 
the variables by defining $(z) = '^2i. m <Pm(p)Yi i (4 ) i^i)- The radial modes satisfy the 
equation |148j 

™ d p [ay- 1 ^ 2 d p <p m ] = m 2 m , (3.34) 



(J7 (n-l)/2 

where m 2 = m 2 , + A/Rq contains the contributions from the angular momentum modes 
Zj and A denotes the Laplacian on 5' n_1 . The differential operator (|3.34j) is self-adjoint 
provided that we impose the boundary conditions 

C(0) = C(oo) = , (3.35) 

where the modes m satisfy the orthonormality condition 

00 

An-i J dpa^-W^ = 5 mn . (3.36) 
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Using the specific solution (|3.24j) . (j3.32|) . the differential operator (|3.34j) becomes 



C - C ^m + = • (3.37) 

This equation is the same as that obtained for the 5d domain wall solution [12], when 
n = 1 and the local string-like solution [126] when n = 2. We see that each extra 
transverse coordinate augments this coefficient by 1/2. When m = we clearly see that 
0o (p) = constant is a solution. Thus we have a zero-mode tensor fluctuation which is 
localized near the origin p = and is normalizable. 

The contribution from the nonzero modes will modify Newton's law on the 3-brane. In 
order to calculate this contribution we need to obtain the wave-function for the nonzero 
modes at the origin. The nonzero mass eigenvalues can be obtained by imposing the 
boundary conditions (j3.35|) on the solutions of the differential equation (j3.37|) . The 
solutions of (|3.37|) are 



= e !("+3)P 



.2,tn c 2,ui c 

CiJi M ( 7 e^)+C 2 y |(n+3) ( 7 e^) 



(3.38) 



where C 1; C 2 are constants and Ji( n+3 ), Yi( n+3 ) are the usual Bessel functions. Imposing 
the boundary conditions (j3.35|) at a finite radial distance cutoff p = p max (instead of 
p = oo) will lead to a discrete mass spectrum, where for k — 1, 2, 3, . . . we obtain 

m k ~ c(k + -) —e 2^ max . (3.39) 

With this discrete mass spectrum we find that in the limit of vanishing mass m^, 
, 9 . . 1 ixc (n + l) 2 /mt\" 

< (0) - A^Yiwrm (t) • (3 - 40) 

where T[x] is the gamma- function. On the 3-brane the gravitational potential between 
two point masses mi and m 2 , will receive a contribution from the discrete, nonzero modes 
given by 

where is Newton's constant. In the limit that p max — > oo, the spectrum becomes 
continuous and the discrete sum is converted into an integral. The contribution to the 
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gravitational potential then becomes 



V(r) ~ G 



mim 2 



n- 



= G 



N- 



1 + 



1 + 



71 + 1 



2c n+i r 2 [(n + 3)/2] 



dm m n e 



T[n + 2] 



2r 2 [(rz + 3)/2] (cr)™+ 1 



(3.42) 
(3.43) 



Thus we see that for n transverse dimensions the correction to Newton's law from the 
bulk continuum states grows like l/(cr) n+1 . This correction becomes more suppressed 
as the number of transverse dimensions grows, because now the gravitational field of the 
bulk continuum modes spreads out in more dimensions and so their effect on the 3-brane 
is weaker. 



3.6 Bulk £>-form field 

The global topological defects considered in the previous section inevitably contain mass- 
less scalar fields - Nambu-Goldstone bosons associated with the spontaneous breakdown 
of the global symmetry. Thus, the stability of these configurations is far from being 
obvious. 

We will now consider the possibility of introducing other types of bulk fields (p-form 
fields), which directly lead to a regular geometry. The stability of the corresponding con- 
figurations may be insured simply by the magnetic flux conservation. The D-dimensional 
action is 

S = J d D xV^ (^M n D +2 R — Ad — \f m _.. Ap+1 F a - a ^ . (3.44) 
The energy-momentum tensor associated with the p-form field configuration is given by 

T A B = - [\s A B F Al ... Ap+1 F A ^ - V -±±F A M ... A F B M - A ^ . (3.45) 

A solution to the equation of motion for the p-form field when p = n — 2 is 

iV..0„-! = Q(sin^_ 1 )("- 2 ) . . .sin0 2 , (3.46) 

where Q is the charge of the field configuration and all other components of F are equal 
to zero. In fact, this "hedgehog" field configuration is the generalization of the magnetic 
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field of a monopole. The stress-energy tensor associated with this p-form field in the 
bulk is 



t; = JJL^^-H, (3.47) 
T{ = -It = J ^i^£ ■ (3AS) 



Let us again assume a solution with the asymptotic form 

a{p) = e~ cp and 7 = const . (3.49) 

With this ansatz we see from the Ricci scalar that the transverse space will have a 
constant curvature and the effective four-dimensional Planck constant will be finite. If 
we substitute this ansatz and also include the contribution of the p-form bulk field to 
the stress-energy tensor, the Einstein equations (j3.7j) - (|3.9|) . with A phys = are reduced 
to the following two equations for the metric factors outside the 3-brane source 

(n-l)\ - — ( n -2)(n + 2) + -^ = 0, (3.50) 



2 = _±_^D + J_( 2 y (351) 

2M™+ 2 4 7 v ; v ' 



We are interested in the solutions of these two equations which are exponential, c 2 > 
and do not change the metric signature, 7 > 0. Remarkably, solutions to these equations 
exist for which these conditions can be simultaneously satisfied. Let us first consider the 
case n = 2, where the solutions simply reduce to 

O 2 

— = -A 6 , (3.52) 

7 

= -^f. (3.53) 

Thus, for A 6 < we see that there is a solution that satisfies ()3.49|) . 

Next we consider the case n — 3. Only the "+" solution to the quadratic equation 
(J3.50|) gives rise to a solution with both c 2 > and 7 > 0. This solution can be written 
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Figure 3.1: The n = 3 solution for Q 2 /( / yMj ) (dashed line) and {Q 2 c 2 )/M^ (solid line), 
as a function of Q 2 A 7 /M^°. Only the branches with c 2 > and 7 > lead to solutions 
that localize gravity. 

in the form 



When Q 2 Aj/Mj° < 25/32 we obtain real solutions which are plotted in Figure l3~Tl In 
fact requiring c 2 > we see that there are solutions not only for A7 < 0, but also for 
A7 > 0, provided that Q 2 At/Mj° < 1/2. Thus, the bulk cosmological constant does not 
need to be negative in order to localize gravity. 

Similarly, solutions for which c 2 > and 7 > exist for values of n > 3. Again, 
we find that solutions with both positive, zero and negative bulk cosmological constants 
exist. 

The nice property of these solutions ()3.49|) is that since 7 is a constant the Ricci 
scalar does not diverge at any point in the transverse space. In particular for the n = 2 




8Q 2 



M v 10 



(3.54) 



(3.55) 
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solution the Ricci scalar is R = 5A 6 / (2Mg), while for n = 3 it is 

(3.56) 



i? = + 7 



2Q 2 




16 ^ M 7 10 8 V M 7 10 



The space is again a constant curvature space, but it is not necessarily anti-de-Sitter. 
This is also confirmed by checking the higher-order curvature invariants RabR AB and 

p tdABCD 

Finally, we see from ()3.34|) that the equation of motion for the spin-2 radial modes 
using the solution (J3.49|) is qualitatively similar to the 5d case. The constant 7 factor 
effectively plays no role in the localization of gravity. Thus, the corrections to Newton's 
law will be suppressed by 1/r 2 for all solutions n > 3. This is easy to understand since 
the geometry of the extra dimensions is simply R + x S^-i, as in the global defect case, 
and there is just one non-compact dimension for all n > 3. In the special case when 
Ap = 0, the metric solution we have found corresponds to the near-horizon metric of a 
class of extremal non-dilatonic black branes 149J . 

Another possible solution for the p-form in the bulk includes 



where 



F^.. 4ln = e w ... Mn «(p) , (3.57) 



k( P ) = Q 1 _ (sin ^_!) (n " 2) ... sin ^ 2 . (3.58) 

<7 A 



In this case the contribution to the stress-energy tensor is 



n) O 2 
4 a 



t; = If = +j^- (3.60) 



This contribution does not appear to make the solution of the Einstein equations any 
easier. 

The above two p-form solutions have only included components in the transverse 
space. If we also require the p-form field to transform non-trivially under the 3-brane 
coordinates, then we can have 



012301. ..e„ 



^ 1 =Q(sm6 n _ i y n - 2 K..sm6 2 , (3.6i; 
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where all other components of F are zero. The components of the stress-energy tensor 
for this field configuration are 



TP = -T e e = + {n + 3)l f - . (3.63) 



Again, there is no simple solution of the Einstein equations with the inclusion of this 
contribution. 

Finally, one can also generalize the solution (J3.57JI 

F fj.ua/3aL.. a„ = ^fiua/3a 1 ...a n K (p) , (3.64) 

where 

K {p) = Qa 2 ^ n - 1)/2 (sm6 n ^) {n - 2) . . .sm9 2 , (3.65) 
and the stress-energy tensor 

rpM (n + 4). 2 rM /q C( ;\ 

TV = 4 <3 <V ) ( 3 - 66 ) 

is a constant. Thus, we see that the addition of a n + 4 form field is equivalent to adding 
a bulk cosmological constant. 



3.7 Conclusions 

In summary we can say that the generalization of the exponential solutions found in |40j 
and |126j in the case of one or two extra dimensions which lead to the localization of 
gravity on a 3-brane is nontrivial. In the case of a strictly local defect, where the stress- 
energy tensor vanishes exponentially outside the core of the defect, we were unable to find 
any solution to the Einstein equations. However, when a scalar field with "hedgehog" 
configuration is added to the bulk, there exist solutions which localize gravity. In the 
case (13.32)1 of exponentially vanishing 7, the transverse space no longer has constant 
curvature, and furthermore develops a singularity at p = 00. The corrections to Newton's 
law on the 3-brane are suppressed by l/(cr) n+1 , where we see that each additional extra 
dimension adds one more power of 1/r to the correction factor. If however 7 is a constant 
fl3.30|) - ()3.31|) then regular gravity localizing solutions exist. 

By adding a n — 1 form field in the bulk, which can be considered as the generalization 
of the magnetic field of a monopole, gravity localizing solutions can be found for positive, 
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zero and negative bulk cosmological constant. In this case, the transverse space has 
constant curvature but is not an anti-de-Sitter space. Since 7 approaches a constant at 
infinity in transverse space, the corrections to Newton's law have the same form as in 
the original Randall-Sundrum II model [IU]. Eventually we saw that the addition of an 
n + 4- form field in the bulk is equivalent to adding a bulk cosmological constant. 

In the next chapter we will consider a particular field theoretical model in 7 dimen- 
sions in which the asymptotic solution (J3.49j) found in this section is realized. More 
concretely, we will discuss numerical solutions to the coupled Einstein and matter-field 
equations for a 't Hooft-Polyakov monopole residing in a 3-dimensional spherically sym- 
metric transverse space. 
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Chapter 4 



A 
't 



numerical example - the 
Hooft-Polyakov monopole 



4.1 



Introduction 



This chapter is dedicated to a detailed numerical analysis of the coupled Einstein and 
matter field equations for a particular topological defect, a 't Hooft-Polyakov monopole, 
residing in a 3-dimensional transverse space of a 7-dimensional space-time. The motiva- 
tions for this are twofold: first, as we outlined briefly in the last chapter, a more involved 
topological background in the extra dimensional transverse space provides richer struc- 
ture to generate chiral fermions on the brane as zero modes of higher- dimensional non- 
chiral spinor-fields. From this point of view, such considerations may help to envisage the 
construction of realistic models. Second, a common feature of the solutions to Einstein 
equations presented in the last chapter was their asymptotic nature, valid only outside of 
the core of the defect. It is not a priori clear, whether the fully coupled system of Einstein 
plus matter field equations allows for singularity free solutions that interpolate smoothly 
between the topological defect at the origin of transverse space and the asymptotic at 
infinity in transverse space. It was shown recently |14Uj (see also [150, 151"! 1152j ) that 
gravity can be localized on an abelian Higgs string in 6 dimensions in the absence of 
any singularity under the assumption of a certain amount of fine tuning between brane 
tension components. The aim of the numerical work presented here is to apply the 
same idea to the case of a monopole in 7 dimensional space-time. This chapter is orga- 
nized as follows: in the following section 14.21 we review the physics of 't Hooft-Polyakov 
monopoles in flat 4-dimensional Minkowski space-time. In section 14.3.11 we present the 
SU(2) invariant Georgi-Glashow toy model (having 't Hooft-Polyakov monopoles as flat 
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space-time solutions |2Zj,j2Hl) coupled to gravity in 7 dimensions. The Einstein equa- 
tions and the field equations are obtained in the case of a generalized 't Hooft-Polyakov 
ansatz for the fields. While boundary conditions and the desired asymptotic of the so- 
lutions are discussed in sections 14.3.21 and 14.3.31 we give the relations between the brane 
tension components necessary for gravity localization in section 14.3.51 A proof of these 
relations is given in appendix [BJ Section I4~4l is dedicated to numerics: after presenting 
the numerical methods used for solving the differential equations in section 14.4.11 we 
discuss the fine-tuning surface (the relation between the independent parameters of the 
model) in section 14.4.21 We give explicit sample solutions in section 14.4.31 and discuss 
their general dependence on the parameters of the model. Sections 14.51 and 14.61 finally 
treat the Prasad- Sommerfield limit and the fine-monopole limit, respectively. While in 
the former case no gravity localizing solutions exist, in the latter case we demonstrate 
the possibility of physical choices of the model parameters which are in agreement with 
observations. 

4.2 The 't Hooft-Polyakov monopole in flat 4d space- 
time 

As we pointed out in the last chapter, from a field theoretical point of view monopoles 
can be characterized as non-singular, point-like (O-dimensional), finite-energy solutions 
to classical field equations. Their interest for particle physics stems mainly from their 
stability which again can be understood on mere topological grounds. Excellent reviews 
and textbooks can be found on this subject, see for example [87| ITUUl I10H I102| [103 . 
In the following we summarize the basic properties of 't Hooft-Polyakov monopoles in 
4-dimensional Minkowski space-time. Monopoles are produced during the process of 
spontaneous symmetry breaking whenever the vacuum manifold M. = G/H of the Higgs 
field after the symmetry breaking contains non-contractible two surfaces ^(G/H) ^ I. 
Making use of the second fundamental theorem (j3.2j) of section 13.1.21 we see that this 
can be rephrased as 7Ti (Hq) ^ I with Hq being the component of H connected to the 
identity. This is the reason why monopoles are a general prediction of grand unified 
theories whenever a single group is to encompass all interactions. At some stage of the 
symmetry breaking U (l) em , the gauge group of electromagnetism has to appear and with 
it the associated monopoles since 7i"i (U(l)) = Z. The existence of monopole solutions 
was established by the pioneering works of 't Hooft |JZ] and Polyakov [§H] taking the 
simple example of a model proposed by Georgi and Glashow 153] in the context of 
electroweak unification. The underlying field theory is a gauge theory based on the 
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group 577(2). It contains a scalar triplet $ a transforming in the adjoint representation 
of 577(2) and the corresponding gauge field W®. Group indices are denoted by a tilde : 
a,b,c = 1 ... 3. The action for this theory is given by 



S 



d x 



A 



-2\2 



■-Gl u G^ u - -£>„$ a £>"$ a - - ($ a $ a - ri 2 ) 



where 
and 



/f/ S -I- p p 5 ' 
•/xi/ ~7»" i/ ^f"^ Tec , IjU , 



(4.1) 



(4.2) 



(4.3) 



The variation of the action (j4.1j) gives the following equations of motion for the gauge 
and scalar fields: 



-e e 



(4.4) 
(4.5) 



In this model the SU(2) symmetry is spontaneously broken down to U(l) which can 
be identified with the gauge group of electromagnetism [§7| . The monopole corresponds 
to the simplest, topologically non-trivial field configuration (with topological charge 1). 
In terms of the parameters rj and A, the Higgs mass is given by m# = rj \^2X. Two of 
the three vector fields acquire a mass given by raw = erj. 

We are interested in static, spherically symmetric solutions to the set of coupled 
equations ()4.4j) - (j4.5j) with finite energy. The energy functional for static configurations 
can be written as: 



E 



d x 



2\ 2 



(4.6) 



From (|4.6|) we can see that the condition of finite energy implies that the fields fall off 
at infinity sufficiently fast: 



Vi<&~ a -> 0, and -> rj 2 for r 



x 



oo . 



(4.7) 



The second condition for the Higgs triplet $ a above shows that in the limit r — > oo 
the field $ can be thought of as a function from physical infinity of topology S 2 to 
the vacuum manifold G/H also of topology S 2 . This argument is at the base of the 
homotopy classification of monopoles solutions in this theory since the above mappings 
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fall into homotopy classes labeled by an integer n according to ^{S 2 ) = Z. The so-called 
't Hooft-Polyakov ansatz respects the above boundary conditions for n = 1: 

= **-—, (4.8) 
r er 

W?(x i )=e aij - 1 ~ K ( T) , Wg = 0, (4.9) 

r er 

provided we have 

H(r) 

lim K{r) = , lim — ^ = 77. (4.10) 

1 — >oo r— >oo er 

Due to the simplicity of the above ansatz, the field equations (|4.4j) - (j4.5|) reduce to the 
following system of ordinary differential equations: 



2 



d 2 H(r) 

r 2 ,\> = 2H(r)K(r) 2 + XH(r) 
dr l 



H(r) 2 2 2 
— 5 r 77 



(4.12) 



No analytic solutions to the system 1)4.11)1 and 1)4.12)1 are known at present except for 
the case A — > 0. This limit is known as the Prasad- Sommerfield limit |154) 1155) and the 
solutions to (j4.11j) and (j4.12j) take the particular simple form: 

K(r)= H(r) = . - 1 . (4.13) 

smh{rerj) tanh(rer/J 

Even though 't Hooft-Polyakov monopoles give rise to plenty of interesting physics, we 
will stop our discussion at this point and come back to our main subject, the gravitating 
monopole in a 7-dimensional space-time. 

4.3 The warped 't Hooft-Polyakov monopole in 7d 
curved space-time 

4.3.1 Einstein equations - field equations 

The action for the setup considered in this section is a straightforward generalization 
of the action of a gravitating 't Hooft-Polyakov monopole in 4 dimensions, objects that 
have been extensively studied in the past )156) 1157) 1158) 1159) 1160) , to the case of a 7 
dimensional space-time: 

c = c _|_ q (a -\A\ 

LJ ^gravity >^ brane • \ / 
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S gravity is the 7-dimensional generalization of the Einstein-Hilbert action: 

S gramty = [ dJxyf^ \R - . (4.15) 



Here, g is the determinant of the metric guN with signature mostly plus — K . . +. Upper 
case latin indices M, N will run over ... 6, lower case latin indices m, n over 4 ... 6 and 
greek indices //, v over ... 3. M 7 denotes the fundamental gravity scale and A 7 the bulk 
cosmological constant. S brme is formally identical to the Georgi-Glashow action ()4.1|) in 
Minkowski space-time apart from the dimensionality: 



Sbrane j d £"r. 



(4.16) 



with 



£m = -\G\ iN G aMN - ^V M &V M & - ~ - rff . (4.17) 

Here, rj is the vacuum expectation value of the scalar field. Following the same conven- 
tions as in the last section, we denote group indices by a tilde: a,b,c = 1 ... 3. T>m $ a 
denotes the gauge and generally covariant derivative of the scalar field 

V M & a = d M ® & + ee al ~ c W l M &. (4.18) 

Again we have 

G\ [N = d M W% - d N Wl + e e sSe W h M W%. (4.19) 

The general coupled system of Einsteins equations and the equations of motion for the 
scalar field and the gauge field following from the above action are 

Run — -gMNR + j^gMN = j^T M n , (4.20) 



=V M (^/=g~V M <S> i ) = A $ 5 (V# - t? 2 ) , (4.21) 
— V M {V=gG &MN ) = -ee^ (v N <fP) $ c ~, (4.22) 



-g 

where the stress-energy tensor T MN is given by 

T MN = - J^f g MN = G mlG% L + V M ^V N ^ + g MN C m . (4.23) 
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We are interested in static monopole-like solutions to the set of equations (|4.2(J|) - (j4.22|) 
respecting both 4D-Poincare invariance on the brane and rotational invariance in the 
transverse space. The fields $ a and (and as a result G a MN ) should not depend on 
the coordinates on the brane x^. The brane is supposed to be located at the center of 
the magnetic monopole. A general non-factorizable ansatz for the metric satisfying the 
above conditions is 

ds 2 = M 2 {p)g^dx^dx u + dp 2 + L 2 {p) (d6 2 + sin 2 9 dp 2 ) , (4.24) 

where g$ is interpreted as the metric of 4D general relativity which satisfies the 4D- 
Einstein equations with an arbitrary cosmo logical constant A p h ya , also appearing in the 
7D-Einstein equations ^3], I n the following we will only consider the case A p h ys = 
and we take g$ to be the Minkowski metric rj^ with signature ( — h ++). Since the 
choice A p hys = is not motivated by any physical argument, a possible solution to the 
cosmological constant problem depends on a mechanism explaining why the value is 
preferred. The 't Hooft-Polyakov ansatz (J4.8j) (see jSJj) for the fields is readily generalized 
to our setup. After transforming to spherical coordinates this ansatz becomes 

W p = 0, W e = - 1 -^^e-, ^i-^sinfle*, W, = , 

e e 

l=^e ? , (4.25) 
ep 

with flashes indicating vectors in internal SU(2) space and with 

e p = (sin 9 cos ip, sin 9 sin ip, cos 9) , (4.26) 
e e = (cos 9 cos if, cos 9 sin ip, — sin 9) , (4.27) 
(— sinyj, cos<^, 0) . (4.28) 



"p 



Using the ansatz ()4.25|) for the fields together with the metric (|4.24|) in the coupled 
system of differential equations ()4.20j) - (j4.22j) gives 



M" C" M'C M C 1 1 



M' £' 


+ 3 


M' 2 




C' 2 




M 2 


+ 




ML 




M' £ 


+ 6 


M' 2 








M 2 


+ 




ML 


C 2 



± = (3(e p -A 7 ) , (4.30) 
M" C" M'C M' 2 n . , . , 

'ir + t +4 w +6 m? = ^- A ''- (431) 
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J " + K 2 ^ + £) J '~ 2 ^ = o.J(.P-l), (4.32) 
K » + K ^-f) +i M K . = fiK . (4.33) 

where primes denote derivatives with respect to the transverse radial coordinate r rescaled 
by the mass of the gauge boson mw' 



r = m w p = r/e p. 



(4.34) 



All quantities appearing in the above equations are dimensionless, including a,/3 and 7. 
We have 



A 1 ( m H 



e 2 2 \m w 



P M 7 5 ' 7 



e 2 rf ' 



fi 



e 2 r) 4 ' 



with 



(4.35) 



T^ = ^/ , TP = f p , T e 6 = f e , T% = f v , C = Lm w , J 



H 

rje p 



(4.36) 



The dimensionless diagonal elements of the stress-energy tensor are given by 



e 



K 



r> 



'1-K 2 



K >2 (1 _ ^2)2 



1 T ,2 , J 2 K 2 
2 J + C 2 



+ ?(^-l) S 



J 2 ^ 2 
^C 2 



£ 2 

(1-K 2 
2XF 



2£ 4 



-J' 



/2 



/2 



2 

a: 



a 



(J 2 -1) : 



(4.37) 

(4.38) 
(4.39) 



The rotational symmetry in transverse space implies that the (6 9) and the (<p<p) compo- 
nents of the Einstein equations are identical (and that e v = eg). Equations (j4.29J) - (j4.31|) 
are not functionally independent 132 . They are related by the Bianchi identities (or 
equivalently by conservation of stress-energy VmT^ 1 = 0). When solving the system 
(|4.29|) - (j4.33|) numerically, eq. (|4.3U|) can be considered as a constraint which is automat- 
ically satisfied. 

The Ricci scalar and the curvature invariants R, Rab R ab , Rabcd R abcd , Cabcd C abcd 
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with Cabcd the Weyl tensor are given by 
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(4.40) 
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Following |132j we can define various brane tension components by 



IH = -J drM{rYC{rfei{r) 
o 



(4.41) 



(4.42) 



(4.43) 



(4.44) 



4.3.2 Boundary conditions 

The boundary conditions should lead to a regular solution at the origin. Thus we have 
to impose 

M| r=0 = l, M'| r=0 = 0, £| r=0 = 0, £'| r=0 = l, (4.45) 
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where the value +1 for M\ r=0 is a convenient choice that can be obtained by rescaling 
brane coordinates. 

Both the vacuum manifold of the scalar triplet $ a and infinity in transverse space 
(r = oo) are homeomorphic to S 2 , so that all possible field configurations with specified 
boundary conditions at infinity can be classified by a topological (winding) number, see 
e.g. jHZj or jHHj. For the monopole solutions we are interested in, this winding number is 
1. The scalar field given by the ansatz ()4.25|) is radial at infinity in the bulk. Unwinding 
such a field configuration is impossible, as is continuously contracting a smooth "identity" 
mapping S 2 —>■ S 2 to a mapping where the whole S 2 is mapped into a single point. The 
radial component of the scalar field J(r) must have at least one zero [HZj which is chosen 
to be at the origin (r = 0) of transverse space. Thus 

J(0) = , lim J(r) = 1 . (4.46) 

r— >oo 

For the gauge field the condition of vanishing covariant derivative of $ a at infinity and 
the regularity of the ansatz ()4.25|) gives 

K(0) = 1 , lim K(r) = 0. (4.47) 

r— >oo 

Up to this point our considerations were similar to those of |14Uj . However, the question 
of localizing gravity differs from the string case in several ways. Most importantly, in 
order to obtain solutions that localize gravity, i.e. solutions for which 

oo 

Mp = f M(r) 2 C{rfdr (4.48) 

"% J 
o 

is finite, it is impossible that M(r) and C(r) should both be exponentially decreasing 
in the bulk. No solutions corresponding to a strictly local defect have been found for 
n = 3 extra dimensions, see |132j and the discussion in the last chapter. This follows 
from the presence of the l/£(r) 2 -term in the Einstein equations. The solutions we 
are looking for are concrete (numerical) realizations of the metric asymptotic ()3.49j) . 
(j3.54j) . ()3.55|) found in the last section for the more general case of a p-Form field: 
M(r) = M e~2 r ; C(r) = const. Another important difference from the 5D case [40J and 
the 6D case |14Uj is that even though the space-time has constant curvature (for r — > oo), 
it is not necessarily anti-de-Sitter. Finally, both signs of the bulk cosmological constant 
are possible. 
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4.3.3 Asymptotic of the solutions - behavior at the origin 

Once boundary conditions at the center of the defect are imposed for the fields and the 
metric, the system of equations (|4.29J) - (|4.33|) can be solved in the vicinity of the origin 
by developing the fields and the metric into a power series in the (reduced) transverse 
radial variable r. For the given system this can be done up to any desired order. We 
give the power series up to third order in r: 



It can easily be shown that the power series of M(r) and K(r) only involve even pow- 
ers of r whereas those of C(r) and J(r) involve only odd ones. The expressions for 
C{r) and J{r) are therefore valid up to 5 th order. One observes that the solutions sat- 
isfying the boundary conditions at the origin can be parametrized by five parameters 
[at, (3, 7, J'(0), K"(0)]. For arbitrary combinations of these parameters the correspond- 
ing metric solution will not satisfy the boundary conditions at infinity. Therefore, the 
numerical task is to find those parameter combinations for which ()4.48|) is finite. For 
completeness, we give the power series expansion at the origin to zero-th order for the 
stress-energy tensor components 




e p | r =o — £e\r=o 



£()|r=0 




(4.53) 



(4.54) 



and the curvature invariants: 





7a + 28 7 + 30 J'(0) 2 + 18fr"(0) 2 



(4.55) 




7a 2 + 56a7 + 112 7 2 + 60a J'(0) 2 + 240 7 J'(0) 2 + 300 J'(0) 4 (4.56) 
+12 (3a- + 12 7 + 70 J'(0) 2 ) iT / (0) 2 + 732ff , (0) 4 , 
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17 a 2 + 136 a 7 + 272 7 2 - 60 a J'(0) 2 - 240 7 J'(0) 2 + 900 J'(0) 4 

(4.57) 
(4.58) 



-36 (9 a + 36 7- 110 J'(0) 2 j K"(0) 2 + 4932 K"(0) 4 
a + 4 7 _ 6 (j'( f + 3K"(0f 12 



4.3.4 Asymptotic of the solutions - behavior at infinity 

The asymptotic of the metric functions M(r) and C(r) far from the monopole is: 

M(r) = M e-2 r and C{r) = C = const , (4.59) 

where only positive values of c lead to gravity localization. This induces the following 
asymptotic for the stress-energy components and the various curvature invariants: 
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(4.61) 
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(4.62) 
(4.63) 
(4.64) 
(4.65) 



c and Cq are determined by Einsteins equations for large r and the corresponding asymp- 
totic was derived in the last chapter in ()3.54j) - (j3.55j) . In the present case these relations 
become: 



32M 1 " 



1 " 1^ 



-1 = JL 1 1 

C\ 8/? 



1 - - 7 /3 2 
25 lh 



(4.66) 
(4.67) 
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The charge of the field configuration can either be determined by comparing the energy- 
momentum tensor with the general expression given in (J3.46)) of section I3~T)1 or by a direct 
calculation of the magnetic field strength tensor following 't Hooft |HZ] : 



MN 



l$l 



e|$| 3 



(4.68) 



The only nonzero component of Qmn is Ge<p = — nfr- Either way gives Q = ~. In 
order to obtain some information about the asymptotic behavior of the fields J(r) and 
K(r) at infinity we insert relations 1)4.59)1 into 1)4.32)1 and (|4.33j) . Furthermore we use 
K(r) = SK(r) and J(r) = 1 — 5J(r) with 5K <C 1 and 5 J <C 1 for 1 <C r to linearize 
these equations: 



5 J" -2c5J'-2a5J 







5K"-2c5K'+ ( -j^- I j 6K 



(4.69) 
(4.70) 



By using the ansatz 5K = Ae kr and 5 J = Be J> we find 



-C ± W C 2 - I ^2 - 1 



ji,2 = -c ± Vc 2 + 2a 



(4.71) 



To satisfy the boundary conditions we obviously have to impose k > and j > 0. We 
distinguish two cases: 



1. 



c > Gravity localizing solutions. In this case there is a unique k for -h > 1 with 



the positive sign in (j4.71j) . 



(b) 
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> 1 




1 

'-o 


< 1 




1 


= 1 





2. | c < | Solutions that do not localize gravity. 

Both solutions of 1)4.71)) are positive. 
In this case there are no real solutions. 



c 2 
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- 1 




c 2 


< 
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- 1 



There is a unique solution with the positive sign in (j4.71)) 
The important fc-value here is k = —2c. 
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It is easily shown that for large enough r the linear approximation (|4.7U|) to the equation 
of motion of the gauge field is always valid. Linearizing the equation for the scalar field 
however breaks down for 2k < j due to the presence of the term oc K 2 J in (J4.32)) . In 
that case J(r) approaches 1 as e~ 2kr . For solutions with c > (the case of predominant 
interest) a detailed discussion of the validity of 2k > j gives: 



1. a = 



Prasad- Sommerfield limit 
scalar field is governed by e~ J> . 



. One has 2k > j. The asymptotic of the 



2. < a < 2 The validity of 2k > j depends on different inequalities between a, 



-ps and c 2 . 



a + -k < 1 



2k > j 



(b) 



a 



4>1 



and 



In this case 2k > j is equivalent to c 2 < 



2 £ ; 



1-a- 



where equality in one 



of the equations implies equality in the other. 

-> 2k < j and the scalar field asymptotic is governed by e~ 2kr . 



f + ^>1 



3. a > 2 



2k < j. Also in this case the gauge field K{r) determines the asymp- 



totic of the scalar field J{r). 



This information is essential for obtaining numerical solutions to the equations ()4.29|) - 
1)4.33)1 since the integration has to be stopped at some finite distance from the origin 
where boundary conditions have to be imposed. 



4.3.5 Fine-tuning relations 

It is possible to derive analytic relations between the different components of the brane 
tension valid for gravity localizing solutions. Integrating linear combinations of Einsteins 
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equations (|4.29|) - (j4.31|) between and oo gives after multiplication with M 4 (r)£ 2 (r): 



If f M 4 

to -He = pjM 4 dr = j(l-K 2 )—dr, (4.72) 



oo 

to -»p- 2to = 2 7 J M A C?dr , (4.73) 

o 

oo 

fx + fi p + 2fi e = a J (l- J 2 ) M 4 C 2 dr. (4.74) 



To obtain the above relations integration by parts was used where the boundary terms 
dropped due to the boundary conditions given above. A detailed derivation of these 
relations is given in appendix iBl 



4.4 Numerics 



4.4.1 General discussion 

As already pointed out, the numerical problem is to find those solutions to the system 
of differential equations ()4.29|) - (j4.33|) and boundary conditions for which the integral 
defining the 4-dimensional Planck-scale (J4.48)) is finite. This is a two point boundary 
value problem on the interval r = [0, oo) depending on three independent parameters 
(a,(3,j). Independently of the numerical method employed, the system of equations 
(I4.29|) - (j4.33|) was rewritten in a different way in order for the integration to be as stable as 
possible. By introducing the derivatives of the unknown functions J(r), K(r), M(r) and 
C(r) as new dependent variables one obtains a system of ordinary first order equations. 
In the case of M(r) it turned out to be convenient to define y-j = M'(r)/M(r) as a new 
unknown function (rather than M'(r)) since the boundary condition for y 7 at infinity 
then simply reads 

limj/ 7 = ~. (4.75) 

r— >oo / 
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With the following definitions, we give the form of the equations which is at the base of 
several numerical methods employed: 
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— ^ - 42/62/7 + 2/12/2 , 

-4y 7 2 - 2</ 7 g - f (2 7 + e - e p - 2e„) , 
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with 



£o _ £p _ 2£ , = | (v? _ ir _(i^I)!_ 2 |, (4.77) 

Z 2/4 2/4 

-4e -e p + 3e e = £ (y? - l) 2 + 4 fl-^ + 5 V M + 3 4 . (4.78) 

2 2/4 yl yl 

This is an autonomous, ordinary system of coupled differential equations depending on 
the parameters (a, (3, 7). The boundary conditions are 

2/i(0) = 0, y 2 (0) = 1, y 3 ( ) = 1, y 4 (0) = 0, 

limyi(r) = 1, limy 2 (r) = 0, y 7 (0) = 0, y 8 (0) = 1. 1 ' 

r— >oo r— >oo x ' v ' 

In order to find solutions with the desired metric asymptotic at infinity, it is useful 
to define either one (or more) of the parameters (a, j3, 7) or the constants J'(0) and 
K"(0) as additional dependent variables, e. g. 2/9( r ) — a with y' 9 (r) = 0, see |161j . 
Before discussing the different methods that were used we give some common numerical 
problems encountered. 

• Technically it is impossible to integrate to infinity. The possibility of compactifying 
the independent variable r was not believed to simplify the numerics. Therefore, 
the integration has to be stopped at some upper value of r = r max . For most of the 
solutions this was about 20 ~ 30. 

• Forward integration with arbitrary values of (a, (3, 7, J'(0),K"(0)) turned out 
to be very unstable. This means that even before some integration routine (e.g. 
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the Runge-Kutta method reached r max , the values of some went out of 

range which was due to the presence of terms jhn or quadratic and cubic (positive 
coefficient) terms in y^. 

• Some right hand sides of ()4.76|) contain terms singular at the origin, such that their 
sums remain regular. Starting the integration at r = 0, is therefore impossible. To 
overcome this problem the solution in terms of the power series (|4.49j) - (j4.52j) was 
used within the interval [0, e] (for e = 0.01). 

• The solutions are extremely sensitive to initial conditions which made it unavoid- 
able to pass from single precision to double precision (from about 7 to about 15 
significant digits). However this didn't completely solve the problem. Even giving 
initial conditions (corresponding to a gravity localizing solution) at the origin to 
machine precision is in general not sufficient to obtain satisfactory precision at r max 
in a single Runge-Kutta forward integration step from e to r max . 

4.4.2 Discussion of solutions - fine-tuning surface 

The method used for exploring the parameter space (a, /3, 7) for gravity localizing so- 
lutions was a generalized version of the shooting method called the multiple shooting 
method jl6H ll62tll63|H64*j . In the shooting method, a boundary value problem is solved 
by combining a root-finding method (e.g. Newton's method |161p with forward integra- 
tion. In order to start the integration at one boundary, the root finding routine specifies 
particular values for the so-called shooting parameters and compares the results of the 
integration with the boundary conditions at the other boundary. This method obviously 
fails whenever the "initial guess" for the shooting parameters is too far from a solution 
such that the forward integration does not reach the second boundary. For this reason, 
the multiple shooting method was used in which the interval [e , r max ] was divided into a 
variable number of subintervals in each of which the shooting method was applied. This 
of course drastically increased the number of shooting parameters and as a result the 
amount of computing time. However, it resolved two problems: 

• Since the shooting parameters were specified in all subintervals the precision of the 
parameter values at the origin was no longer crucial for obtaining high precision 
solutions. 

• Out of range errors can be avoided by increasing the number of shooting intervals 
(again at the cost of increasing computing time). 
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Figure 4.1: Fine-tuning surface for solutions with the metric asymptotic The 
bold line separates solutions that localize gravity (c > 0) from those that do not (c < 0). 
Numerically obtained values 0/7 = -^h are plotted as a function of a = \ and (3 = . 



Despite all these advantages of multiple shooting, a first combination of parameter values 
(a, j3, 7, J'(0) , K"(0)) leading to gravity localization could not be found by this method, 
since convergence depends strongly on how close initial shooting parameters are to a 
real solution. This first solution was found by backward integration combined with 
the simplex method for finding zeros of one real function of several real variables. For 
a discussion of the simplex method see e.g. |161j . This solution corresponds to the 
parameter values (a = 1.429965428, (3 = 3.276535576, 7 = 0.025269415). Once this 
solution was known, it was straightforward to investigate with the multiple shooting 
method, which subset of the (a, /3, 7)-space leads to the desired metric asymptotic. We 
used a known solution (cti, (3\, 71) to obtain starting values for the shooting parameters 
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of a closeby other solution (a,\ + 5a, Pi + 5(3, 71 + 5j). This leads in general to rapid 
convergence of the Newton-method. Nevertheless, 5a, 5(3, £7 still had to be small. By 
this simple but time-consuming operation the so-called fine-tuning surface in parameter 
space was found, see Fig. 14.11 A point on this surface (ao, Po, 70) corresponds to a 
particular solution with the metric asymptotic ()4.59|) . The bold line separates gravity 
localizing solutions c > from solutions that do not localize gravity c < 0. The parameter 
space has been thoroughly explored within the rectangles Aa x A/3 = [0, 10] x [0.5, 4.0] 
and Aa x A/3 = [0, 1] x [3.5, 10] in the (a, /3)-plane. To illustrate the dependence on 
the parameter (3 (strength of gravity), we present in the following subsection a series of 
solutions corresponding to decreasing values of (3 for a fixed value of a = 1. 

4.4.3 Examples of numerical solutions 

Fig. 14.21 and Fig. 14.51 show two different numerical solutions corresponding to positive 
and negative bulk cosmo logical constant, respectively. Both solutions localize gravity 
c > 0. The corresponding parameter values are (a = 1.0000000, {3 = 5.50000000, 7 = 
-0.05434431) and (a = 1.0000000, j3 = 3.50000000, 7 = 0.02678351). Figs. IP and IP 
and Figs. 14.41 and 14.71 show the corresponding components of the stress-energy tensor 
and the curvature invariants. For high values of (3 the metric function C(r) develops a 
maximum before attaining its boundary value. Gravity dominates and the volume of the 
transverse space stays finite. Furthermore, it can be seen from Fig. 14.11 that for every 
fixed a there is a particular value of (3 such that c equals zero, which is the case for all 
points on the solid line shown in Fig. 14.11 By looking at eq. (|4.66|) . we immediately see 
that c = is equivalent to (3 2/ y = 2^10 = §■ We will discuss the limit c 1 in more 
detail in section 14.61 For lower values of (3 the transverse space has infinite volume and 
gravity can not be localized. See Figs. 14.81 and 14.91 for a solution that does not localize 
gravity (c < 0) and corresponds to (a = 1.0000000, (3 = 1.80000000, 7 = 0.04053600). 
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Figure 4.2: Gravity-localizing solution with negative bulk cosmological constant corre- 
sponding to the parameter values (a = 1.00000000, = 5.50000000, 7 = -0.05434431) 
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Figure 4.3: Stress- energy components for the solution given in Fig. 
(a = 1.00000000, (3 = 5.50000000, 7 = -0.05434431) 
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Figure 4.5: Gravity-localizing solution with positive bulk cosmological constant corre- 
sponding to the parameter values (a = 1.00000000, (3 = 3.50000000, 7 = 0.02678351) 
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Figure 4.6: Stress-energy components for the solution given in Fig. \4-5\ 
(a = 1.00000000, = 3.50000000, 7 = 0.02678351) 
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Figure 4.7: Curvature invariants for the solution given in Fig. |^.<5| 
(a = 1.00000000, (3 = 3.50000000, 7 = 0.02678351) 
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Figure 4.8: Example of a solution that does not localize gravity (c < 0) corresponding to 
the parameter values (a = 1.00000000, (3 = 1.80000000, 7 = 0.04053600) 
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Figure 4.9: Stress-energy components for the solution given in Fig. \4-8[ 
(a = 1.00000000, (3 = 1.80000000, 7 = 0.04053600) 
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4.5 Prasad- Sommerfield limit (a = 0) 



The Prasad- Sommerfield limit (a = 0) was explored numerically for /3-values ranging 
from 0.4 to about 70. The corresponding intersection of the fine-tuning surface Fig. 14.11 
and the plane a = is given in Fig. 14.101 There exist no gravity localizing solutions 
as the separating line in Fig. 14.11 indicates. The point in Fig. 14.101 corresponds to the 
solution shown in Figs. 14.111 and 14.121 One sees that 7 tends to zero for (5 going to 
infinity and that 7 tends to —00 for (3 going to zero. 
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Figure 4.10: Section of the fine-tuning surface Fig. \4-l\ corresponding to the Prasad- 
Sommerfield limit a = 0. None of the shown combinations of (3 and ^-values correspond 
to solutions that lead to a finite transverse volume. The point indicates the sample 
solution given below. 
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Figure 4.11: Sample solution for the Prasad-Sommerfield limit corresponding to the pa- 
rameter values (a = 0.00000000, (3 = 3.60000000, 7 = 0.02040333) 
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Figure 4.12: Stress-energy components for the solution given in Fig. \4-ll 
(a = 0.00000000, (3 = 3.60000000, 7 = 0.02040333) 
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4.6 Fine-monopole limit - physical requirements on 
solutions 



The fine monopole can be characterized by c < 1. As anticipated in section 14.4.31 
c = in eq. ()4.66|) immediately leads to (3 2/ y = 2 ^w — §• Hence, we deduce that 
the fine-monopole limit can not be realized for a negative bulk cosmo logical constant. 
In addition, from (|4.fi7|) it follows that C = ^J~(3. The c <C 1 limit is qualitatively 
different from its analogue in the 6-D-string case j!4(Jj . The solutions do not correspond 
to strictly local defects. Einstein equations never decouple from the field equations. Due 
to the particular metric asymptotic (|4.59|h the stress-energy tensor components tend 
to constants at infinity in transverse space. In the 6-dimensional string case the fine- 
string limit was realized as a "strictly local" defect (in the language of chapter OJ) having 
stress-energy vanishing exponentially outside the string core. Despite this difference, 
the discussion of the physical requirements which we have to impose on our solutions is 
very similar. In the following we show that in the fine-monopole limit the dimensionful 
parameters of the system A7, M7, raw, A, e can be chosen in such a way that all of the 
following physical requirements are simultaneously satisfied: 

1. M 2 P equals (1.22 ■ 10 19 GeV) 2 . 

2. The corrections to Newtons law do not contradict latest measurements. 

3. Classical gravity is applicable in the bulk. 

4. Classical gravity is applicable in the monopole core (r = 0). 

To find solutions with the above mentioned properties, it is possible to restrict oneself 
to a particular value of a, e.g. a = \. This choice corresponds to equal vector and 
Higgs masses my/ = run- Even though extra dimensions are infinite, the fact that 
M(r) decreases exponentially permits the definition of an effective "size" vq of the extra 
dimensions: 

M = M e _ 2 r = M e 2 m w r = . (4.80) 

cm w 

In order to solve the hierarchy problem in similar lines to [2H we parametrize the 
fundamental gravity scale as follows: 

M 7 = Kl0 3 GeV. (4.81) 

k = 1 then sets the fundamental scale equal to the electroweak scale. 
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1. The expression (J4.48|) for the square of the Planck mass Mp can be approximated 
in the fine-monopole limit by using the asymptotic (|4.59|) for the metric in the 
integral (|4.48jl rather than the exact (numerical) solutions. This gives 

m 6 w c 
By using one of Einstein's equations at infinity 

Cl = - 1 - , (4.83) 

4c 2 + 2/3 7 ' v ' 

and by developing to lowest order in c one finds 



„ 4vrM 7 5 M n 2 1/1 , 

M ^^rwAc +0{c)] - (484) 



Numerical solutions for c — > and a = \ converge to the following approximate 
parameter values 

P = 3.266281, 
7 = 0.0468665, 
(3 2 1 = 0.4999995, 

M = 0.959721. (4.85) 

The above values were obtained by extrapolation of solutions to c = 0, see Figs. 14. 131 
to 14.161 Therefore, the errors are of order 10~ 6 which is considerably higher than 
average errors from the integration which were smaller than 10 -8 . 

Note that (3 2r j tends to | in the fine-monopole limit (see Fig. 14.15)) . Neglecting all 
orders different from i in the above expression for Mj, and using 1)4.80)1 one has 

5/2 

m w = 1.13 ■ 10~ 5 ^yjGeV , (4.86) 
where r has been parametrized by r = 0,2 ^ mm ~ 10 12 GeV -1 • |. 



2. Since Newton's law is established down to 0.2 mm j22j, we simply need to have 
£>1. 
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3. In order to have classical gravity applicable in the bulk, we require that curvature 
at infinity is negligible with respect to the corresponding power of the fundamental 
scale. By looking at eqs. ()4.62|) - (j4.65|) . we see that we have to impose 



m 



2 

c 2 m 2 w ^M 2 and < M 2 . (4.87) 

Using ()4.83j) and the first of the above relations, the second one can immediately 
be transformed into 

/? 7 m^<M 7 2 . (4.88) 

One then finds 

r s 3^ 1/2 

-15 S „ , „„ j , 1n -3 I K 



2 ■ W~ ^ < 1 and 4 ■ 10~ d j < 1 , (4.89) 

where (3 and 7 have again been replaced by their limiting values ()4.85|) for c — > 0. 

4. Classical gravity is applicable in the monopole core whenever the curvature invari- 
ants R 2 , RabR AB , RabcdR ABCD and CabcdCabcd are small compared to the 
forth power of the fundamental gravity scale. Since these quantities are of the 
order of the mass 77% (see (|4.55J1 - (|4.58J1 ) we have 

m 4 w < M 7 4 . (4.90) 

Using (gSD and (gSBt one finds 

' 3\ 1/2 
-8 ' K 



1.1 • 10 s ( — ) <i. (4.91; 

Finally, the fine-monopole limit requires 



2 



r m w 

Again ()4.86|) implies 



< 1 . (4.92) 



£3 \ 1/2 

1.8 - 10 7 ( ^- <1. (4.93) 



K 5 y 

It is now easy to see that for a wide range of parameter combinations all these require- 
ments on £ and k can simultaneously be satisfied. One possible choice is £ = 100 and 
K — 1. This shows that already in the case a = | there are physical solutions correspond- 
ing to a fine monopole in the sense of eq. ()4.92|) . respecting all of the above requirements 
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4.7 Conclusions 



As we have seen in this chapter, it is possible to realize the idea of warped compactifi- 
cation on a topological defect in a higher-dimensional space-time with n = 3 transverse 
dimensions by considering a specific field theoretical model. Numerical solutions were 
found in the case of a 't Hooft-Polyakov monopole residing in the 3-dimensional trans- 
verse space of a 7-dimensional bulk space-time. As discussed in the general case in 
section 1331 of chapter El the transverse space approaches a constant curvature space at 
spatial infinity, however, not necessarily an anti-de-Sitter space. Both signs of the bulk 
cosmological constant are possible in order to localize gravity. We also considered a 
fine-monopole limit in the case a = 1/2 {my/ = ran) and verified explicitly that the 
model proposed is not in conflict with Newton's law, that it leads to a possible solution 
of the hierarchy problem and that classical gravity is applicable in the bulk and in the 
core of the defect. 

In the next and last chapter of this work, we are going to study a brane-world scenario 
which in many respects resembles the Randall-Sundrum II model. We are going to treat 
the case of a singular 3-brane embedded in a 5-dimensional anti-de-Sitter bulk space-time. 
Contrary to the Randall-Sundrum II model, however, we consider a brane with induced 
metric given by Einstein's static universe rather than by 4-dimensional Minkowski space- 
time. One of the advantages of such a setup, compared to the Randall-Sundrum II case, 
is that the corresponding space-time manifold turns out to be geodesically complete. 
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Chapter 5 



Einstein static universe as a brane 
in extra dimensions 

5.1 Introduction 

Recent proposals of large [2HI EH] and infinite ^11231101 extra dimensions provide new 
opportunities for addressing several outstanding problems of modern theoretical physics 
like for example the hierarchy problem j2Sj _ |ZZ| or the cosmo logical constant problem 
[T3]-[THj. These proposals triggered an immense research activity in theories involving 
3-branes with interests ranging from elucidating the global space-time structure of brane 
world scenarios, properties of gravity, cosmology and brane cosmological perturbations, 
generalizations to higher dimensions etc. While the possibilities are rich, realistic sce- 
narios remain rare. For example, the simple Randall-Sundrum II model jlU] faces the 
problem of being geodesically incomplete |lfci5t I16b| llbTj . It is therefore reasonable to 
look for adequate alternatives or generalizations to the Randall-Sundrum II model which 
avoid amongst others the above mentioned problem while sharing its pleasant feature of 
an effective 4-dimensional low energy gravity on the brane. 

In this final chapter, we present a 5-dimensional brane-world model which solves 
the geodesic incompleteness of the Randall-Sundrum II model while preserving its phe- 
nomenological properties concerning the localization of gravity. 

We try to illustrate our motivation for considering a particular geometry by using the 
simple picture of a domain structure in extra dimensions, resulting from a spontaneously 
broken discrete symmetry. The associated Higgs-field takes different values in regions 
separated by domain walls, which restricts the possibilities of combining domain walls 
depending on the global topology of the space-time under consideration. For simplicity 



75 



we restrict ourselves to five dimensions. Let us analyze a few simple cases: 

• Non-compact extra dimension: y £ (—00, +00). In this case, by choosing the origin 
y = to coincide with the position of the brane, we obtain two non-overlapping 
regions (—00, 0) and (0, +00) and we can thus consistently have one brane in such 
a theory. An example for this configuration is provided by the Randall-Sundrum 
II model gO]. 

• Compact extra dimension: y £ [0,27r]. We now consider two possibilities, depend- 
ing on the spatial topology of our manifold: 

— (a) I 3 x S 1 : If the ordinary dimensions are supposed to be non-compact, it 
is not possible to consistently put only one brane in the extra dimensions. At 
least two branes are needed. 

— (b) S A : If all spatial coordinates are part of a compact manifold, the simple 
picture shown in Fig. 15.11 seems to suggest that it is possible to consistently 
put a single brane in the bulk space-time. 

To the best of our knowledge, case (b) has not yet been considered in the literature and 
our aim is to present such a construction. 

This chapter is organized as follows: in Section 15.21 we discuss the basic geometric 
and topological properties of our brane-world scenario like Einstein equations, junction 
conditions and the distance hierarchy between the extra dimension R and the observable 
universe Rjj. Section 1531 is dedicated to the study of geodesies and the demonstration, 
that our model does not suffer from being geodesically incomplete. In section 15.41 we 
present a detailed computation of the propagator of a massless scalar field in the given 
background, serving as a easy, phenomenological approach to the study of gravity. We 
discuss the behavior of the two-point function in three different distance regimes. It 
turns out that the computations are rather technical and we therefore collect large parts 
of it in four appendices. We eventually draw conclusions in section 1531 

5.2 The background equations 

The aim of this section is to present the topology and geometry of the brane-world 
model which we set out to study but also to motivate how this particular model emerged 
through imposing several physical conditions on the more general set of solutions. 
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Upper Hemisphere 




Figure 5.1: 5* 3 -brane embedded in S 4 . +1 and —1 symbolically represent the different 
vacuum expectation values taken by the Higgs field in different regions of space. 

5.2.1 Einstein equations 

We would like to construct a space-time with all spatial dimensions being part of the same 
compact manifold. One more motivation for this lies in a possible solution of the strong 
CP-problem within theories with extra dimensions The overall topology would then 
be given by R x K, where R represents the time-coordinate and K any compact manifold. 
For simplicity we will restrict ourselves to the case of one extra dimension and an induced 
metric on the brane characterized by a spatial component of geometry S 3 . The idea is to 
combine two 5-dimensional regions of space-time dominated by a cosmological constant 
A in such a way, that the border of the two regions can be identified with a 3-brane, 
constituting our observable universe. As pointed out in [21] the manifold K has to be 
highly anisotropic in order to single out the small extra dimension from the three usual 
ones. It is not a priory clear whether the Einstein equations allow for such solutions at 
all and if so whether gravity can be localized on the brane in such a setup. 

We choose the following ansatz for the 5-dimensional metric consistent with the above 
requirements: 

ds 2 = g MN dx M dx N = -a 2 (6) dt 2 + R 2 V 7 2 (0) dQ 2 3 + R 2 d6 2 , (5.1) 

where 6 G [— | , | ] denotes the extra coordinate and R and Ru are constants represent- 
ing the size of the extra dimension and the size of the observable universe at present, 
respectively, dfl 2 denotes the line element of a 3-sphere: 

dQ 2 = dfl + sin 2 ifi dip\ + sin 2 if\ sin 2 <p 2 dtp\ , (5-2) 

where (fi, (p 2 belong to the interval [0, ir], <p 3 to [0, 2rr]. Capital latin letters M,N, .. will 
range from to 4. In order to obtain a compact space, we require 7(±|) = 0. 
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Metrics which can locally be put into the form (J5.1|) are sometimes referred to as 
asymmetrically warped metrics due to two different functions (a and 7 in our nota- 
tions) multiplying the temporal and spatial part of the 4-dimensional coordinate dif- 
ferentials. Their relevance in connection with 4-dimensional Lorentz-violation at high 
energies was first pointed out in |118j . More recent discussions of this subject can be 
found in |ltj8j . [TBU] or [165j and references therein. A common prediction of theories of 
this kind is that dispersion relations get modified at high energies. For a field-theoretical 
discussion of Lorentz violating effects in the context of the Standard Model of particle 
physics see |17Uj . 

The Einstein equations in 5 dimensions with a bulk cosmological constant A and a 
stress-energy tensor T MN take the following form: 



Rmn — 2 ^ 9m n + A gMN 



87T 



T, 



MN 



(5.3) 



where M is the fundamental scale of gravity. We choose to parametrize the stress-energy 
tensor in the following way, consistent with the symmetries of the metric: 



T° = eo (£)<0, T 1 =5)e{9) ) T% = e e (9) 



(5.4) 



where, as indicated, the above diagonal components depend only on the extra coordi- 
nate 9. Lower case latin indices i,j label the coordinates on S 3 . 

Using the metric ansatz (j5.1j) together with the stress-energy tensor (J5.4J) the Einstein 
equations ()5.3|) become 



1 

#2 



3 



+ 



y a 1 y 
+ 2- — + 2— 



7 a 



7 



7 
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a 



3 



7 



7 a 
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Ru 



r 



v 



r 



A 



A 



A 



8tt 



87T 

M3 e ' 



87T 



(5.5) 
(5.6) 
(5.7) 



where ' denotes differentiation with respect to 9. The conservation of stress-energy, or 
equivalently, the Bianchi-identities lead to the following constraint relating the three 
independent components eo, e and eq: 



e« + I — + 3— 

a 7 



o Q 7 
— e - 3 — e 

a 7 



0. 



(5. 
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5.2.2 Vacuum solution 



We do not intend to provide a field theoretical model which could generate the geometry 
we are about the describe. Our aim is to study a singular brane, located at 9b = 
separating two vacuum regions in the bulk. Of course the solutions to the Einstein 
equations in vacuum with a cosmological constant term are nothing but the familiar 
de-Sitter (A > 0), Minkowski (A = 0) and anti-de-Sitter (A < 0) space-times. To see 
how these geometries can be recovered using the line element (J5.1)) we set the right hand 
sides of (|5.5j) - ([5.7)) equal to zero and subtract eq. ()5.7|) from eq. ()5.5|) to obtain a = cy', 
where c is a constant. Putting this back in eqs. ()5.5|) - (j5.7|) . we are left with only one 
independent differential equation 



3 

R 2 



+ 



r 

7 



R 

Ru 



r 



A = 



(5.9) 



The solution of ()5.9|) is trivial and by means of a = cy' we find 



sinhKi |,|)] coshK! |,|)] 

smh(o;|) cosh(u;|) 

provided that Ru = Rsmhiui^) /oj and oj 2 = —AR 2 /6. Notice that the solution (j5.1()j) 
does not contain any integration constant because we already imposed the boundary 
conditions 7(±f ) = and 7(0) = <r(0) = 1. Moreover, we chose the whole setup to be 
symmetric under the transformation 9 — > —9. From (j5.10|) it is now obvious that locally 
(for 9 > and 9 < 0) the line element (J5.1j) correctly describes de-Sitter, Minkowski and 
anti-de-Sitter space-times for imaginary, zero and real values of u, respectively. For our 
purposes only the AdS solution will be of any interest as we will see in the following 
section. A simple change of coordinates starting from ()5.1|) and ()5.10j) shows that the 
AdS-i&dms in our notations is given by RacIS = R/uj = a/— 6/A. 1 At this stage, the 
validity of the vacuum solution (|5.1U|) is restricted to the bulk, since it is not even 
differentiable in the classical sense at 9 = 0. In order to give sense to (J5.1U)) for all values 
of 9 we will have to allow for some singular distribution of stress-energy at 9 = and 
solve ()5.5|) - (j5.7|) in the sense of distributions. 



1 See also footnote El 
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5.2.3 The complete solution for a singular brane 

As announced, we now refine our ansatz for the energy-momentum tensor (|5.4|) to allow 
for a solution of eqs. (|5.5|) - (|5.7J1 in the whole interval — 7r/2 < 9 < n/2: 

m = »'m, m=e m, £ , w= c,f). (5.H) 

After replacing the above components (|5.11j) in eqs. (|5.5|l - (|5.7|) and integrating over 9 
from —i] to T], followed by the limit i] — > we find: 

-■■ =ZM^ (5 ' l2) 

:^: + 2-X: = (5,3) 

= c e , (5.14) 

where the symbol : . . . : is used to denote the jump of a quantity across the brane defined 
by: 

:/: Elimf/W-ZH)]. (5.15) 
Note that in the above step we made use of the identity 

" ^ ' ^ (5.16) 



7 \7/ V7 

together with the continuity of 7 on the brane. Specifying (j5.l2j) and (|5.l3j) to (J5.l0|) we 

have: 

c = -|-M 3 | coth (^) , (5.17) 



— M — 
4tt R L 



tanh ( ijj— ) + 2 coth \ uj — 



(5.18) 

c e = 0. (5.19) 



Eqs. (I5.l7jl and (|5.l8jl relate the energy-density and the pressure of the singular brane 
to the bulk cosmological constant A (via lu), the size of the extra dimension R and 
the fundamental scale of gravity M. With the above relations (j5.l7|) - (j5.l9|) we can 
now interpret (j5.l0|) as a solution to the Einstein equations ()5.5|) - ([5.7)) in the sense of 
distributions. Also the stress-energy conservation constraint ()5.8|) is satisfied based on 
the identity S(x) sign(x) = again in the distributional sense. 
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Note that Co/c — > 1 in the limit u — > oo. Moreover, for larger and larger values of 
u), Co and c approach the brane tension of the Randall- Sundrum II model and the above 
eqs. (15. 17(1 and (j5.18j) merge to the equivalent relation in the Randall-Sundrum II case. 
This is no surprise since taking the limit u — > oo corresponds to inflating and flattening 
the 3-brane so that we expect to recover the case of the flat Randall-Sundrum II brane. 

We finish this section by the discussion of some physical properties of our manifold. 
We first observe that its spatial part is homeomorphic to a 4-sphere S 4 . This is obvious 
from the metric ()5.1|) and the explicit expression for 7 given in ([5. 10)1 . Geometrically, 
however, our manifold differs from S* 4 due to the high anisotropy related to the smallness 
of the extra dimension. The ratio of typical distance scales in the bulk and on the brane 
is given by 

— = — %r • (5-20) 
Ru sinh(^f) V ; 

It is now clear that the above ratio (|5.2Uj) can only be made very small in the case of real 
uj (Ac/S-space-time). For the size of the observable universe we take the lower bound 
Ru > 4Gpc ~ 10 28 cm while the size of the extra dimension is limited from above : 
R < 10~ 2 cm, leaving us with to > 50. 

Finally we would like to point out a very peculiar property of the manifold under 
consideration: as it can immediately be deduced from the line element (J5.1|) and (|5.1U|) . 
any two points on the brane are separated by not more than a distance of the order of R, 
regardless of their distance as measured by an observer on the brane using the induced 
metric. 



5.3 Geodesies 

It is well known, that the Randall-Sundrum II model is time-like and light-like geodesi- 
cally incomplete |lfi5| I167j which means that there exist inextendible time-like and 
light-like geodesies. 2 An inextendible geodesic is a geodesic parametrized by an affine 
parameter r such that by using up only a finite amount of affine parameter the geodesic 
extends over infinite coordinate distances. In a more physical language one could re- 
formulate the above statement by saying that it takes only a finite amount of affine 
parameter r in order to reach the infinities of the incomplete space-time. As we will il- 
lustrate later in this chapter, the reason why the Randall-Sundrum II setup ceases to be 
geodesically complete is simply due to a specific way of gluing two patches of AdS 5 . One 

2 For a precise definition of a geodesically incomplete space see e.g. |171| . 



81 



of the main motivations for considering this setup was to provide an alternative to the 
Randall-Sundrum II model that has the advantage of being geodesically complete while 
conserving the pleasant phenomenological features of the latter. We divide the discus- 
sion of geodesies in two parts: in subsection 15.3.11 we illustrate the effects of incomplete 
geodesies in the Randall-Sundrum II setup for time-like and light-like geodesies. In the 
following subsection 15.3.21 we demonstrate why our setup is geodesically complete by 
looking at corresponding geodesies. Finally, we complement the discussions by illustrat- 
ing the physics with the use of the Penrose-diagram of (the universal covering space-time 
of) AdS 5 . 

5.3.1 Geodesies in the Randall-Sundrum II setup 

Our discussion of geodesies in this chapter is in no sense meant to be complete. Without 
going into the details of the computations we merely intend to present the solutions of 
the geodesic equations in certain cases. For more general and more complete discussions 
of this issue we refer to the literature, see e.g. 166, 172J and references therein. We 
first consider light-like geodesies in the Randall-Sundrum II background metric given in 
appendix O eq. (|C.1|) . Let us suppose that a photon is emitted at the brane at y — 
in the positive ^/-direction at coordinate time t = then reflected at y = y\ at the time 
t = t\ to be observed by an observer on the brane at time t = £2- In this situation t 
corresponds to the proper time of an observer on the brane at rest. A simple calculation 
reveals 



A brane bound observer will therefore note that it takes an infinite time for a photon 
to escape to y = 00. However, parameterizing the same geodesic by an affine parameter 
r reveals the light-like incompleteness of the Randall-Sundrum II space-time. Let the 
events of emission, reflection (at y = y\) and arrival on the brane again be labeled by 
t = 0, r = 7~i and r = r 2 , respectively. Using the geodesic equation 



£2 = 2ti = - 





d 2 x^ 



1 Ui 



dx v dx p 




dr 2 



vp dr dr 







an easy computation shows: 



r 2 = 2n = — 




(5.23) 
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The constant c is a remnant of the freedom in the choice of an afhne parameter. 3 From 
the last equation we see that now 

lim n = \, (5.24) 

j/i^oo ck 

meaning that in order to reach infinity (y = oo) in the extra dimension it takes only 
a finite amount l/(ck) of affine parameter r, the expression of incompleteness of the 
Randall-Sundrum II space-time with respect to affinely parametrized light-like geodesies. 

In the case of time-like geodesies, the inconsistency is even more striking. As shown in 
e.g. |lb7| , a massive particle starting at the brane with vanishing initial velocity travels 
to y = oo in finite proper time given by r p = n/(2k) while for a brane-bound observer 
this happens in an infinite coordinate time. To summarize, the Randall-Sundrum II 
brane-world model is geodesically incomplete both for null and for time-like geodesies. 
As we will see in the next section, the geodesic incompleteness is a direct consequence 
of the use of a particular coordinate system in AclSc,, the so-called Poincare coordinate 
system, which covers only a part of the full AdS§ space-time. We will also see that the 
problem of incomplete geodesies is absent in the setup we propose here. 



5.3.2 Geodesies in the background ( 15.11) and Penrose-diagram 

We would now like to answer similar questions to the ones considered in the previous 
section for the background (|5.1|) . For example, we would like to know what time it takes 
for light to travel from the brane (at 8=0) in the ^-direction to a given point in the 
"upper hemisphere" with 9 coordinate Q\, to be reflected and to return to the brane. As 
in the last section, t and r will denote the coordinate time (proper time of a stationary 
observer on the brane) and the affine parameter used for parameterizing the geodesies, 
respectively. Again we choose t = (r = 0) for the moment of emission, t\ fa) for the 
reflection and t 2 fa) for the time where the photon returns to the brane. Due to the 
enormous hierarchy of distance scales in our model, one might wonder whether photons 
(or gravitons) are able to carry information from an arbitrary point on the brane to 
any other point on the brane connected to the first one by a null geodesic in the extra 
dimension. For an observer on the brane such a possibility would be interpreted as 4- 
dimensional causality violation. However, as we will see in the following, none of these 

3 In general, two affine parameters A and r are related by A = ct + d, since this is the most general 
transformation leaving the geodesic equation (|5.22(l invariant. The choice of the origin of time to coincide 
with the emission of the photon only fixes d but does not restrict c. 
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possibilities exist in our model. Omitting all details we find 



t 2 = 2ti = 4R V coth 



— arctan 
V 2 / 



sinh (^) 



cosh [f (tt - 9 l 



so that an observer on the brane will see that the photon reaches the 
9 1 = 7i / 2 at finite time 



(5.25) 



'north pole" 



ti = Rjj coth J arctan 



sinh ( — ) 
V 2 J 



TT 



Rij— 



(5.26) 



Note that due to the warped geometry, it is Rjj entering the last relation and not R, so 
that even though the physical distance to the "north pole" is of the order of R it takes 
a time of the order of Rjj for photons to reach it, excluding causality violation on the 
brane as discussed above. If the same geodesic is parametrized using an affine parameter 
we obtain 



r 2 = 2ti 



2R 



c to 



i f UJ7l \ 
tanh — 

V 2 / 



sinh [u (| - X )] 



cosh^f 



(5.27) 

tt/2 starting from the 
(5.28) 



such that the amount of affine parameter needed to reach 9\ 
brane is: 

1R , /Co>7T\ IR 

At = tanh — a; . 

c uj V 2 / ecu 

Here again c reflects the freedom in the choice of the affine parameter. 

The results formally resemble those of the Randall- Sundrum II case. However, the 
important difference is that in our case each geodesic can trivially be extended to ar- 
bitrary values of the affine parameter, a simple consequence of the compactness of our 
space. Once the photon reaches the point 9 = it continues on its geodesic, approaching 
the brane, entering the "southern hemisphere", etc. It is clear that it needs an infinite 
amount of affine parameter in order to travel infinite coordinate distances. Therefore, 
the null geodesies in our setup which are the analogues of the incomplete null geodesies in 
the Randall-Sundrum II setup turn out to be perfectly complete due to the compactness 
of our space. The situation for time-like geodesies is fully analog to the case of the null 
geodesies. 

To end this section about the geometric properties of our model, we would like to 
discuss the conformal structure of our space-time and point out differences to the Randall- 
Sundrum II setup. Let us review briefly the basic properties of AdS 5 space-time to the 
extend that we will need it in the following discussion. 4 AdS 5 space-time can be thought 



4 We mainly follow |17rij . 
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of as the hyperboloid defined by 



X 2 + Xl - X\ - X 2 - X\ - X 2 = a 2 (5.29) 

embedded in a flat space with metric 

ds 2 = -dXl - dX 2 + dX\ + dX 2 + dXl + dX 2 , (5.30) 

a being the so-called AdS-radius. The global coordinates of AdS$ are defined by 

X = a cosh x cos r, X 5 =a cosh x sin r, 

Xi = a sinh x cos , X 2 =a sinh x sin (pi cos <£> 2 , 

X3 = a sinh x sin y9i sin cos ^3, X4 =a sinh x sin yji sin sin y93, (5.31) 

where the coordinates are confined by < x ? — 71 ^ r ^ t 1 "; < <^i < it, < (^2 < tt, 
< </?3 < 27r and r = — tt is identified with r = it. These coordinates cover the full 
hyperboloid exactly once. Allowing r to take values on the real line without the above 
identification of points gives the universal covering space CAdS$ of AdS$. 5 In these 
coordinates, the line element (J5.3U)) can be written: 6 

ds 2 = a 2 (- cosh 2 x dr 2 + d X 2 + sinh 2 x dttf) . (5.32) 

Another coordinate system can be defined by 

X = — [l + u 2 (a 2 + x 2 - i 2 )] , X 5 =a u i, 

X i = aux\ i = 1,2,3 , X 4 =— [l - u 2 (a 2 - x 2 + t 2 )] , (5.33) 

with 11 > 0, t G (—00,00) and x % G (—00,00). In these Poincare coordinates the line 
element (15.30)) takes the form 



ds = a 



du 



2 



— + u 2 (-dt 2 + dx 2 ) 
u ' 



(5.34) 



In contrast to the global coordinates, the Poincare coordinates do not cover the whole 
of the AdS 5 and CAdS 5 space-times |173j . From (|5.34|) . after changing coordinates 

5 Whenever we used the word AdS^ so far in this chapter we actually meant CAdS§. For reasons of 
clarity, however, we will brake with this common practice in the rest of this section. 

6 It is obvious that the metric l|5.1|l for < 6 < n/2 (— ir/2 < 8 < 0) reduces to the above line-element 
(|5.32|l under the following coordinate transformation: \ = =F 8), r = tujj [i?cosh(^)], together 
with a — R/lu. 
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according to dy = —adu/u and rescaling t and x l by the AdS'-radius a, we recover the 
original Randall-Sundrum II coordinate system given in (|C.1|) . The restrictions on y in 
the Randall-Sundrum II setup further limit the range covered by their coordinate system 
to the < u < 1 domain of the Poincare coordinates. 

Coming back to the global coordinates, we introduce p by 

7T 

tanp = sinhx with < p < — , (5.35) 

so that ()5.32j) becomes 

ds 2 = i-dr 2 + dp 2 + sin 2 p dVll) . (5.36) 

cos^p v 

The Penrose- diagrams of AdS^ space-time and its universal covering space-time CAdS§ 
are shown in Fig. 15.21 see j!7 llTl 75 . While the AdS§ space-time contains closed time-like 
curves (denoted by 7 in the figure), its universal covering space-time CAdS$ does not. 
The arrows in the left diagram indicate that the lines t = —71 and t = it should be 
identified. The symbol J* stands for the time-like surface p = tt/2 (spatial infinity). It is 
this surface which is responsible for the absence of a Cauchy surface in AdS-spa.ce. We 
will concentrate in the following on the CAdS^ diagram. First we note that each point in 
the diagram corresponds to a 3-sphere. The shaded region indicates the patch covered by 
the Poincare (and Randall-Sundrum) coordinates. Note that the position of the Randall- 
Sundrum brane cannot be represented in a simple way in the Penrose diagram of CAdS$. 
The reason is that the u = const, hypersurfaces of the Poincare coordinates generate 
a slicing of flat 4-dimensional Minkowski space-times while the points in the diagram 
represent (curved) 3-spheres. From the Penrose diagram it is immediately clear that the 
Randall-Sundrum II space-time is geodesically incomplete. The time-like curves denoted 
by 7 emanating from the origin (p = 0, r = 0) will all eventually exit the shaded region 
after a finite coordinate time r. These geodesies appear inextendible from the point 
of view of the Randall-Sundrum space-time. The problems arise due to the arbitrary 
cutting of a space-time along the borders of a given coordinate patch which covers only 
a part of the initial space-time. Similar conclusions can be drawn for null geodesies, 
represented by lines making angles of 45 degrees with the vertical lines in Fig. 15.21 The 
vertical dotted line with topology R x S 3 at coordinate p^ = arctan [sinh f^)] close to 
p = tt/2 in the right diagram corresponds to the location of our curved 3-brane, r being 
proportional to our time coordinate t (see footnote • We note that since u> > 50, this 
line is drawn by far too distant from p = it/2 as a simple expansion shows: 

--2e-- + 0(e~ — ). (5.37) 



Pb = arctan 
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P = ■§ p = Pi,p=§ 

((pi = 7r) (ip! = 0) 



Figure 5.2: Penrose diagram of space and its universal covering space CAdS^ 

The slice of CAdS$ to the right of pb (pb < p < 7r /2) is discarded and replaced by another 
copy of the slice to the left of pb (0 < p < pb) ■ 

From the Penrose diagram we can also deduce that our space-time is geodesically 
complete. Time-like geodesies (emanating from (p = 0, r = 0)) cross the brane pb at 
some point or return to the origin depending on the initial velocities of the particles that 
define them. In both cases (as in the case of null geodesies) there is no obstacle to extend 
the affine parameter to larger and larger values. 

As we already mentioned, the absence of a Cauchy surface in CAdS space-time is due 
to the existence of the time-like surface J* at spatial infinity. Even though the family 
of space-like hypersurfaces r = const, cover the whole CAdS space-time, none of these 
hypersurfaces is a Cauchy surface. By considering the two null-geodesics which cross 
the origin (r = p = 0) in the right diagram of Fig. 15.21 we note that they will never 
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reach for example the hypersurface r = n. Since the affine parameter A corresponding 
to these geodesies lies within the full interval (—00, 00), the hypersurface r = n can not 
be a Cauchy surface. Put in another way, there are points on J? which lie neither in the 
future nor in the past Cauchy development of the hypersurface r = n. Similar arguments 
hold for any other r = const, hypersurface. 

By construction, our space-time excludes J* . A corresponding null-geodesics in our 
setup which starts for example at the "north pole" 9 = n/2 will reach the brane after 
using up only a finite amount of affine parameter and can be extended to the "southern 
hemisphere" where it will eventually reach the "south pole" 9 = —n/2 (again by using up 
only a finite amount of affine parameter), etc. It is therefore clear that all inextendible, 
non space-like lines in our manifold eventually have to cross any given r = const, hyper- 
surface. We conclude that the family of r = const, hypersurfaces indeed constitutes a 
family of Cauchy surfaces for the manifold considered in this chapter. 

5.4 Gravity localization on the brane 

There are at least two equivalent ways of addressing the problem of the localization of 
gravity in brane-world scenarios. Both ways have their advantages and disadvantages. 
The first way is based on a detailed study of the Kaluza-Klein excitations of the graviton. 
After integrating out the extra dimension(s) one obtains an effective 4 dimensional La- 
grangian involving the full tower of Kaluza-Klein gravitons. Considering the low energy 
scattering process of two test particles on the brane via exchange of Kaluza-Klein exci- 
tations allows to relate the non-relativistic scattering amplitude to the static potential 
felt by the two test particles. Since the coupling of each individual Kaluza-Klein particle 
to matter on the brane is proportional to the value of its transverse wave-function at the 
position of the brane, this approach necessitates a proper normalization of all Kaluza- 
Klein modes. An advantage of this approach is the possibility of distinguishing between 
the contributions to the potential coming from the zero mode (Newton's law) and from 
the higher modes (corrections thereof). 

The second approach is based on a direct calculation of the graviton two-point func- 
tion in the space-time under consideration, having the obvious advantage of bypassing all 
technicalities related to the Kaluza-Klein spectrum and the wave-function normalization. 
However, a physical interpretation of the effects of individual Kaluza-Klein modes from 
the point of view of a 4-dimensional observer is, to say the least, not straightforward. 
Finally, due to the equivalence of the two approaches, it is clearly possible to fill this gap 
by reading of the Kaluza-Klein spectrum and the (normalized) wave-functions from the 
two-point function by locating its poles and determining corresponding residues. 
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In the setup considered in this chapter, we choose to work in the second approach, 
the direct evaluation of the Green's functions, due to additional difficulties arising from 
the non-Minkowskian nature of the induced metric on the brane. The lack of Poincare 
invariance on the brane clearly implies the absence of this symmetry in the effective 4- 
dimensional Lagrangian as well as non-Minkowskian dispersion relations for the Kaluza- 
Klein modes. 

5.4.1 Perturbation equations and junction conditions 

We would like to study the fluctuations H^ of the metric in the background (j5.1J) defined 
by: 

ds 2 = -o 2 dt 2 + R 2 vl 2 [rjij + 2K i3 \ d^d V 3 + R 2 d6 2 . (5.38) 

transform as a transverse, traceless second-rank tensor with respect to coordinate 
transformations on the maximally symmetric space with metric r]^, a 3-sphere in our 
case: 

rflHij = , V ij ViH jk = , (5.39) 

where V denotes the covariant derivative associated with the metric rjij on S 3 . The gauge 
invariant symmetric tensor Hij does not couple to the vector and scalar perturbations. Its 
perturbation equation in the bulk is obtained from the transverse, traceless component 
of the perturbed Einstein equations in 5 dimensions: 

5R MN -^A5g MN = 0. (5.40) 
Specifying this equation to the case of interest of a static background we find: 

h 6 - 4f AH " - h H « - h (7 + 3 ~i) H -> + 4f H -> = ' (5 - 41) 

where A denotes the Laplacian on S 3 . The last equation is valid in the bulk in the 
absence of sources. We can now conveniently expand Hy in the basis of the symmetric, 
transverse, traceless tensor harmonics on S* 3 [H3J I176| 1177] : 

00 

ff« = E£* <u) M)r«*\ ( 5 - 42 ) 

1=2 X 

where the T^- satisfy 

Atp + kffW = , kf = 1(1 + 2) - 2 , I = 2, 3, ... 

V ij Vif^ X) = , rflfM = , fj$ = . (5.43) 



N9 



The sum over A is symbolic and replaces all eigenvalues needed to describe the full 
degeneracy of the subspace of solutions for a given value of /. Introducing the expansion 
(|5.42|) into (j5.41j) and using the orthogonality relation of the tensor harmonics |177j 

jVvV tk rf l fpf^d^ = 5 ll '5^, (5.44) 



we obtain: 



u R 2 R 2 \a 7 J Rjjl 



(5.45) 



This equation has to be compared to the equation of a massless scalar field in the back- 
ground (j5.1j) . After expanding the massless scalar in the corresponding scalar harmonics 
on S 3 , we recover (J5.45|) with the only difference in the eigenvalue parameters due to 
different spectra of the Laplacian A for scalars and for tensors. Motivated by this last 
observation, and in order to avoid the technicalities related to the tensorial nature of the 
graviton Hij, we confine ourselves to the study of the Green's functions of a massless 
scalar field in the background (j5.1|) . 

The differential equation ()5.45|) alone does not determine uniquely. We have to 
impose proper boundary conditions for Hy. While imposing square integrability will 
constitute one boundary condition at 9 = ±tt/2, the behavior of on the brane will 
be dictated by the Israel junction condition 



•K^: =-^(r^-±T«g^y (5.46) 

Here K^ u and g^ v denote the extrinsic curvature and the induced metric on the brane, 
respectively. T^ v is the 4-dimensional stress-energy tensor on the brane 7 and the symbol 
: . . . : is defined in (I5.15J) . In our coordinate system the extrinsic curvature is given by: 



= (5 - 47) 



Hence, the non-trivial components of the Israel condition become: 



SttR f2 

-c - c , (5.48) 



' a ' M 3 \3 

V 87ri? 
: -L ( Vij + 2H l3 ) + E[ 3 : = ^c ( V ij + 2ffy) . (5.49) 



7 T° = cq, and T\ = 5 l c in the notation of section l5.2.3l 
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Separating the background from the fluctuation in (|5.49j) we find: 



8nR 



(5.50) 



' 7 ' 3M 3 



(5.51) 



where we only used the continuity of Hij on the brane. While eq. ()5.50j) directly coincides 
with eq. (|5.12|) of section 15.2.31 eq. ()5.48|) turns out to be a linear combination of ()5.12|) 
and (lo~T3j) . 

Eq. (|5.51|) taken alone implies the continuity of H'- on the brane. If in addition the 
fluctuations are supposed to satisfy the 6^—6 symmetry, this condition reduces to a 
Neumann condition, as for example in our treatment of the scalar two-point function in 
the Randall- Sundrum II case (see appendix IU|) . 

In the next section we are going to find the static Green's functions of a massless 
scalar field in the Einstein static universe background. This will serve as a preparatory 
step for how to handle the more complicated case of a non-invertible Laplacian. More 
importantly, it will provide the necessary reference needed for the interpretation of the 
effect of the extra dimension on the potential between two test masses on the brane. 

5.4.2 Gravity in the Einstein static universe 

As announced in the previous section, we will concentrate on the massless scalar field. 
Our aim is to solve the analog of Poisson equation in the Einstein static universe back- 
ground. Due to its topology R x S* 3 we will encounter a difficulty related to the existence 
of a zero eigenvalue of the scalar Laplacian on S 3 necessitating the introduction of a 
modified Green's function. 8 Note that in the more physical case where the full tensor 
structure of the graviton is maintained no such step is necessary since all eigenvalues of 
the tensorial Laplacian are strictly negative on S 3 , see (|5.43|) . 

Definition of a modified Green's function 

We take the metric of Einstein's static universe in the form: 



with dVt\ being the line element of a 3-sphere given in (|5.2j) and A being its constant 
radius. The equation of a massless scalar field then becomes: 



ds 2 = g^dx^dx" = -dt 2 + A 2 dQ 2 3 



(5.52) 




(5.53) 



For an elementary introduction to the concept of modified Green's functions see e.g |178| . 
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which in the static case reduces to 

Vu(x) = -^Am(x) = j(x) , (5.54) 
A 2 

where A is the scalar Laplacian on S 3 and we choose to write x for the collection of 
the three angles on S 3 . All functions involved are supposed to obey periodic boundary 
condition on S 3 so that in Green's identity all boundary terms vanish: 



g {Vv (x)) u(x) d x = j \J —g v{x) V u(x) d x = j ^—gv(x)j(x) d x . (5.55) 

Here and in the following, bars denote complex conjugate quantities. The operator T> 
trivially allows for an eigenfunction with zero eigenvalue (the constant function on S 3 ): 

T>u (x) = , J ~gu (x) Uq(x) d 3 x = 1 . (5.56) 

Then, since the homogeneous equation has not only the trivial solution (u = 0) which 
satisfies the periodic boundary conditions in the angular variables, the operator T> cannot 
be invertible. Therefore, there is no solution to the equation ()5.54j) for an arbitrary 
source. In order to still define a "modified" Green's function we have to restrict the 
possible sources to sources that satisfy the following solvability condition: 

/—gj(x) uo(x) d 3 x = . (5.57) 

Formally this condition can be obtained by replacing v(x) by the non-trivial solution of 
the homogeneous equation Uq(x) in (|5.55|) . We now define the modified Green's function 
as a solution of 

V x Q(x, x') = 5 ^ X ^ - " Uo(x) %(£') • (5.58) 

Replacing v(x) by Q(x,x') in Green's identity ()5.55j) we obtain (after complex conjuga- 
tion) the desired integral representation for the solution of ()5.54j) : 

u(x) = Cu (x) + J ^/^gj(x') G{x\ x) d 3 £ , (5.59) 

where C is given by 

C = y V~ 9 u{x) uq{x) d 3 x . (5.60) 
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Note that the source 



/ — *\ fi (x x) /-*\- /-^ 



-g 



(5.61) 



trivially satisfies the solvability condition ()5.57|) and represents a point source at the lo- 
cation x = x' compensated by a uniform, negative mass density To fix the normalization 
of u (x) we write: 



(5.62) 



u (x) = N $ioo(x) with / <yrj$ 100 (x)$ioo(x)d x = 1, 



where $ioo = — and so 

7TV ^ 



-^u (f)n (f)(i 3 f = |iV |M 3 / v ^$ 10 o(x)$ 10 o(x)^= \N \ 2 A S = 1, (5.63) 



so that 



^3/2 



$100 (5). 



(5.64) 



In order to solve the differential equation ()5.58|) defining the modified Green's function 
Q(x,x'), we expand in eigenfunctions of the Laplace operator on S 3 , the so-called scalar 
harmonics $AZm with properties (see e.g. [%3~1 1179] ): 



A$ Xlm = (1 - \ 2 )$Xlm 



Our ansatz therefore reads: 



A = 1, 2, . . . ; / = 0, . . . , A - 1; m = -I, . . . , I . (5.65) 



oo A— 1 I 



0&*) = £££* 



Xlm(X CximiX 



A=l (=0 m=-Z 

Introducing this in (|5.58|) we obtain 



oo A— 1 I 

EEE 



A 2 -l 

,4 2 



CAim(^) 



$A*m(z) = ~ ^, ^ - - Uq(x) Uo(x') 



-9 



A=l 1=0 m=-l 

If we now multiply by $\'i> m >(x) y/r\ and integrate over S 3 we find 



A' 2 — 1 1 — — 



I 

S 3 



(5.66) 



(5.67) 



(5.68) 
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where we made use of the orthogonality relation of the scalar harmonics 

v /ry$ Aim (f) ^x'i'm'(x) d 3 x = Sxx>6 U 'S mm t. (5.69) 



In the case A' = 1 (and vanishing I' and m') the above equation is identically satisfied 
for all values of Ci 00 (x'). 9 

In the case A' 7^ 1 the coefficient c\>i> m i(x') follows to be 



^X'l'm'(x') 



A{1 - A' 2 ) 

so that the formal solution for the modified Green's function can be written as 

00 A— 1 / 



(5.70) 



Q(x, x') = E E E • (5.71) 



A=2 1=0 m=-l 



A(l - A 2 ) 



Due to the maximal symmetry of the 3-sphere, the Green's function Q(x,x') can only 
depend on the geodesic distance s(x, x') G [0, 7r] between the two points x and x': 



coss = cos^! cos^ + sinyjx sin^ cos/3 , 
cos/5 = cos</>2 cos</?2 + siny9 2 sin<// 2 cos(<^3 — ip' 3 ) 



(5.72) 



which in the case <p 2 = ¥2 an d ^3 = f'3 clearly reduces to s = (fi — (p' v Indeed, the sum 
over I and m can be performed using 10 



A-l I 



E E ^ m (£)$ 



A sin (As) 
2tt 2 sin s 



1=0 m =-l 

such that even the remaining sum over A can be done analytically: 



Q(s) = Q{x,x* 



1 



1 



1 ) cot s 

7T 



8n 2 A AttA 

To interpret this result, we develop Q{s) around s = obtaining 11 



s E [0, tt] . 



CM = 3 



1 3 S r 

47TS 87T 2 127T 



(5.73) 



(5.74) 



(5.75) 



9 This means that we have the freedom to choose c\qq(x') freely. Our choice is c\qq(x') = without 
restricting generality since from eq. I|5.59|l we immediately conclude that for sources satisfying i|5.57f) 
there will never be any contribution to u(x) coming from cioo(x'). 

10 See e.g. PHI]. 

n Note that the constant term in the expansion (|5.75|) is specific to our choice of c\qo{x'). 
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By introducing the variable r = sA and by treating Q (s) as a gravitational potential we 
find 

We notice that for short distances, r < A, we find the expected flat result whereas the 
corrections to Newton's law become important at distances r of the order of A in the 
form of a constant attracting force. 12 



5.4.3 Green's function of a massless scalar in the background 
(EH} 

We are now prepared to address the main problem of this chapter namely the com- 
putation of the modified Green's function of a massless scalar field in the background 
space-time 1)5.1)) . Since our main interest focuses again on the low energy properties 
of the two-point function, we will limit ourselves to the static case. The main line of 
reasoning is the same as in the previous section. Due to the fact that our space has the 
global topology of a 4-sphere S 4 , we again are confronted with a non-invertible differential 
operator. We would like to solve 



Vu(<fi,9) = 



Au((fi 



+ 



1 



d 



Rll(e) 2 R 2 a(9) 1 (9fd9 



a{9) 1 {9) 



3 du(^j,9) 
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j(<Pi,9), (5.77) 



where the independent angular variable ranges are < ip± < ir; < Lp 2 < < Lp% < 
2tt] — 7r/2 < 9 < 7r/2. Every discussion of Green's functions is based on Green's identity 
relating the differential operator under consideration to its adjoint operator. Since T> is 
formally self-adjoint we have 



TT 2tT 



dipi J dif>2 J df3 J d9y/—g (Vv)u-v(Vu) 

I) 



(5.75 



7T 27T 



dip! J dip 2 I dip 3 — ^ 











+ . . . 



12 There exists numerous articles treating the gravitational potential of a point source in Einstein's 
static universe (see [1801 1181j and references therein) . We only would like to point out here similarities 
between our result (|5.74|l and the line element of a Schwarzschild metric in an Einstein static universe 
background given in |18()| . 
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where we dropped the arguments of u and v for simplicity. The dots in (|5.78|) refer to 
boundary terms in the variables </?j and since we again employ an eigenfunction expansion 
in scalar harmonics on S 3 , these boundary terms will vanish. In order to find an integral 
representation of the solution u(<pi, 6) of ()5.77|) we have to impose appropriate boundary 
conditions on u and v at 9 = ±7r/2. For the time being we assume this to be the case 
such that all boundary terms in (J5.78|) vanish and proceed with the formal solution of 
(J5.77|) . We will address the issue of boundary conditions in 6 in detail in appendices IU1 
and 

In the following we will collectively use x instead of ((pi, 9). The homogeneous equa- 
tion T>u(x) = does not have a unique solution under the assumption of periodic 
boundary conditions. In addition to the trivial solution (u = 0) we also find 

Vuo(x) = 0, J \J—g uq{x) uq{x) d 4 x = 1 . (5.79) 

Therefore, the corresponding inhomogeneous equation (|5.77|) does not have a solution 
unless we again restrict the space of allowed sources: 

f—gj(x) uo(x) d 4 x = . (5.80) 

As in the last section, this condition can be obtained by replacing v by the non-trivial 
solution of the homogeneous equation u in (|5.78|) . We now define the modified Green's 
function by 

V x G(x, x') = 6 - u (x)u (x') . (5.81) 

From Green's identity (|5.78|) . with v (x) given by Q(x, x'), we again obtain after complex 
conjugation the desired integral representation: 



u(x) = Cu (x) + I y/^jj(x')Q(x', x) d A x' , (5.82) 

with 



C= / V~9 u(x') uq(x') d x . (5.83) 



As before, the source 



Jo (x) = [X * } - u (x) uo(x') (5.84) 
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satisfies the solvability condition (|5.8U|) by construction. The normalization of the con- 
stant mode Uq(x) is slightly more involved than before due to the nontrivial measure 
a(8) / -f(8) 3 in the 9 integration. By inserting 



u (x) = iVo$ioo(v\)Xi(0) with xi(0) = 1 
in the integral in ()5.79|) we obtain 



(5.85) 



N 



2lu 



RRfjt&nh(±f) 
For the solution of eq. (|5.81j) we use the ansatz 



(5.86) 



oo A— 1 I 



C\lm{tp'i, 9, 9') 



(5.87) 



A=l 1=0 m=-l 



where from now on we decide to write all arguments explicitly. After inserting this in 
(I5~£T1) we find 



oo A-l I 
A=l 1=0 m=-l L 



+ 



A 2 — 1 



1 1 d 



R 2 a(9)^(9f 89 
5^-^8(9-9') 



dc xlm (ri,6,6') 
89 



$\lm(<Pi) 



-U (<fi,9) Uq^O'). 



(5.88) 



After multiplication by ^^yi'm'ifi) and integration over S 3 we obtain 

^dcx'Vm'(^'i,9,9') 



C\>l'm'{<Pi,9,9 ) + 



1 1 



R 2 a(9) 1 (9fd9 _ 



a(9) 1 (9Y 



89 



(5.89) 



R 3 u Ra(9) 1 (9f 
We now have to distinguish the cases A' = 1 and A' ^ 1. 



0')®\'i'm'{<p'i) ~ u {<p'i,9') / y/rju (<pi,9) $xi, m ,((pi)d 3 <pi. 
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1. (A7W) = (100). 

In this case, the last term on the right hand side of (|5.89|) will give a non- vanishing 



contribution: 
1 



8 



R 2 a{9) 1 {9) 3 89 
1 1 



°{0)i{e) 



oe 

5(9 - 0')$ioo(^) 



=1 



By defining g^'(9, 9') by the relation 



RfjR 



we obtain the following differential equation for g^(9, 9'): 



8 



R 2 a(9)~/(9) 3 89 
where we used 



a{9) 1 {9) 



> dgV(9,9') 
89 



9') 



a{9) 1 (9f 



Xi(0) = 



2oo 



tanh (!f) 

Note that we defined Xi(9) in such a way that 



xi(9), (xi(e) = i) 



(5.90) 



(5.91) 



Xi(9)Xi(9'), (5.92) 



(5.93) 



a{9) 1 {9fxi{9)Xi{9)d9 = l. 



(5.94) 



2. (A'Z'm') ± (100). 

Due to the orthogonality of $i o and $A'J'm' on S 3 , the last term on the right hand 
side of ()5.89|) vanishes. We therefore have 



8 



Rf J Ra(9) 1 (9y 



R 2 a(9)~/(9) 3 89 
5{9-9')$ x , Vm >{<p'i)- 



a{9) 1 {9) 



,dCx'Vm'{^'y9,9') 

89 



(5.95) 
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Introducing g( x '\9,9') again via 



RfjR 



(5.96) 



we see that g( x '\9,9') has to satisfy 

1 1 

& a{9) 1 {9) 3 09 



^ {X '\0,9') + 1-—^ ° 



Rli{0) 



-9 y 



a{9) 1 {9) 



3 dgW(9,9') 



89 



5(9 



a(6)y(6y 



(5.97) 



Combining the above results for A 
solution of (|5.81j) as 



1 and A / 1 we are able to write the formal 



G(<p i ,<p' i ,e,o') 



i 



A-l 



Rf r R 

U A=l 1=0 m 



(5.98 



We obtain a further simplification of this formal solution by employing the spherical 
symmetry on S 3 , see eq. ([5.73)1 . leaving us with a representation of the two-point function 
by a Fourier sum, which is natural for a compact space without boundaries: 



i 



oo 



RRl ^ 2tt 2 



A sin(As)^ (A)( 



,9'\ 



A=l 



sin s 



(5.99) 



Finding the formal solution ()5.99|) was straightforward apart from minor complica- 
tions inherent to the use of a modified Green's function. To solve the differential equa- 
tions (|5.92|) and ()5.97|) by imposing appropriate boundary conditions again is a routine 
task without any conceptual difficulties, though slightly technical in nature. Replacing 
back this solution in (|5.99|) we are left with a Fourier sum which at first sight looks 
intractable, anticipating the fact that the solutions g( x >(9,9') (for A ^ 1) are given by 
hypergeometric functions. Nevertheless it is possible to extract the desired asymptotic 
information from the sum ()5.99|) . 

In order not to disturb the transparency and fluidity of the main discussion, we 
provide large parts of the technical calculations in four appendices. In appendix O we 
report similarities and differences in the evaluation of the two-point functions between 
our case and the case of Randall- Sundrum, since the latter served as a guideline for 
handling the more difficult case under consideration. The solutions of eqs. ()5.92|) and 
(15.97)1 are presented in appendix \D\ and the evaluation of the Fourier sum ()5.99j) at 
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distances exceeding the size of the extra dimension in appendix |E| Eventually, appendix 
|F] contains the evaluation of the Fourier sum (|5.99|) for distances smaller than the extra 
dimension. In this way, we can offer the reader less interested in the details of the 
computations to have the main results at hand. 

From its definition (J5.81|) we understand that the Green's function (J5.99|) can be 
considered as the response of the scalar field to the combination of a point-like source 
located at coordinates (cp^, 9') and a delocalized, compensating negative contribution. 
Since we would like to see the response to a point-like particle on the brane, we put 
9' = in eq. ()5.99|) and explicitly write the A = 1 term: 

^• 0) = CS+4^^ SM ^ (5 ' 100) 

with S[s, 9, uj] given by (|K.1|) of appendix|E|(see also (|D.13J) of appendix^- The general 
result for the sum S[s, 0,u] obtained in appendix lEl is 

S[s, 0, oj] = lim S[s, 9, lo] = ( ~ cos s si* 1 s J (5.101) 



o — >o ujI- z(0) V 2 4 

- — *(0) - * l n [i _ 2(0)] sin s - -*(0)5 lim R[s, 9, u] , 

where we have z(0) = tanh 2 (^) and where we refer to (jE.9|) for the definition of 
R[s,9,u]. The first term in (J5.101|) is the zero mode contribution 13 and we see that it 
reproduces exactly the 4-dimensional static Green's function of Einstein's static universe 
given in ([5.74)1 . The other two terms are the contributions from the higher Kaluza-Klein 
modes. 

We first concentrate on the case where s ~ 1 or what is equivalent t ~ Rjj . Since one 
can easily convince oneself that i?[s,Q,o;] ~ [1 - z{0)Y in this regime, the contributions 
of the higher Kaluza-Klein modes are strongly suppressed with respect to the zero mode 
contribution. This means that as in the case of the Randall-Sundrum II model it is the 
zero mode which dominates the behavior of gravity at distances much larger than the 



13 At this point we have to explain what we mean by "zero mode" in this context, since due to the 
asymmetrically warped geometry, the spectrum of Kaluza-Klein excitations will not be Lorentz-invariant 
and strictly speaking, different excitations cannot be characterized by different 4-dimensional masses. 
It would be more accurate to say that for a given value of the momentum eigenvalue A, there exists a 
tower of corresponding Kaluza-Klein excitations with energies given by with n = 0, 1, . . .. In our 
use of language the "zero mode branch" of the spectrum or simply the "zero mode" is defined to be the 
collection of the lowest energy excitations corresponding to all possible values of A, that is by the set 
{ (Eq , A) , A = 1, 2, . . .}. The definition is readily generalized to higher branches of the spectrum. 
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extra dimensions R. The main difference to the Randall-Sundrum II case is that the zero 
mode of our model not only gives rise to the typical 4-dimensional 1/r singularity but 
also accounts for the compactness of space by reproducing the Einstein static universe 
behavior ([5.74)1 . This result is somewhat surprising given the extreme anisotropy of our 
manifold. Due to the fact that the distances between two arbitrary points on the brane 
are of the order of R, one might intuitively expect that the extra dimension can be 
effective in determining also the large distance behavior of gravity (on the brane). As 
we could show by direct calculation the above expectation turns out to be incorrect. 

Next, we consider physical distances r much larger than the extra dimension r ^> R 
and much smaller than the observable universe r -C Ru in which case the results for 
R[s,0,u] can be seen to be: 



7i 7i 1 — z(0) r 
R[s,0,uj}~— + - g4 l ; {8-6hi2-61n 



8 [1 - Z(0)] 



-1/2 



(5.102) 



+ o 



[1 - z(0)f 



In 



s[l-z 



-1/2 



valid for [1 — ^(O)] 1 ^ 2 C s C 1. The zero mode contribution (to S[s, 0, lu]) in this regime 
is simply the constant obtained by setting s = in the first term of 1)5.101)1 so that we 
obtain: 



S[s, 0, lu) 



Ti z(0)2 f 1 1 r , , m ( ns 

K - 1 T -r ll + ^r-, r + -r4-31n2-31ns+0 ^- 

lu 1 - z(0) 1 2s 2 s 4 L J V s e 



(5.103) 



where we introduced s 



s\l-z 



physical distance r = R v s instead of s we obtain 

£M,o) = 
s (1) (o,o) 



1 ^ 2 . Inserting this result in ()5. 100)1 and using 

(5.104) 



1 wcoth(^) f 1 1 r , , . ^ / 

, ^ 1 + ^ + -4-3ii2-3es+G 

27i 2 RRl 4vrr R \ 2s 2 s 4 L V 

g m (0,0) 1 ^coth(^f) f tanh 2 (if) tanh 4 (if) ^ 
27T 2 RRfj ~ R I + 2f2 + V 4 X 



In s 



X 



4-31n2-31n 



tanh (f) 



+ o 



tanh 6 (<f) 



In 



tanh (if) 



Since we would like to compare our result with the corresponding correction in the 
Randall-Sundrum II case, we introduced the dimensionless distance variable f = ru/R, 



101 



the physical distance measured in units of the AdS-radius, in the last line of the above 
result. We see that in complete agreement with the Randall-Sundrum II scenario, our 
setup reproduces 4-dimensional gravity at large distances with extremely suppressed 
corrections. The only remnant effect from the different global topology manifest itself 
through the factors of tanh {jrf) and coth f^J which are very close to 1. We furthermore 
emphasize that apart from these deviations, the asymptotic we obtained coincides exactly 
with the asymptotic for the case of a massless scalar field in the Randall-Sundrum II 
background, see ijCTijl and e.g. jTTSl HHl HEI EES- From the factor cucoth (*f) /R 
in eq. (|5.104j) we see that also the relation between the fundamental scale M and the 
Planck-scale Mpi gets modified only by the same factor of tanh {jrf)'. 

M^ = M 3 ^tanh(^) , (5.105) 

where we remind that R/u is nothing but the AdS-Radius. 

Eventually, we treat the case of distances inferior to the extra dimension r <C R- 
After using the result (JF.12j) of appendix |F] in (j5.100jl . a short calculation reveals 

Q ^ °' 0) = 2^RRj - + W (5 ' 106) 

a result that has to be compared to the characteristic solution of the Poisson equation 
in 4-dimensional flat space. In n-dimensional flat space one has: 



A 



(n - 2)V S n-ir 



n-2 




6(r), r= > x\\ , (5.107) 



with Vgn-i = 27rt/r[|] denoting the volume of the n — 1 sphere. Specifying to n = 4, 
we recover the correct prefactor of — 1/47T 2 in (|5.106|) multiplying the 1/r 2 singularity. 

Finally, we mention that in none of the considered cases we payed attention to the 
additive constant in the two-point function on the brane. The arbitrary constant entering 
the solution g^{9, 9') can always be chosen in such a way that <7^(0, 0) vanishes on the 
brane (see appendix ID| . 

5.5 Conclusions 

In this chapter we considered a particular brane world model in 5 dimensions with the 
characteristic property that the spatial part of the space-time manifold (including the 
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extra dimension) is compact and has the topology of a 4-sphere S A . Similar to the original 
Randall-Sundrum II model, the 3-brane is located at the boundary between two regions 
of AdSs space-time. The coordinates of AdS$ used by Randall and Sundrum are closely 
related to the so-called Poincare coordinates of AdS$. While the extra dimension in this 
set of coordinates provides a slicing of AdS 5 along flat 4-dimensional Minkowski sections 
(resulting in a flat Minkowskian induced metric on the brane), their disadvantage is that 
they do not cover the whole of AdS§ space-time. 

The coordinates we used are the global coordinates of AdS§ known to provide a global 
cover of the AdS§ space-time. In this case the "extra" dimension labels different sections 
with intrinsic geometry Rx 5 3 , the geometry of Einsteins static universe. The induced 
metric on the 3-brane in our setup is therefore also given by K. x S* 3 . 

As we illustrated with the use of the Penrose-diagram of AdS$, the incompleteness 
of the Poincare patch is at the origin of the incompleteness of the Randall-Sundrum II 
space-time with respect to time-like and light-like geodesies. Moreover, we were able 
to demonstrate that the setup considered in this chapter provides an alternative to the 
Randall-Sundrum II model which does not suffer from the drawback of being geodesically 
incomplete. The latter point was part of the main motivations for the study of this setup. 

The spatial part of our manifold is characterized by an extreme anisotropy with 
respect to one of the coordinates (the extra coordinate) accounting for thirty orders of 
magnitude between the size of the observable universe and present upper bounds for the 
size of extra dimensions. 

Another interesting property of our manifold, related to the anisotropy of its spatial 
part, is the fact that any two points on the brane are separated by a distance of the order 
of the size of the extra dimension R regardless of their distance measured by means of the 
induced metric on the brane. Despite the difference in the global topology, the properties 
of gravity localization turned out to be very similar to the Randall-Sundrum II model, 
though much more difficult to work out technically. We computed the static (modified) 
Green's function of a massless scalar field in our background and could show that in 
the intermediate distance regime R <C r <C Ru the 4-dimensional Newton's law is valid 
for two test particles on the brane, with asymptotic corrections terms identical to the 
Randall-Sundrum II case up to tiny factors of tanh(^). We could also recover the 
characteristic 5-dimensional behavior of the Green's function for distances smaller than 
the extra dimension r ^ R. Eventually we saw that in the regime of cosmic distances 
r ~ Ru, somewhat counterintuitive given our highly anisotropic manifold, the Green's 
function is dominated by the behavior of the corresponding static (modified) Green's 
function in Einstein's static universe. 

In the simple setup considered in this chapter, the 3-brane is supposed to be mo- 
tionless. In the light of recent progress in the study of 4-dimensional cosmic evolution 
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induced by the motion of the brane in the bulk, it would be interesting to explore this 
possibility and see what kind of modifications of our results we would have to envisage. 
Finally, a related, important question which would be interesting to address would be 
the question of stability of our setup. 
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Concluding remarks 



Topological defects or topological solitons have important applications not only in 4- 
dimensional condensed matter physics, particle physics and cosmology, but also more 
recently in theories with extra dimensions, as sketched in this work. There are various 
reasons to study topological defects in the latter context: to start with, we mention 
their classical stability. Strictly speaking, the topological arguments in favor of classical 
stability are valid only in the case of a flat space-time. However, in certain cases of 
topological defects weakly coupled to gravity, we can expect to obtain good candidates 
for stable configurations. Another reason why topological defects are interesting for 
model building in theories involving higher dimensions is that they explain the chiral 
nature of the observed 4-dimensional fermions. The higher dimensional Dirac operator 
in a topologically non-trivial background contains a certain number of zero modes which 
form a 4-dimensional point of view can be identified as Standard Model fermions. It 
would be desirable to reproduce in this way the field content of the Standard Model 
from a higher dimensional theory. 

The introduction of extra dimensions in any theory always has to be accompanied 
by an explanation for why our world appears to be 4-dimensional at low energies. One 
crucial test of such theories when gravity is concerned is whether they are able to predict 
the correct 4-dimensional Newton's law in the distance regime explored by latest obser- 
vations. At this point, non-factorizable metrics play an important role since it has been 
show recently that backgrounds involving such metrics do pass the test of reproducing the 
correct 4-dimensional behavior of the gravitational potential. It is the interplay between 
these two ingredients, topological defects and non-factorizable geometries in higher di- 
mensions, which lies at the center of interest of this work. After an introductory chapter 
on non-factorizable metrics, we presented some possibilities of topological defect solu- 
tions to the higher dimensional Einstein equations. The results presented in chapter |3] 
were mainly of an asymptotic nature in the region of transverse space far from the core of 
the defects. In chapter 0] we gave a numerical realization of one of the metric asymptotics 
found in chapter |3] in the case of a magnetic monopole residing in transverse space, as a 
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particular example of a bulk p-form field. The main motivation for the study presented in 
chapter was the incompleteness of the Randall- Sundrum II space-time with respect to 
time-like and light-like geodesies. As we showed, the brane-world model we discussed is 
geodesically complete and still reproduces the correct 4-dimensional behavior of gravity 
on the brane. 

Needless to say that numerous open questions remain. As mentioned above, the sta- 
bility of all solutions of the Einstein equations given in this work has to be analyzed 
carefully. A related problem would be the study of the dynamics, in a cosmological 
context for example, since all cases considered correspond to static solutions. The ver- 
ification of Newton's law on the brane is only one necessary consistency check. More 
generally it should be verified that 4-dimensional general relativity can be recovered on 
the brane. 

To summarize, higher dimensional brane-world scenarios provide numerous oppor- 
tunities for addressing many unresolved problems of theoretical particle physics and 
cosmology. However, the introduction of extra dimensions always has to be accompa- 
nied by an explanation for why they are not observed at accessible energies. Despite 
numerous attempts in this direction, no theoretically and phenomenologically satisfying 
model exists at present. It is interesting to see in which way new collider experiments 
will change our present view of extra dimensions. 
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Appendix A 

A generalization of the metric ( 12.11 ) 



In this appendix we plan to consider a particular metric ansatz which can be considered 
as a generalization of the metric (|2.1|) given in chapter |21 As we will see, the Einstein 
equations with a cosmological constant term in otherwise empty space-time are very 
restrictive in the case we are going to study. The metric generalization which we will 
briefly describe is used in theories of higher dimensions characterized by two scale-factors 
(in the cosmological sense). In this context of dynamical compactification |18Hj - [T9T] . the 
cosmological expansion of our 4-dimensional universe is thought to be accompanied by 
a shrinking of the size of an internal manifold. 

The form (|2.1j) of the metric can still be generalized in the following sense 

ds 2 = g MN dx M dx N = a (x a ) g^ (x p ) dx» dx u + rj (a") g ab (x c ) dx a dx b , (A.l) 

introducing a second warp factor 77 (x p ) depending only on the coordinates x^. In this 
case the connections become: 

f c ab = t c ab , f° 7 = ^-V 7 r^, ?r b = -±V'Tig ttb , (A.2) 
Similarly, the Ricci tensors Rmn, R^v, Rab and Ricci scalars R, R and R corresponding 
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to Qmn i 9»v and g ab generalize to: 



R 



2a 



do 



~ — V M V„ T) + — V M T]V U T]. 

2rj 4r]^ 



Rab — Rab ~ ~ — 9ab 



d 2 - 2. 
2r] 



2 a 4 a 2 



(A.3) 



(A.4) 



R 



d x + d 2 -2 
4arj 



K=-K Va V cr — V a (TV a 

rj ar) Aa z r] 

1 „ rfo _ _„ do (do — 3)„ _„ 

cr 77 a Ar] z a 



(A.5) 
(A.6) 



In the case 77 = 1 the relations (|A.2J) - (|A.6J) obviously reduce to (J2.2j) - (J2.6j) . 

For the generalized metric (|A.1|) the Einstein equations ()2.7|) turn out to be more 
restrictive, at least in the case of vanishing energy-momentum tensor, see for related 
work in 4 dimensions. Apart from the physically uninteresting case d\ + d 2 = 2, the off- 
diagonal components ^ a of (|2.7jl imply that either a or 77 has to be constant, as a direct 
consequence of (IA.5J) . To some extend, this brings us back to the original metric (|2.1jl . In 
particular with the choice 77 = 1, (jA.l|) reduces to ([2.1)1 . The case a = 1 has been studied 
in the context of dynamical compactification of extra dimensions |183j - [T!?T] . Due to the 
formal invariance of the metric under the following interchange of "conventional" and 
"internal" dimensions 



sr 



Qab [X 



(Tlx 



7/(x A 



(A.7) 



the Einstein equations for this case can be obtained from the Einstein equations 
(12.10)1 by generalized replacements of this kind. However, it is more useful to present the 
equations ()2.7j) in a somewhat more suggestive way: 



Rfj,u ~ ~ Qfiu R + A g^ u 



di d 2 _ — 

— V M V„ 77 - — V M 7/ V„ 77 
2 7/ 4 77^ 



R 



ab 



g ab R + A g ab = , 



(A.8) 
(A.9) 
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where 



A = d2 ( t 2 = 3) V^V^ + ^V^ + Ab-^, (A. 10) 

} Tj 2 7] 2 7] 



A = A B n + {d2 ~ ^ ~ 4) v», v>, + v», - 1 a. (a.h) 

Due to the choice a = 1 the off-diagonal components Ma are identically satisfied. We 
are now interested in solutions to the above equations (|A.8|) and (jA.9|) in which the 
dependencies on the coordinates are separated. In other words, we would like to eliminate 
the dependence of A in (jA.10|) on the "internal" coordinates x a via R. This is equivalent 
to imposing that the internal space should be a constant curvature space R ^ R(x a ). 
By taking the trace of eq. (jA.9|) . this by itself implies that also A has to be a constant. 
It is now an algebraic exercise to take the trace also of eq. (|A.8jl in order to eliminate 
the dependencies of A and A on R and R, respectively. The result is 

A = . d !~ 2 9 Ab + £ V A - A Va ^ r, , (A.12) 
d\ + d 2 — 2 4 T) 8r] z 

A = /:; 2 - A B V + V A „V* v + V A VS . (A.13) 

«i + «2 — 2 8 77 4 

The inverse powers of 77 in (|A.12|) suggest the following redefinition of the warp-factor: 

■q(x?) = e 2p(xM) , (A.14) 
so that the "effective" 4-dimensional Einstein equations (|A.8|) take the form: 



^ ~ 2 g ^ uR+ d l + d 2 -2 Ab9 ^ v = T ^ u ' (A.15) 
where we introduced by 

T^^^jtV^pV.p + V^V.pl-i^tV^pV^ + V^V^j . (A. 16) 

Eq. ()A.15|) has to be supplemented by the highly non-linear "field" -equation (jA.13|) for 

the warp factor p: 

V A VV + *V,pV^- s H_ e -*A + ;fT i- 5 A B = 0, (A.17) 
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where A is related to the Ricci scalar R of the internal constant curvature space by 

A-^A. (A,8) 

From eq. (|A.15J) we see that the effective 4- dimensional cosmological constant is a fraction 
of the bulk cosmological constant and that from a 4-dimensional point of view the warp- 
factor plays the role of a scalar field coupled to Einstein-equations via the effective energy- 
momentum tensor (|A.16[) . Note that in addition to the standard energy-momentum 
tensor of a massless scalar field given by the two first terms in the square brackets of 
eq. (jA.16j) there are also two non-standard contributions given by the two last terms in 
the square brackets. Since we consider a pure gravity theory it is needless to say that 
there is no Planck-scale involved in the coupling of gravity to the scalar field p. The 
constant A is a free parameter of the theory. 
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Appendix B 

Fine-tuning relations for the 
't Hooft-Polyakov monopole 



By taking linear combinations of Einstein equations ()4.29|) to (|4.31|) one can easily derive 
the following relations: 

[M(r) 3 M'(r)£(r) 2 }' « 

M{rYC{ry = -5(27 + eo-e p -2e e ), (B.l) 

{M(r)*£(r)£'(r)]' 1 (3 

WWW) " W? " -5 (27 - 4e °-^ + 3ee) - (R2) 

By multiplying with M(r) 4 £(r) 2 , integrating from to oo and using the definition of the 
brane tension components ()4.44|) we obtain 



oc 



M(r) 3 M'(r)£(r) 2 \™ = J M(r) A £(r) 2 dr (B.3) 



M(r) 4 £(r)£'(r)\™ = jM(r)*dr -^2 jM(r) 4 £(r) 2 dr (B.4) 

o o 

~ (4fi + fi p - 3 fi e ) . 

Using the boundary conditions for the metric functions (|4.45jl and ()4.59|) . and taking 
the difference of the eqs. flB.3|) and (jB.4|) then establishes the first part of eq. (|4.72j) . To 
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prove the second part of ([4.72)1 . we start from the general expressions of the stress-energy 
components e i; relations ()4.37j) - (j4.39j) : 



jlQ - fig 



drM{rfC{r) 2 {e e - e ) 



(B.5) 



drM{r) 4 C{rf 



K'{rf [1-K{r) 2 ] 2 J{r) 2 K{rf 



C(r) 



C(r) 



Cr) 



Multiplying the equation of motion for the gauge field (|4.33|) by K(r) and substituting 
the J(r) 2 K(r) 2 term gives 



Ho ~ He 



drM{rf { K\r 



-I i„\2 



l-K(r) 2 [M(r) 4 K'(r)]' K(r) 



C(r) 2 M(r) 4 
o 

Integration by parts in the last term of the above equation leads to 

oo 

"1 -K(r) 2 



vB.6) 



Ho ~ He 



drM(r) 



C(r) 



+ [K(r)M(r) 4 K'{ 



(B.7) 



which together with the behavior of the gauge field at the origin K'(0) = (see (|4.52Jl ) 
finishes the proof of relation (|4.72|) . 

The proof of relation ()4.73|) is simply obtained by rewriting f)B.3|) . 

To establish (|4.74|) we start directly from the definitions of the stress-energy tensor 
components, relations ()4.37|) - (j4.39j) : 



Ho + H P + 2fi e = J drM(r) 4 C(rf 



2 J(r) 2 K(r) 
£(r) 2 



+ a(j(r) 2 -l)' 



(B.8) 



Collecting derivatives in the equation of motion for the scalar field (|4.32|) and multiplying 

by J{r) gives 

2J(r) 2 K(r) 2 [M (r) 4 C(r) 2 J 1 '(r)]' 



-J(r)-aJ(r) 2 (j(r) 2 - l) . 



£(r) 2 M(r) 4 £(r) 2 
Eliminating now the second term in the equation (|B.8|) leads to 

[M(r) A C(r) 2 J'(r)] f 



(B.9) 



Ho + H P + = J drM{r fC{rf 



J'(r) 2 + 



M(r) 4 £(r) 



J(r) + a (1 - J(r) 2 ) 

(B.10) 



114 



If we now expand and integrate the second term by parts we are left with 

oo 

H + H P + 2 f i 9 =a j dr (l - J(rf) M(rfC(rf + [M (rf £(rf J (r) J' (r)]™ , (B.ll) 
o 

which reduces to (|4.74|) when the boundary conditions for the metric (|4.45|) and ()4.59j) 
are used. 
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Appendix C 



Computational parallels in the 
corrections to Newton's law 



The purpose of this appendix is to review briefly the calculations of the static two-point 
function of a scalar field in the Randall-Sundrum II background |113H114] and to compare 
each stage with the corresponding stage of calculations in the background considered in 
chapter This serves mainly for underlining similarities and differences between the 
two calculations. Let us begin by writing down the metric of the Randall-Sundrum II 
model 

ds 2 = e-^ri^dx^dx" + dy 2 , -vrr c < y < vrr c . (C. 1) 

Here, y stands for the extra dimension while k and rj^ denote the inverse radius of AdS^ 
space and the (4-dimensional) Minkowski metric with signature — h ++. Due to the 
orbifold symmetry, the allowed range of y is < y < isr c } As it is well known, each 4- 
dimensional graviton mode in this background satisfies the equation of a massless scalar 
field. Therefore, for the study of the potential between two test masses on the brane 
we confine ourselves to solving the equation of a massless scalar field with an arbitrary 
time-independent source: 

Vu{x, y) = j(x, y) with V = e 2 ^A x . - 4fc|- + |^ . (C.2) 

Since the operator T> is formally self-adjoint, the corresponding Green's function will 
satisfy: 

Vg(x,x;y } y) = -j= . (C.3) 

^^We allow for a finite r c only to impose boundary conditions in a proper way. Eventually we are 
interested in the limit r r — > oo. 
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We are now able to write down the usual integral representation of the solution of (|C.2|) : 



u(x, V) = J V^9J(x, y') G(x, y, y') (fx'dy', (C.4) 

where the ^/-integration extends from to irr c and the x' integrations from — oo to 
oo. The absence of boundary terms in (jC.4|) is of course the result of an appropriate 
choice of boundary conditions for u(x, y) and G(x, x'; y, y'). In the above coordinates of 
the Randall- Sundrum II case the orbifold boundary conditions together with the Israel 
condition imposed on the fluctuations of the metric give rise to a Neumann boundary 
condition at y — 0. One can easily convince oneself that in the limit r c —>■ oo the resulting 
Green's function is independent of the choice of the (homogeneous) boundary condition 
at y = irr c . We therefore follow |H3j and use also a Neumann boundary condition at 
y = irr c . Another important point is that the Green's function we are considering is 
specific to the orbifold boundary condition and so describes the situation of a semi- 
infinite extra dimension. We decided to carry out the calculations in the semi-infinite 
case as opposed to [10], where eventually the orbifold boundary conditions are dropped 
and the case of a fully infinite extra dimension is considered. 2 An immediate consequence 
of this will be that the constant factors of the characteristic short distance singularities of 
the solutions of Laplace equations also will be modified by a factor of 2. Finally, we note 
that this factor of 2 can be accounted for by an overall redefinition of the 5-dimensional 
Newton's constant leaving the two theories with equivalent predictions. We conclude 
the discussion of boundary conditions by noting that at infinity in x' we suppose that 
u(x',y') and G(x,x';y,y') vanish sufficiently rapidly. 

The solution of eq. (|C.3J) can be found most conveniently by Fourier expansion in the 
x coordinates and by direct solution of the resulting differential equation for the "trans- 
verse" Green's function. After performing the integrations over the angular coordinates 
in Fourier space we obtain: 

oo 

G(x-x',y,y') = -L f Pg^igzjl) W ( , )d (c . 5) 

Z7T J \X — X'\ 



where g^ p \y,y') satisfies 

_ p V V P W) = , (0<y<vrr c ) (C.6) 

2 However, as long as matter on the brane is considered, the Green's function for the two setups differ 
only by a factor of 2. 



d 



e -iky Q y 



-Aky 



d_ 

dy 



g {p) (y,y') 
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together with Neumann boundary conditions at y = and y = irr c . The solution of 
(|C.6|) is straightforward with the general result 



s (p W) 



a 2k(y > +y < ) 



x 



(C.7) 



x 



X 



p 



P 



fc 



Ji ( ^e fcwrc 



l, 2 (| e ^ 



x 
A; 



where ?/> (?/<) denote the greater (smaller) of the two numbers ?/ and y' and Ii(^), ^2(2), 
K 2 (z) are modified Bessel functions. Since we want to study sources on the brane, 
we now set y' = and take the well-defined limit r c — > 00, as can easily be verified from 
the asymptotic behavior of I n and K n for large values of z. We are therefore able to 
write the static scalar two-point function as 



G(x-x;y, 0) 



00 

e 2fc f /s in (plx-^l) AT 2 (fe fc?/ ) 

Mir 



2vr 2 



(C.8) 



We would like to point out a particularity of the Fourier integral (jC.8|) . By using the 
asymptotic expansions for the modified Green's functions Ki(z) and K 2 (z) we have 



K2 (fg) 
Ml) 



~ e 2 



l^p-f(^-i) 



O 



(C.9) 



showing that for y = the integral over p in (|C.8Jl does not exist. The correct value of 
the Green's function on the brane is therefore obtained by imposing continuity at y = 0: 



Q(x - f ; 0, 0) = lim Q(x - x; y, 0). (CIO) 

2/->0 

Note that this subtlety is still present in our case of the Fourier sum representation 
of the modified Green's function f)5.99|) . as discussed in appendix lEl 

The representation (jC.8|) is very suitable for obtaining the short distance behavior of 
the Green's function. By inserting the expansion (jC.9|) into (jC.8|) we can evaluate the 
integral which will give reasonable results for distances much smaller than l/k: 

«*-** >~-£jl^+7- (cm) 
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We see that after replacing the exponential factor by 1 (valid for y <C 1/&) we recover 
(up to a factor of 2) the correct 5-dimensional behavior of the Green's function in flat 
4-dimensional space ()5.107|) . 3 

For large distances the Fourier representation (jC.8|) is less suited for obtaining correc- 
tions to Newton's law. This is due to the fact that all but the first term in the expansion 
of K%jK\ in powers of p around p = lead to divergent contributions upon inserting in 
()C.8|) . We therefore seek another method which will allow us to obtain the corrections 
to Newton's law by term-wise integration. 

The idea is to promote (|C.8|) to a contour integral in the complex p-plane and to 
shift the contour in such a way that the trigonometric function is transformed into an 
exponential function. First we introduce dimensionless quantities by rescaling with k 
according to X = \x — x'\k, Y = yk, z = p\x — x'\: 

U2 2Y °r XT (-P Y ) 

3 ^-I^S^^tk) iz - (C12) 



=I[X,Y] 



Following |114j (see also |192p we now use the relation 

K 2 [w]=K [w] + ^K 1 [w] (C.13) 

to separate the zero mode contribution to the Green's function (Newton's law) from the 
contributions coming from the higher Kaluza-Klein particles (corrections to Newton's 
law). Using (ITIT31) in I[X,Y] we find 



I[X,Y) 



. K (f e Y ) j 2X 

sin z ; — T^-az H — tt- 

K X (f ) 



sin zKi (f e y ) 
~ Ki (f ) 



dz . 



(C.14) 



=h[X,Y] 



=HX,Y] 



The additional factor of 1/z in the integrand of ^[X, Y] allows us to take the limit Y — » 
with the result: 

UmI 2 [X,Y]=TrX. (C.15) 

3 As discussed above, the factor of 2 is a direct consequence of the boundary conditions corresponding 
to a semi-infinite extra dimension. 
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Ci 



Figure C.l: Contour used in Cauchy's theorem to evaluate the integral (jC.16|) . 



We still need Y > for convergence in I\[X, Y\. The next step is to use sin z = Q {e lz } 
and to exchange the operation S with the integration over z: 4 



h[X,Y] 



3 



K Ue Y ) 
e lz V f J dz 

Mi) 



(C.16) 



The integrand in (|C.16|) is a holomorphic function of z in the first quadrant (see [192 , 
p. 377 for details). We can therefore apply Cauchy's theorem to the contour de- 
picted in Fig. IC.ll As we will show shortly, the contribution from the arc C*2 = 
{z| z = Re lip , < <f < |} vanishes in the limit R — » oo (the limit we are interested 
in). Symbolically we therefore have 



lim 

R^oo 



Ci 



0, 



(C.17) 



where . . . replace the integrand in (|C.16J) and the sense of integration is as indicated in 
the figure. This means that we can replace the integration over the positive real axis in 
(jC.16|) by an integration over the positive imaginary axis without changing the value of 



4 This is justified since both the real and the imaginary part of the resulting integral converge. 
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I\[X, Y\. Substituting now z in favor of n according to n = 



—iz we obtain: 



h[X,Y] = $i 



K (f e Y ) 
e- n ° V f M dn 
Kr(f) 




K ] (f e Y ) 
HP (f ) 



^(f)] z +^(f)]' 



(C.18) 



In the first line in (jC18|) we replaced modified Bessel functions with imaginary arguments 
by Hankel functions of real argument. In the second line we explicitly took the imaginary 
part after replacing the identities relating the Hankel functions and the Bessel functions 
J and Y of the first kind. The benefit from the rotation of the contour of integration is 
obvious at this stage. First, the last integral in (jC.18|) converges also for Y = due to 
the presence of the exponential function in the integrand. In this case we have 



h[X,0] 



oo 

2X fe- 



ll 



dn . 



(C.19) 



Note that we used the Wronskian relation for Bessel functions to simplify the numerator, 
the general reference being once more |192j . Second, it is straightforward to obtain the 
large X asymptotic (1 X) of the integral (jC.19j) by a simple power series expansion 
in n around n = of the integrand (apart from the exponential factor), followed by 
term- wise integration. For the results see |115j . where the above integral (|C.19|) has 
been found in the wave-function approach: 



h[X, 0] ~ ^ + ^ (4 - 31n2 - 31nX) + O 



lnX 
~X*~ 



(C.20) 



The result for Q(X, 0, 0) is therefore: 



g(x,o,o) 



2vrX 



1 + 



2X 2 X 4 



' (4-31n2-31nX) + /,U A 



X 6 



(C.21) 



However, we still have to verify that the contribution from the arc Ci in Fig. IC.ll 
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vanishes in the limit R — > oo. Parametrizing 2 by z = Re lifi , < tp < tc/2 we find: 



^2 #1 (f J 



< R 



M^) 



iR(cos ip+i simp) I 



k (B^ e y) 



(C.22) 



Making now use of the asymptotic properties of the modified Bessel functions to estimate 
the second modulus in the last integrand in (|C.22|) we can write: 



< c 



e 2 e 



(C.23) 



so that 



c 2 



Mi) 



dz 



< CRe 



y f K[ S in^+(^l)co S ^] 



dip. 



(C.24) 



with C being a constant of the order of unity, independent of R and (p. To finish the 
estimate we observe that 



sin ip + 
and obtain: 



X 



cos (p > e = min 



e Y - 1 
X 



71", 



>0,V^G[0,-], (F>0) (C.25) 



c 2 



Mf) 



< 



-e 2 e 



-fie 



for i? — > 00. 



(C.26) 



This establishes the result that in the limit R — > 00 the arc C2 does not contribute 
to the integral in (|C.16J) and completes our discussion of the two-point function in the 
Randall- Sundrum II setup. 

We are now going to summarize briefly the complications arising when the general 
scheme of computations outlined above for the Randall-Sundrum II case are applied to 
the brane setup considered in chapter El First and foremost, the main difference is due 
to the fact that the induced metric on our brane is not the flat Minkowski metric but 
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the metric of the Einstein static universe. This implies the use of a modified Green's 
function. In the eigenfunction expansion, the discrete scalar-harmonics on S 3 will replace 
the continuous plane-wave eigenf unctions. We therefore expect the formal solution for 
the two-point function (the analog of ([(15)) ) to take the form of a Fourier sum. 

The solution to the transverse Green's function (the analog of eq. ()O.6|0 as presented 
in appendix \D\ turns out to be again straightforward and completes the formal solution. 

When trying to distinguish between the zero mode contribution and the contribution 
coming from higher modes (in the language of the Kaluza-Klein approach) a relation 
similar to (jG.13|) proves useful. While the evaluation of the zero mode contribution poses 
no problems, the corrections coming from higher modes are much more involved this 
time. We face the following major technical difficulties: there is no analog of Cauchy's 
theorem in the discrete case of the Fourier sum. One solution to this problem is to use 
a variant of Euler-Maclaurin's sum rule in order to replace the Fourier sum by an 
analog Fourier integral and a remainder term also in the form of an integral. Now the 
procedure again is similar to the one outlined for the Randall- Sundrum case. Another 
technical problem, however, in connection with large distances corrections has its origin 
in the simple fact that by large distances in our setup we mean distances large with 
respect to the extra dimension R but also small with respect to the size of the observable 
universe Ru- The detailed computations of the corrections in this distances regime can 
be found in appendix lEl 

What concerns the computation of the two-point function at distances smaller than 
the extra dimension, there is no conceptual difference to the Randall-Sundrum case. 
The Fourier sum can be calculated analytically after inserting the corresponding large 
momentum asymptotic. We give the details in appendix IFl 
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Appendix D 



Solutions to the differential 
equations (E92D and (I5TTTT1 ) 



The aim of this appendix is to obtain the solutions to the eqs. ()5.92|) and (|5.97|) and 
hence to complete the formal computation of the modified Green's function. We will 
concentrate mainly on ()5.97|) for two reasons: first and foremost, because in the formal 
expansion (J5.99)) g ( - 1 \9,9') multiplies a constant function in s and since we are mainly 
interested in the behavior of the Green's function on the brane, this constant is of no 
relevance for our considerations. Secondly, because the construction of the Green's func- 
tions g^(6, 9') and g^'(9, 9') follows the standard procedure, so it is not necessary to go 
into details twice. We will merely state the result for g^'(9, 9'). The equation we would 
like to solve is 



X- 



h 0') + 



1 



1 







R 2 a(9) 1 (9fd9 



a{9) 1 {9) 



3 dg {x) (9,9') 



09 



5(9 - 9') 
a{9) 1 {9f 



(d.i; 



where A = 2, 3, 4 . . . and the independent variables 9 and 9' are restricted to the intervals 
— 7r/2 < 9 < 7r/2, and — n/2 < 9' < ir/2. The homogeneous equation associated with 
flD.l|) can be reduced to the hypergeometric equation. Its general formal solution in the 
interval 0<#<7r/2is given by 



<^+)(A, 9) = aS + Vi(A, 9) + 4 + V 2 (A, 9), < 9 < - (D.2) 
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where 



, tanl/" 1 [w (f - 0)1 

y V ; cosh 4 [u (f-0) 



A + 3 A + 3 x 1 u2 
— — , — — ; A + l;tanh 

A+3 A+3 



3; cosh 



TT 



(D.3) 



Due to the symmetry of the background under 9 — > —9 we can immediately write down 
the general solution of the homogeneous equation obtained from (jD.lj) in the interval 
-tt/2 < 9 < 0: 



p<->(A, 0) = 4"Vi(A, -0) + a^V 2 (A, -0) 



TT 



< < 

2 - ~ 



(D.4) 



We now have to decide whether we want to restrict our Green's function to describe 
perturbations which also possess the symmetry 9 — > —0 of the background (as in the 
Randall-Sundrum case) or not. In the first case one can limit the solution of (jD.lj) 
to 9 > and continue the result symmetrically into 9 < 0. Note that imposing the 
symmetry 9 — ► —0 on the perturbations goes hand in hand with imposing 9 — > —0 for 
the source which means that one should add a corresponding delta-function 5(9 + 9') to 
the right hand side of (jD.lj) . The Israel-condition (j5.51j) on the brane with 9 = then 



gives -fc- 



0, as in the Randall-Sundrum case. 



Since in our setup we do not have any convincing argument for imposing the 9 — > —9 
symmetry also on the metric fluctuations (the scalar field), the second case, which allows 
a breaking of this symmetry will be the case of our choice. The Israel-condition (j5.51j) 
now merely indicates that the perturbation (the scalar field) should have a continuous 
first derivative at the brane. 

In order to obtain a uniquely defined solution we have to impose two boundary 
conditions on g ( - x '(9,9') one at 9 = n/2 the other at 9 = —n/2. It turns out that the 
requirement of square integrability of our solutions (with the correct weight-function) 
forces us to discard one of the two fundamental solutions at the boundaries 9 = ±7r/2, 
namely (^(A, 9). 

Before going into details, we introduce the following abbreviations: 



p(9) 



R 2 



S (9) 



a{9) 1 (9) 



Eq. (jD.lj) then becomes 

[p(9)g^'(9, 9')}' - (A 2 - 1) s(9)g^(9, 9') = 6(9 - 9') , 



(D.5) 



(D.6) 
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where ' denotes derivatives with respect to 9 here and in the following. We start by 
taking 0' > 0. The case 9' < is fully analogous to the one considered up to minus signs. 
Our ansatz for g( x \6,6') is 



A(6') V [ x \-6) 
g W(9,9') = { B(0>5 A) (0) + C(0>2 A) (0) 

D{e')<pf\e) 



-f <e< 



o 

B' 



< 
< 



0, 

< & 

< 5 . 



(D.7) 



We impose the following boundary and matching conditions on g^ x \ 
uniquely determine the coefficients A, B, C, D: 

1. continuity of gM(0, 6') at = 0. 

2. continuity of g w '(6, 9') at B = 0. 

3. continuity of #( A )(0, 0') at = 0'. 

4. jump condition of g {x) \e,6') at = 0'. 

Using the ansatz (jD.7|) we obtain from the above conditions: 

= 0, 
= 0, 

= o, 



, B') which will 



[A(6')-B(6')]^\0) - C(6>)^\0) 

[A(tf') <^ A)/ (0) + C(^)^2 A)/ (0) 

[B(B')-D(B'))cp[ x \B') + C(B')^(B') 

[B(&>)-D(B')}p[ X) '(&') + C(9')^ x) '(d') 



vie 1 ) 



the solution of which is easily found to be: 

R 2 



A(6') 
B(9') 

C(6') 

Dm 



2 ^S A) (0)^S A) '(0) ' 

r 2 ^\d') (0)^(0) + ^ A) (ovi A) '(o) 



-R 2 - 



(A) (A) „ 
<Pl ,<P2 >° 



^S A) (o)^r ;/ (o) 



w 

£(0") - # 2 



(A) (A) n 



w 



^ A) ,o 



(D.8) 



(D.9) 
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tp[ , </?2 j denotes the Wronskian of tp± and </4 at 6 1 = 0. In obtaining 



where W 
(ID. 91) we used the relation 



(D.10) 



R 2 p{6) 

As already mentioned, the case 9' < can be treated in complete analogy to the case 
9' > 0. Combining the two results gives the final expression for the solution g^ x \9, 9') of 
(JEH) (for A = 2,3,...): 

[ A{\9'\)^\\9\) -f <0sign(0')< 0, 

gM(6,tf) = I B(\9>\)<p[ x) (\9\) + C(\9>\)^\\9\) <0sign(0')< \9'\, (D.ll) 
[ D(|#'|)^ A) (|#|) <#sign(£')< f- 

It is easily verified that the Green's function (|D.11|) has the property 

gM(9,9') = gW(9',9), (D.12) 

as expected for a self-adjoint boundary value problem. We will not need the result (jD.llj) 
in its full generality. We focus on the case where sources are located at the brane (9' = 0) 
so that (jD.ll|) reduces to: 

R 2 A x \\o\) 



9 W (9,Q) 



(A = 2,3,...). 



2 ^ A) '(o) ' 

For the sake of completeness, we also give the solution g {1) (9,9') to flgH 



d 



R 2 a(9) 1 (9fd9 _ 



a{9) 1 {9) 



9') 

09 



6(9 - 9') 
a(#) 7 (#) 3 



(D.13) 



(D.14) 



The general procedure of finding the solution is fully analogous to the case of g^(9, 9') 
(for A = 2, 3, . . .), the only difference being that in regions where 9^9', ()D.14)1 reduces 
to an inhomogeneous differential equation. In the interval < 9 < ir/2 the general 
solution is given by 



ij(9)=M0) + AM9)+Bi; 2 (9), 



with 



MO) 
MO) 



R 2 



In 



2cutanh(^f ) 
1, M0)=2\n 



cosh 
tanh 



(D.15) 
(D.16) 



+ cotlr 
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With these definitions, it is easy to verify that g^\0, 9') is given by: 

f ^ 2 (0) -f < 0sign(0') < 

gW(6,6') = const.+M\0\)+MW\) + c{ < 0sign(0') < \6'\ , (D.17) 

I -02 ( 1 6>' |) <0sign(0') < f 

with 

c =-mra- (D ' 18) 

By W [ipi, ip2, 0] we again mean the Wronskian of ipi and tp 2 evaluated at 9 = 0. 
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Appendix E 



Detailed computation of the 
corrections to Newton's law 



In this appendix we give a detailed computation of the sum 

s[ s ,M = E Asin ( As )^o7V 0<e -l (K1) 

A=2 Pi (0) 

with 

y[ X) (0) _ z{0)*? [1 ~ <0)f 2F1 Pfc jrj A + 1; *(fl)] 
<//i A) (0) w(A - l)z(O)** 2 [1 - z(0)] 2 2 F, [*±1 ±±3; A + 1; *(0)] ' 

where we explicitly excluded 8 = since the Fourier sum (|E.1|) does not converge for 
this value of 0. 1 We also used the definition z{9) = tanh ho (~ — 8j] . 

Note that the sum S[s,9,u] contains the contribution to the two-point function 
coming from zero mode and higher Kaluza-Klein modes. Inspired by the Randall- 
Sundrum case, we try to separate the two contributions by functional relations be- 
tween contiguous Gauss hypergeometric functions. Using (15.2.15) of reference |192j 



1 Note that this mathematical delicacy about the divergence of the Green's function on the brane also 
applies to the Randall-Sundrum case of the Fourier integral representation of the two-point function. 
The obvious remedy is to keep 6 strictly greater than zero during the whole calculation and eventually 
take the limit 8 — ► 0. 
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for a = b — l = (A + l)/2, c = A + lwe obtain 



[1-2 
2 

ATI 



A+3 A+3 



i F 1 



2 2 
A+l A+3 



A + 1; z{8) 



(E.3) 



A+l; 2(g) 



A + l 



A + l A+l 



;A + l;z(g) 



In this way we are able to rewrite S[s, 8, u>] in the form: 



S[s, 8, uj] = Si [s, 6, u] + S 2 [s, 8, u] , 



(E.4) 



where 



S x [s,8,u) 



2 2(0) 1 - z(8) 
[l-2(0)f 



X 



A sin (As) 



5 2 [s,0,o;] 



A=2 
1 2(0) 



A 2 -l 



2F1 [^.^; A + l; z(g)] 
U(0)J a fi[a±l,*±5;A + l;*(0)] 



2(g) 



(E.5) 



1 - 2(g) 



W 2(0)3 [1 - z(0)f 

A sin (As) 2(6*) 



2F1 [ 



A+l A+l . 

2 ' 2 ' 



A + l; 2(g)] 



U(0)J 2^1 [ 



A+l A+3 . 
2 ' 2 ' 



A + l; 2(0)] 



(E.6) 



From the asymptotic formula (|F.1J) we infer that now the sum S\[s,8,u] is convergent 
even for 8 = due to the additional power of A in the denominator. We therefore obtain 



Si[s,0,a>] 



2 2(0)5 ~ A sin (As) 2 2(0)5 



^1_ 2(0 )^ A 



E 



col -2(0) 



7T 



1 



COS s 



sm s 



(E.7) 



Since we recognized in the last sum the two-point function of Einstein's static universe, we 
attribute the contribution coming from Si[s, 8, u] to the zero mode in the Kaluza-Klein 
spectrum. We also should mention that in S^s, g, u;] we still need 8 > for convergence 

A-l 

since only then the factor [2(g)/2(0)] 2 provides an exponential cutoff in A for the sum. 

The main challenge in the evaluation of the two-point function is therefore to tame 
the sum S2 [s, 8,uj}. Our strategy is the following: first we extend it from A = to 00 by 
adding and subtracting the A = and A = 1 terms. Next, we replace the sum over A by 
the sum of two integrals employing a variant of the Euler-Maclaurin sum rule called 
the Abel-Plana formula (see. e.g. [193J, p. 289-290). Finally, we shall see that the 
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resulting (exact) integral representation will allow us to extract the desired asymptotic 
of the two-point function on the brane in the distance regime of interest (i?Cr< Ru)- 
As announced, we start by extending the range of the sum from to oo. Since the 
addend with A = vanishes, we only need to subtract the A = 1 term with the result: 



r f „. 1 z(0) 1 - z{9) D| 
S 2 [s, 9, uj\ = — — - m [1 — z(9)\ sin s tttt^ 77vT^H s > oj\ 



2u z{9) 1 - z(0) 
where R[s,9,cu] is defined by 

A sin (As 



R[ s ,e,u] = J2 



A=0 



A + l 



z(e) 

z(0) 



u z (6)-2 1-^(0) 



(E.8) 



We used the identity 
"A + l A + 3 



;A + 1;*(0) 



[i- Z {o)r i 2 F 1 



A + l A 



A + l; z(0) 



(E.9) 



(E.10) 



in the last step (see 15.3.3 of [192]). The next step consists of extending the sin function in 
(|E.9J) to an exponential and taking the imaginary part out of the sum (as in the Randall- 
Sundrum case). Now we make use of the Abel-Plana formula, first considering the 
partial sums: 



R in) \s,9, 



-f(0, s, 9, cu) + -f(n, s, 6,u) + / /( A. .s. 9. } <IX 



+ i 



f(iy, s, 9, u)-f(n + ly, s, 9, u)-f(-iy, s, 9, u)+f(n - ty, s, 9, u) 



o 



e 2iry _ I 

where for the sake of clarity we introduced 

Xe lXs \z{9)~ X 



dy 



;e.ii) 



f(\,8,0,u) = 



A + l 



A+l A+l . 

2 > 2 ' 



X + l;z(9)] 



In the limit n — > oo we find: 

oo 

R[s,9,u] = $\ [ f(X,s,9,tu)dX + i 



U(0)J a F a [*±l,*=i;A + l;*(0)] ' 

oo 

f(iy,s,9,u) - f(-zy,s,9,uj) 



(E.12) 



dy 



(E.13) 



V 

Ei?i[s,e,w] 



sR 2 [a,e,u] 
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E.l Leading order asymptotic of ^ 0, uj}} 

We focus first on R\ [s, 9, uj] and observe that the function /(A, s, 9, to) is holomorphic (in 
A) in the first quadrant. 2 Therefore, in perfect analogy with the Randall-Sundrum case 
treated in appendix O we can use Cauchy's theorem to replace the integration over the 
positive real axis by an integration over the positive imaginary axis. 3 After substituting 
A with xy we obtain : 



OO ly 



' l + y 



|*(o) 



2 



2Fl [l+3i i= Hs. 1+J!/ . z(0) ] 



The upper deformation of the path of integration is advantageous for at least two reasons. 
Firstly, through the appearance of the exponential in the integrand, convergence on the 
brane does no longer rely on 9 > and we can set 9 equal to zero in (|E.14|) to obtain: 



RM = _ / vjl^u- ^[y-y i+ '^ ( °>w (E.i5) 







Secondly, the form ()E.14|) is perfectly suited for obtaining asymptotic expansions for 
large distances. Note that since s G [0, 7r], we carefully avoided saying "for large s" since 
this would not correspond to the case of our interest. We try to compute the two-point 
function in a distance regime R <C r <C Rjj, distances clearly far beyond the size of 
the fifth dimension but far below the size of the observable universe. Rewritten in the 
geodesic distance coordinate s this becomes uj/ sinh ( l rr) C s<C 1. It is this last relation 
which complicates the computation of the asymptotic considerably. We are forced to look 
at an asymptotic evaluation of (|E.15|) (or ()E.14jl ) for intermediate s values such that we 
cannot make direct use of Laplace's method. The correct way of extracting the above 
described asymptotic is to expand the ratio of hypergeometric functions in (jE.15|) in a 
power series of the small parameter 1 — z(0). Using rel. (15.3.10) and (15.3.11) (with 



2 Everything but the ratio of the hypergeometric functions is clearly holomorphic. The dependence of 
Gauss's hypergeometric function on the parameters is also holomorphic so the only danger comes from 
the zeros of 2-F1 [2^1 A + 1; z(0)l . These, however, turn out to be outside of the first quadrant. 

3 To be mathematically correct, one has to apply Cauchy's theorem to the quarter of the disk of 
radius R bounded by the positive real axis, the positive imaginary axis and the arc joining the points 
R and tR, see Fig. IC.l| in appendix |0 The above statement will then be correct if in the limit R — > 00 
the contribution from the arc tends to zero. We will verify this in section TE. 41 of this appendix. 
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m = 1) of |192| we find after some algebra 



2 F l [^,=^;l + ty;z(0)] 2 { y> 



2 7 + 2 ^(i^) + ln [1-2(0)] 



+ 



1 + zy 



2^(2) - 2^ ( ) - In [1 - 2(0)] 



1 + zy 



-1 + zy 



2^(1) - 2^ 



1 + zy 



x 



In [1 - 2(0)] - V(l) - ^(2) + ^ ( ) 



In [1 - 2(0)] 
1 + zy 



X 



[1-2 



+ O {In [1 - 2(0)] [l-2(0)] 2 } . 



(E.16) 



In the above relation 7 denotes Euler-M ascheroni's constant and ip(z) = T'[z]/T[z] 
the so-called Digamma-f unction. Note that we developed up to linear order in 1 — 2(0) 
since we also want to compute the next to leading order term in the asymptotic later in 
this appendix. 

We are only interested in the imaginary part of R\[s, 0, u] and it turns out that the 
first term in the expansion (JE.16J) after insertion in (|E.15J) can be integrated analytically: 



I 7 

s{4 0) [ S ,0,o;]} = 3M 



ye-*'® 



l + iy 



dy. (E.17) 



The last step is justified since also the real part of the first integral in (jE.17|) converges, 
as follows immediately from the asymptotic ip(z) ~ In 2 — j- + 0(-^). Using 



9 



1>( 



1 +ty, 



7T {Tiy 
— tanh — 
2 V 2 



(E.18) 



( P2| . p.259, 6.3.12) and introducing the /3-function (see [Ml, p.331, 3.311, 2.) 



(E.19) 
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we find 




.0 



By /3' in (|E.20|) we mean the derivative of (3 with respect to its argument. 



E.2 Leading order asymptotic of ^ {i?2[s, 0, uj}} 



We can treat B,2[s,9,uj] starting from (|E.13|I in very much the same way as Ri[s,9,u]. 
Noting that we can again put 9 equal to zero since the exponential factor e 2ny guarantees 
the convergence of the integral, i^fs, 0, a>] becomes explicitly: 



R 2 [s, 0,u] 



-y + w 
i + y 2 ' 

-y - w 2 
i + y 2 



-ys_ 



+iy l+w . 
2 ' 2 ' 



l + iy-z(0)] 



l F 1 [^,=^L;l + ty;z(0)] 



oV S . 



i-w i-iy . i 
2 ' 2 ' 



iy,z(0)] 



l—iy —1—iy. 
2 ' 2 ' 



1 - iy;z(0 



dy. (E.21) 



If we now employ the expansion (jE.16|) two times in (jE.21|) we find after some algebra: 



s{4 0) [ S) o,c]} 



9 



dy } (E.22) 
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and using (|K.18|) we have 

s{4 0) [ s ,o,^]} = 



7T y 



2 e 2n y - 1 



-y* 



tanh {^j + e ys tanh dy 



(E.23) 



-7T 



— tanh sinh(ys) 



-7T 



-7T- 





d 
ds 



o 



'1 + e^) 



sinh(ys) dy 



cosh(ys) 
(1 + e 7 ^) 2 



dy 



1 d 
2ds 



it 



1 d 



2 ds 2 - 



2^1 



2,2 ;3 

7T 



S 

; 

Ti 



+ 



1 



du 



u 







(l + u) : 



d-u 



2 + - 



s s 
2,2 + -;3 + -: 

7T 7T 



While we substituted y by it according to m = e in the third line, we used (3.194, 



p. 313) of |194j with u — 2, u = (3 = 1 and \x = 2^f s/n (valid for s < 2ir) in the last 
line. In order to simplify the hypergeometric functions, we first use another relation 
between contiguous functions, namely eq. (15.2.17) of |192j with a = 1, b = k, c = k + 1 
and z = —1 and then the formulae (15.1.21) and (15.1.23) of 192J together with the 
duplication formula for the T-function to obtain: 

i 2 F l [2, k-k + l- -1] = 2 F 1 [1, k; k; -1] - 2 F X [1, k; k + 1; -1] 



2- h V^ 



T[k] 



r [§] r p±i] 



(k-i)p(k) 



+ (i-k)p{k). 



(E.24) 



After making use of this in ()E.23j) we finally obtain for the contribution to lowest order 
in 1 - z(0): 

1 d 



9 



0,uj] = - — 



7T/ 7T 



1 + /3(2 + -) 

7T/ 7T 



Summarizing the results to lowest order in 1 — z(0), we therefore have 
9f{i4 0) [s,0 



2S^ 7T 



(E.25) 



(E.26) 



9 



3 /5(2 -^ )+ ( 1 + ^) /3(2 + ^ 

7T/ 7T V 7T/ 7T 
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Expanding (|K.26|) around s = we find: 

3 0,"]} = - ^ + O(s) , (E.27) 

3 {4 0) [ S ,0,.]}=(i-^) s + O( s =), 

with £ denoting Riemann's (^-function. As expected the contribution from [s, 0, cu] 
is sub- leading with respect to the one from r{ [s, 0,lj]. 



E.3 Next to leading order asymptotic of ^ {R\ [s, 0, lj] } 



We would now like to obtain the next to leading order term in the asymptotic of 
S {Ri[s, 6, u]}, that is the term obtained from (jE.15|) by taking into account the lin- 
ear contribution in 1 — z(0) of the expansion (jE.16|) . Since the calculation of the in- 
tegral obtained in this way is rather cumbersome and anyway we do not need the full 
s-dependence, we do not intend to evaluate it fully. Extracting the leading s-divergence 
at s = is sufficient for our purposes here. To start, we insert the linear term of (|E.16Jl 
in (lE~T5j) : 



9f{j# ) [a,0 ) a;]} 



(E.28) 




1 + iy 



In [1-2(0)] -V(l) -V>(2) + V 



1+ty 



3 + iy 



2^(2) - 2ip 
2 

+ ip 



3 + ty 
2 

-In [1 

1 + iy 



-In [1-2(0)] 



[1 - z(0)] 



After some algebra and by using the functional relation of the Digamma-function ip(z) 
(see jliEp 



ip (2 + 1) = ip (2) + 



we find 



Of \B$\s,U, id] 



2(0) 



ye- ys ^{...}dy 



(E.29) 



(E.30) 
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with 



9f{...} = 3 



+ 



2 

4 

x 



2^(2) - 2^ ( ) - In [1 - ^(0)] 



2^ (1) - 2ij) 



2 

1 + n/ 



In [1 - z 



x 



+ 



1-iy 



In [1 -^(0)] - ip (1) - ip (2) + 2ip 
1 + iy 



1 + iy 



2ip (1) - 2*/j 



- In [1 - 2(0)] 



(E.31) 



We will now keep only the terms proportional to y 3 in (jE.30|) (y 2 in (|E.31(l ) since only 
these will contribute to the leading 1/s 4 singular behavior. After a few lines of algebra 
we find 



1 - 2(0) 



(E.32) 



[l-2 7 - In [l--z(O)]] / y s e- ys Q 



l + iy 



dy 



+ [1-2(0)] / y 3 e- ys % ' 1 1+W 
o 

+ terms involving lower powers of y. 



ft 



if, 



1 + iy 



dy 



The first integral can be reduced to the integral (jK.17|) simply by replacing each power 
of y by a derivative with respect to s and by taking the derivatives out of the integral: 



y 3 e- ys % 




(E.33) 



We could not find an analytic expression for the second integral in (jE.32|) . However, all 
we need is the first term of its small s asymptotic 4 and this can be easily obtained from 



4 In the sense s -C 1. 
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the large y asymptotic of 3? [if) ((1 + ly) /2)]. Since asymptotically 

if) (z) ~ in z — + C( — ) , [z — > oo in | argz| < 7r) 



In z + 0(-r 

2z 12z 2 z 4 ' 



(E.34) 



we expect that this term will contribute logarithmic terms in s and we find after straight- 
forward expansions 



3? 



if) 



1 + iy 



ln y --— 2 +0\^ 

2 6y 2 



(E.35) 



From the last formula (jE.35|) we learn that within our current approximation of keeping 
only highest (cubic) power terms in y, it is sufficient to take into account only the 
contribution coming from ln|. Therefore, we need to calculate the following integral 



J yV'tanh^) \nydy 



7T d 3 

2d? L 



-ys 



tanh 



Try 
2 



liaydy 



(E.36) 



o 

s — 



where we used (|E.18|) once more. Note that the term proportional to In 2 can be ac- 
counted for by adding a contribution of the type of the first integral in (jE.32|) . One 
way to find the asymptotic of 3~{s) for small s is to integrate the following asymptotic 
expansion of tanh x 



tanhx = 1 - 2e~ 2x + 2e~ 4x - 2e~ 6x 



+ . . . = 1 + 2 £(-l)"e 



v — Ivx 



(E.37) 



which converges for all x > 0. Inserting this in the definition of J{s) we find: 



J{s) 



Iny 



l + 2^(-ire 



u=l 



dy 



(E.38) 



7 + Ins | 2 y^ ^ u+1 -f + \n(s + un 

e ^— ' S + V1X 



7 + In s 2 



+ -(7 + lnvr)/3(- + l)+-^(-l) 



v + 



where we used the series expansion of the /3-function given e.g. in [T94] . The last two 
terms in the above formula are finite in the limit s — > and therefore will play no role 
in the asymptotic (since they are multiplied by a factor of 1 — z(0)). 
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We are now in a position to assemble all contributions to the small s asymptotic of 



{r?[s,0,u>}}: 



3{itf [ S ,0,o;]} 
1 - z(0) 



(E.39) 



[1 - 2 7 - In [1 - z 



, 7T d 3 

2 In 2 — — 

J 2rfs 3 



1 

-- + -P - 

S 7T \7T 



[1 - *(0)] 



7T ^ 

2rfs 3 



7 + Ins 2 



+ - (7 + ln7r)/3 (- + 1 

7T \7T 



| 2 y, (-ir 1 ln(z,+ f) 

i/=l i" 



ttI -z{0) 



8 - 61n2 - 6 In 



In s N 



where we used the following series expansions for the /3-function and for the last term in 
flE~38j) : 5 



\7T/ S 

+ =ln2 + 0(s), 

^ (-l)^ +1 ln(z/ + ^) 1 2 
E 77Ts = o [ ln ( 2 )] - 7 ln(2) + 0(a) 



i/=i 



(E.40) 
(E.41) 

(E.42) 



Note that only the term /3(s/n) gave a contribution to the leading 1/s 4 divergence. 
Relations (jE.26|) and (jE.39|) constitute the main result of this appendix. They provide 
the desired asymptotic behavior of the two-point function in the regime of intermediate 
distances u / sinh < s < 1. A couple of remarks are in place. Clearly (|E.26|) is 
not only an asymptotic result, but is valid for all s G [0,7r]. Secondly, we were writing 
0(\ns/s 3 ) in (|E.39|) to indicate that we dropped all contributions coming from terms in 
the integrand (jE.31|) lower than cubic order. Finally we point out that the development 
of the ratio (fFTTfijl in power of 1 — z(0) indeed generates an asymptotic in the correct 
distance regime. This is clear from the fact that each new power in l — z(0) is paired with 
an additional power of y 2 in the integrand of 3 {Ri[s, 6, 00]}. This by itself implies upon 

5 At this point we would like to thank Elmar Tcufl for the evaluation of the leading order term of the 
sum in IE. 421). 
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multiplication by e ys and integration over y an additional power of 1/s 2 . We therefore 
effectively generate an expansion in powers of y/l - z(Q)/s =1/ [scosh(^f)]. It is also 
easy to understand why this rough way of counting powers works. The reason is that 
the asymptotic properties of our integrals 3 |i?[ n ^[s, 0, to] j are determined only by the 

behavior of the corresponding integrands for large values of the integration variable y. 
The terms 9ft [ip ((1 + ly) /2)] and 3 [ip ((1 + xy) /2)] do not interfere with the above power 
counting since their behavior for large y is either logarithmic in y (for 9ft) or constant in 
y (for 3). 

Finally we want to point out that the expansion in powers of 1/ [scosh (^)l is 
not useful for obtaining information about the Green's function for distances s <C 
ioj cosh(^). In this limit the asymptotic expansion of the hypergeometric functions 
in (jE.2|) for large values of A proves the most efficient way to recover the properties of a 
Green's function in a 4-dimensional space. See appendix |F] for the detailed form of the 
Green's function at distances smaller than the extra dimension. 



E.4 Estimate of the arc contribution to Ri[s,6,uj] 

This final subsection is dedicated to the verification that the arc denoted C2 in Fig. IG.ll of 
appendix IH1 gives a vanishing contribution to the integral Ri[s,9,u>] and thus justifying 
the representation (|E.15[) . We need to know the large A asymptotic behavior of the 
hypergeometric functions entering the definition of Rx[s,9,u] in (jE.12|) . The relevant 
formula ()F.1|) can be found in appendix |F] where we discuss the short distance behavior 
of the Green's function. After canceling all common factors from the ratio of the two 
hypergeometric functions we find: 



2^1 [ 



A+l A+l . 

2 ' 2 ' 



X+l;z{9)] 



aft P±i,*fi;A+l;*(0 



(E.43) 



A + l 



-ZJ9Y 



A+l 
2 



-1/(0) 



-1/(0) 



A+l 

2 



[1 + e 



V 1 



-1/(9)1 



[ 1 + e -K0)]2 [! 



-(o)] 



l + O 



'1 



where e v ^> is defined in (|F.2|I of appendix IF1 (see also (|F.5|) and ()F.6jl ). 
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Substituting A = Re ltp , we now estimate Ri[s, 9, uj\. 



Xe lXs \z(9) 
<h X + l |i(0) 



2 F 1 [^,^;X + l;z(9)} 
2 F 1 [Ml,^i;A + l;z(0 



(E.44) 



< 



C-2 



Ae 



< c 



A + l 

zjoy 1 



[ml 2 F l [^,^;X+l;z(0)] 
[l + e -"W]-3 [l- e -"W]- 



r/A 



-"(«) 



'(0)] 



3 

Ml 



-i/(0)] 



A' 



^Rse 1 ^ I 



Re*? - 1 

-_R(ssin<£— In a cos ip) 



2 COS </3 i ] 



rdy? < / e 



-its (sin p+ — — cos <p] 



dip. 



where C is a constant independent of R and the estimate of the ratio of the two hyper- 
geometric functions is valid for large enough R. We also introduced a as in appendix IF1 
by (|F.7[) . Note that a < 1 for < 9. In the last line we estimated the inverse of the 
square root by an arbitrary constant (e.g. 2), certainly good for large enough R and for 
all ip. All what is left is to observe that 



r In a 1 In a 1 
< e = mm|l, J < sin (p H cos ip . 



s > 0, a 1 > 1 



(E.45) 



to obtain the desired result: 



...dX 



C-2 



2 

< K'R J , 



-Rse 



dip 



71 



K'Re~ Rse ^ for R 



oo. 



;e.46) 
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Appendix F 

The two-point function at very short 
distances (r <C R) 



The aim of this appendix is to calculate the behavior of the two-point function at dis- 
tances smaller than the extra dimension r <C R. We will do so by using an appropriate 
asymptotic expansion for large parameter values of the hypergeometric functions in (|E.2|) . 
Such an asymptotic can be found for example in j!95j or |196j . The general formula is 1 



2-Fi [a + n, a — c + 1 + n;a — 6+1 + 2n; z] 



2 a+b T[a - b + 1 + 2n] 2 
T[a-c + l +n]T[c-b + n] 



-v(n+a) 
v(n+a) 



c—a—b- 



1 + I - 

n 



with v defined via the relation 



2(1 



Z) 2 



;f.2) 



We need to calculate 



S[s, 9, u] = A sin (As" 



(A) 



(0) 



X=2 



;f.3) 



1 Since we noticed that the formulae given by Luke, p. 452 eq. (20)-(23) |196j and the formula given 
by Erdelyi |195| differ by a factor of 2 a+b we performed a numerical check. Its outcome clearly favored 
Erdelyi's formula which we reproduced in IF. II 
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with the ratio (p[ X \9) / (p'[ X \o) given in (jE.2|) . Using the above expansion (jF.l|) two times 
we find after numerous cancellations 



S[s, 9, uj] ~ — - 



4 (l-z{9)\ 2 [l + e-"W]-' [l-e-"W] 



£ 

A=2 



1 - *(0) 

A sin (As 
A- 1 



^(0)2 e 2 



[l + e -K0)]a[l_ e -i/(0)]-5 2^) 



"(0) 

e 2 



x 



-1/(9) 



;f.4) 



where we used the following abbreviations 



2(0) = tanh 2 



7T 



-1/(0) 



tanh 



2 



2 V2 



(F.5) 
(F.6) 



Note that (jF.6|) is obtained after some lines of algebra by inserting (jF.5|) in (jF.2j) . Still 
the sum in (jF.4|) diverges for 6 = but it turns out that it can be performed analytically 
for 9 > so that after the summation the limit 9 — > exists. For the sake of a lighter 
notation we introduce the symbol a by 



-"(0) 



-"(0) 



tanh [| (j - 9)] 
tanh (^) 



(F.7) 



It is now straightforward to evaluate the sum over A: 



E A sin (As) j , x *^ sin (As) , 



A=2 



A=2 



A=2 



, sin s + sin(As) a A + 



sin [(A + l)s\ A+1 



A=l 



A=l 



A 



■ — a sin s + ■ 



1 



sm s 



2 l( a + z)- coss 2 
+ a cos s arctan 



sin s In (l — 2a cos s + a 2 ) 



a sins 



1 — a cos s 



(F.8) 
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The full result can therefore be written as 



S[s, 9,lj] ~ — - 



4 (1 - z(9)\ 2 [l + e~ v ^] ~ 2 [l-e- v ^] 5 z(0)§ e 



3 3./(fl) 
2 



u\l-z(0)J [i + e -Ko]l[i_ e -(o)]-f e -f 



x 



x a sm s < 



1 1 1 

2a | (a + i) - coss 2 

a sin s 

+ cot s arctan 



In (l — 2a cos s + a 2 ) 



a cos s 



;f.9) 



Although the last result is given in closed form, we have to emphasize that its validity 
is restricted to distances smaller than the size of the extra dimension R. 
It is now save to take the limit 9 — > in (IF.9I) . The result is 



lim S\s, 9, ui] 



sinh(^f) 



sm s 



CO 



4sin 2 (|) 



In 



2 sin ( | 



IT — S 



cot s — 1 > . (F.10) 



We note that it is the first term in the curly brackets on the right hand side of (jF.9j) 
which is responsible for reproducing the characteristic short distance singularity of the 
two-point function. The main result of this appendix is therefore: 

sinh ( U7r S \ ( 1 <7r ^ 

\imS[s,9,u] 5lA2 J_ + - + 0(s]ns) } . (F.ll) 

e^o oj { s 2 J 

The 5-dimensional short distance singularity can even be obtained including the extra 
dimension. Apart from the logarithmic term in the curly brackets of (jF.9|) . all but the 
first one are finite at short distances (in 9 and s). Developing the denominator of this 
term, we recover the usual euclidean metric in IR 4 : 



M 1 

a H OSS ~2 



sinh 2 (if) 

where the dots denote higher terms in 9 and s. 



9 2 + s 2 



+ 



(F.12) 
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